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Abstract 

A review is given on the theory of vortex-glass phases in impure 
type-II superconductors in an external field. We begin with a brief 
discussion of the effects of thermal fluctuations on the spontaneously 
broken U{1) and translation symmetries, on the global phase diagram 
and on the critical behaviour. Introducing disorder we restrict our- 
selves to the experimentally most relevant case of weak uncorrelated 
randomness which is known to destroy the long-ranged translational 
order of the Abrikosov lattice in three dimensions. Elucidating possi- 
ble residual glassy ordered phases, we distinguish between positional 
and phase- coherent vortex glasses. The study of the behaviour of iso- 
lated vortex lines and their generalization - directed elastic manifolds 
- in a random potential introduces further important concepts for the 
characterization of glasses. The discussion of elastic vortex glasses, 
i.e., topologically ordered dislocation-free positional glasses in two and 
three dimensions occupy the main part of our review. In particular, 
in three dimensions there exists an elastic vortex-glass phase which 
still shows quasi-long-range translational order: the 'Bragg glass'. It 
is shown that this phase is stable with respect to the formation of 
dislocations for intermediate fields. Preliminary results suggest that 
the Bragg-glass phase may not show phase-coherent vortex-glass or- 
der. The latter is expected to occur in systems with weak disorder 
only in higher dimensions (or for strong disorder, as the example of 
unscreened gauge glasses shows). We further demonstrate that the 
linear resistivity vanishes in the vortex-glass phase. The vortex-glass 
transition is studied in detail for a superconducting film in a parallel 
field. Finally, we review some recent developments concerning driven 
vortex-line lattices moving in a random environment. 
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1 Introduction 



Since its discovery by Kammerlingh Onnes in 1911, superconductivity has 
attracted generations of physicists. In 1935 Fritz and Heinz London (see e.g. 
London (1950)) developed a very successful phenomenological theory which 
describes both the perfect conductivity as well as the perfect diamagnetism 
of superconductors. As discussed later by London (1950) this theory can 
be motivated by considering superconductivity as a phenomenon character- 
ized by long-range order of momentum p. Bohr-Sommerfeld quantization 
/ Psds — nh on a torus gives fluxoid quantization (London 1950). Ginzburg 
and Landau (1950) combined London's electrodynamics of a superconductor 
with Landau's theory of phase transitions, creating a powerful phenomenolog- 
ical description of superconductivity. The transition to the superconducting 
phase corresponds here to the breaking of the U(l) symmetry of the com- 
plex order parameter ^ and the appearance of off-diagonal long-range order 
(ODLRO). 

In a pioneering work Abrikosov (1957) showed the existence of a second 
type of superconductors which (for sufficiently strong external magnetic field) 
allows for a penetration of quantized magnetic flux in the form of vortex lines, 
which form a triangular lattice, reducing the perfect diamagnetism and cre- 
ating a source for dissipation due to the motion of the vortex-line core driven 
by the transport current. In this case the continuous translational symmetry 
of the system is broken in addition to the U{1) symmetry. Effects from ther- 
mal fluctuations, although studied already since 1960, were considered to be 
extremely small because of the large correlation length and low transition 
temperatures of conventional superconductors (Ginzburg 1961). 

To keep the superconducting properties, vortices have to be prevented 
from moving by pinning centers. An early theory of pinning of isolated vor- 
tex lines (Anderson and Kim 1964) shows the absence of dissipation only at 
zero temperature. Thermally activated hopping leads to a small but finite 
dissipation at low temperatures (as compared to the height of energy barri- 
ers). Larkin (1970) extended this theory to the Abrikosov vortex-line lattice, 
showing the destruction of its translational long-range order. Although - 
in principle - the Abrikosov phase could thence be considered not to really 
differ from a pinned vortex- line liquid (and hence from the normal phase), 
the generic phase diagram of conventional superconductors was assumed to 
practically be that of Abrikosov (1957) with a now finite correlation length 
of the vortex-line array. 

In the 1980s this picture was changed by two initially unrelated devel- 
opments: the discovery of high-Tg superconductors by Bednorz and Miiller 
(1986) and the much better understanding of random system, in particular 
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of spin glasses and of random-field systems (for a recent review see Young 
(1998)). In the high-Tc superconductors with their elevated transition tem- 
peratures and pronounced anisotropy, fluctuation effects became now very 
important as can be seen for instance from the observed melting of the 
vortex-line lattice (Cubitt et al. 1993, Zeldov et al. 1995). Moreover, for pure 
systems it was demonstrated (Moore 1989, Moore 1992, Ikeda et al. 1992) 
that thermal fluctuations - which prohibit true long-range translational or- 
der of the vortex-line lattice (VLL) only in < 2 dimensions - destroy the 
ODLRO of the gauge invariant order parameter in the Abrikosov phase even 
in higher dimensions. Thus, in three dimensions thermal fluctuations restore 
the U{1) symmetry of the Ginzburg-Landau Hamiltonian but nevertheless 
allow for the existence of a vortex lattice. This finding is paralleled by an 
earlier observation of Schafroth (1955) that an external magnetic field above 
a critical strength destroys Bose-Einstcin condensation of an ideal Bose gas 
which still shows some remanent diamagnetic moment. 

For systems with disorder the idea emerged that despite of the destruc- 
tion of true translational long-range order the system could show a phase 
with some kind of glassy long-range order, the "vortex glass" (Fisher 1989). 
Because of the residual rigidity in the vortex-line array, arbitrarily large en- 
ergy barriers now exist, leading to a highly non- linear resistivity (Feigel'man 
and Vinokur 1990, Fisher 1989, Feigel'man et al. 1989, Nattermann 1990, 
Fisher et al. 1991a) 

p(j) ~ e-^^*/^)" (1.1) 

where denotes an exponent < // < 1 and jt a threshold current. Since 
the linear resistivity vanishes, the system is truly superconducting. 

In the following years vortex pinning and dcpinning as well as flux creep 
under the action of an external current was investigated by many researchers 
to a great extent. A brilliant summary of the results of these efforts till 
1994 is given in the extensive review article by Blatter et al. (1994) (see also 
Brandt (1995), Gammel et al. (1998), Giamarchi and Le Doussal (1998)). 
It is not the intention of the present article to provide an updated version 
of these reviews by discussing the results obtained since then. Instead, we 
want to focus here mainly on one particular aspect of the theory, namely on 
the discussion of the equilibrium phase diagram of weakly disordered type-II 
superconductors in an external magnetic field. We want to demonstrate that 
the notion of a 'vortex glass' is not a blurred expression for the hopelessly 
intricate situation in a disordered system (as some physicists may still claim), 
but that it has a well defined meaning. The knowledge of the equilibrium 
properties is also important for the proper understanding of situations close 
to equilibrium, e.g. for the discussion of flux creep under the influence of a 



5 



small external current. 

Despite of the conclusion common to all authors mentioned above of ex- 
pecting a non-analytic current-density dependence of the resistivity in the 
glassy phase it seems to be indicated to refer here also to the differences 
between these approaches. Fisher (1989) and Fisher et al. (1991a) in their 
definition of the glassy phase started from correlation functions measuring 
ODLRO and focused on the long-range (glassy) order of phases (see chapter 
2.4). On the other hand, Feigel'man et al. (1989), Nattermann (1990) and 
subsequently Korshunov (1993) and Giamarchi and Le Doussal (1994) con- 
sidered primarily the glassy order of the vortex-line array, i.e., they focused 
on the positions of the vortex lines. In both cases, the expression 'vortex 
glass' was used. The above mentioned differences in the lower critical di- 
mensions for the breaking of the U{1) and of the translational symmetry in 
pure systems however suggests, that there may be also different lower crit- 
ical dimensions for a phase- coherent and a positional vortex glass. In this 
review we concentrate mainly on the positional glass. If, in particular, the 
positional vortex glass is free of large dislocation loops such that an elastic 
description of the vortex array is possible, the positional vortex glass will be 
called elastic vortex glass. The most prominent example of an elastic vortex 
glass is its three-dimensional version: the so-called 'Bragg-glass' (for a recent 
brief review see Giamarchi and Le Doussal (1998)). Whether non-elastic 
vortex-glass phases exist is still unclear. 

Systems with columnar disorder which leads to the formation of the so- 
called 'Bose-glass' phases (see e.g. Tauber and Nelson (1997)) will be com- 
pletely neglected in this review since their physics is substantially different 
from that for systems with point disorder, to which we restrict ourselves here. 
However, we also include the discussion of various phases found for driven 
systems far from equilibrium, which share some features with the equilibrium 
phase diagram. 

The article is organized as follows. In chapter 2 we present a brief sum- 
mary of the Ginzburg-Landau theory of type-II superconductors and a short 
discussion of the influence of thermal fluctuations and of the effect of point 
disorder in the critical region. We also define the different types of vortex- 
glass order and give a brief account of results obtained for models with strong 
disorder - the so-called gauge glasses. In chapter 3 we review the behaviour 
of a single vortex line and its generalizations - D-dimensional directed man- 
ifolds - in a random potential. This simple, but not at all trivial system 
allows for a discussion of different aspects of glassiness of a system. Chap- 
ter 4 is devoted to the superconducting film in a parallel field, a geometry 
which allows for a very detailed description of the vortex glass phase as well 
as of the transition to the normal phase both for the static and dynamic 
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quantities. In chapter 5 we discuss an impure superconducting film in a field 
perpendicular to the film plane. It turns out that dislocations destroy the 
positional vortex-glass phase in this geometry. The 'Bragg-glass' phase of 
a bulk superconductor as well as its stability with respect to dislocations is 
considered in chapter 6. A short account of recent activities on driven vortex 
lattices in impure superconductors is presented in chapter 7. We close the 
paper with a brief summary of the results of this article (chapter 8). The 
appendix contains some technicalities and a list of recurrent symbols. 



7 



2 Ginzburg-Landau description 



In this chapter wc give a very brief introduction into the mean-field theory 
and the effects arising from thermal and disorder fluctuations in type-II su- 
perconductors in the framework of the Ginzburg-Landau theory. Since there 
is extensive (and partially contradicting) literature on thermal effects it is im- 
possible to include all related references. However, we attempt to include the 
most recent articles on the subject which may serve as more comprehensive 
guides to further references. 



2.1 The Ginzburg-Landau model 

In 1950 Ginzburg and Landau proposed a phenomenological description of 
superconductors by introducing a two-component order parameter ^'(r) = 
|^^i.^|g«</'(r) which couples in a gauge- invariant form to the magnetic field 
described by the vector potential A(r) (Ginzburg and Landau 1950). The 
density ns(r) of superconducting charge carriers (i.e., of the Cooper pairs), 
which is a central quantity of the earlier London theory (London and London 
1935), is related to ^'(r) by ns(r) = |^'(r)p. The Ginzburg-Landau (GL) free 
energy is given by 

+-\(tV-^A)^ +±(VAA-H)2j, (2.1) 

where $o = hc/2e denotes the fiux quantum, a{T) oc (T — Tco), Tco is the 
mean-field transition temperature, H is the external field, and m denotes 
the mass of a Cooper pair. The GL free energy is characterized by two basic 
length scales, the coherence length ^ and the penetration depth A, which are 
related to the parameters of Hgl by 

^^(^' = SWF- 

Here |^oP = \c(\/P denotes the saturation value of j^^p in a homogeneous 
current free state for T < T^o and H = 0. For our further discussion it is 
convenient to use the following rescaling to introduce dimensionless quantities 
A' and r' 

* = A = ||A', r = r'e. (2.3) 
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This leads to 



+ |(zV' - A')*f + K^(V' A A' - U/H^'^yy (2.4) 
Here we introduced the Ginzburg number 

the Ginzburg-Landau parameter k = \/^ and the reduced temperature Tcq = 
(Tco — T)/Tco. i!^^^ = $o/27r,^^ is the mean-field upper critical field and 
— 47r/3|\E'o|^ is the thermodynamic critical field. 



2.2 Mean-field theory 

Within mean- field (MF) theory, the GL free energy has to be minimized with 
respect to the fields \l'(r) and A(r). The resulting GL-cquations 

(iV- A')'*' + A*'+I*'r *' = (2.6) 

\a\ 

and (for a < 0) 

k'^V A(V AA') + A' = -V> (2.7) 

then have to be solved with the appropriate boundary conditions. As is clear 
from equation (2.4), the only two parameters which will enter the solution 
in the bulk, are the GL-parameter which plays the role of an inverse 
effective charge of the ^ field, and the strength of the external magnetic field 
H' = YL/H^^. 

For K < l/\/2 (type-I superconductors), mean-field theory yields for 
T < Tco and H < Hc{T) a phase with vanishing resistance and perfect 
diamagnetism. The transition to the normal phase at Hc{T) is first order. 

For K > I/V2 (type-11 superconductors), on the other hand, perfect 
diamagnetism exists only up to the field H^l^ pa (i7^^/2«;^)(ln«;-F0.08). For 
larger fields, magnetic flux penetrates the sample in the form of quantized 
vortex lines, each carrying a flux quantum $o- The energy per unit length 
of the vortex line is therefore given by e\ — ($o/4vr)ifci — SqIhk with £0 — 
(<l>o/4vrA)^ as the important energy scale per unit length. 

The vortex lines form a triangular 'Abrikosov' lattice (Abrikosov 1957, 
Kleiner et al. 1964) of spacing = {2/^/3y/^a, where a = {^/Bf/'^. The 
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Abrikosov lattice (or 'mixed') phase shows both broken translational symme- 
try and off-diagonal long-range order (ODLRO), i.e., broken U(l) symmetry 
of the order parameter. Both broken symmetries vanish simultaneously if 
H reaches . One should however take into account that the correlation 
function for ODLRO (see chapter 2.3.2) (\I'*(r) ^l'(r')) - even if calculated in 
MF approximation - shows strong spatial variations due to the rapid change 
of the phase. Indeed, for a system of radius R the tangential phase gra- 
dient on the boundary is of the order R/2a^ which corresponds to a phase 
change of 27r on a distance ^ Ana^/R (with a ^ lOOnm and R ~ 1cm, 
Ifj) ~ lO^^nm which is smaller than an atom, see Brandt (1974)). Thus 
(\l/*(r) ^{r')) cannot be very meaningful as a physical observable. Quanti- 
ties with a physical significance should be in particular gauge invariant. In 
the GL-description these are the amplitude |^| of the order parameter, the 
'super- velocity' V0 — -^A = V0 and the magnetic induction B = V A A. 
All other physical quantities can be expressed in these fields, e.g. the current 
density can be written as 

j = -^|*rV0. (2.8) 
m 

In treating the vortex system, two main approximations have been used: 
the lowest Landau level (LLL) approximation, which is valid sufficiently close 
to {H > |if^^), and the London approximation, which is valid at 

intermediate and small fields B < 0.2H^^, where ^ -C a. The precise range 
of applicability of the LLL is still under debate (see e.g. O'Neill and Moore 
(1993), Li and Rosenstein (1999) and references therein). 

In the London approximation one neglects amplitude inhomogeneities, 
|\E'| = 1^0 1 5 which leads to a diverging energy density at the vortex cores. 
The position rj(s) of these cores is parameterized by the label i of a vortex 
and the variable s along the contour of the vortex lines. The Ginzburg- 
Landau Hamiltonian takes then the form 

^London = ^ / d'v { o\\V + ^ Af + ^{V A A - uA . (2.9) 

This functional has to be regularized near the vortex cores, e.g. by excluding 
tubes of radius ^ around the vortex cores from the volume integration. The 
phase ct)(r) of the complex order parameter is a multivalued function since 
changes by 27i along a path surrounding a vortex line. We decompose 
now into a vortex part (p^ and a "spin- wave" part 0sw The vortex part is 
assumed to fulfill the saddle point equation apart from the position of the 
vortices rj(s). With the London gauge V • A = this yields 

V'0v(r) = O, r^r,{s) (2.10) 
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and 

V A ( V0v) = 27rm(r) . (2.11) 
Here m(r) denotes the vortex-density field 

m(r) = J^m^l ds^5('\r,{s) - r) , (2.12) 

i 

where the integration is along the vortex line i which carries the vorticity 
mj = ±1. If the spin- wave part 0sw(r) vanishes on the surface of the sample, 
(psw and (f)y decouple. Since the vector potential A appears only quadratically 
in i^GL it can be integrated out by using the saddle-point equation which is 
the second GL-equation in the phase-only approximation 

A + A^V A (V A A) = -^^V0v (2.13) 

Taking the curl of (2.13) gives the modified London equation 

A^V^B(r) - B(r) = $om(r), (2.14) 

where B = V A A can now completely be expressed in terms of the vortex 
degrees of freedom given by the vortex density field m(r), equation (2.12). 
The London Hamiltonian then takes the form 

^London = ^ / d'^r | - 1 *o ( V0sw)' + ^(V A B)^ + -i(B - Uf 
1 J ym 47r 47r 

(2.15) 

In most parts of this review we will use the London picture, since it remains 
valid in the vortex phases we will describe, in particular in the elastic glass 
phases. 

The elasticity theory of the Abrikosov lattice for an isotropic supercon- 
ductor was worked out by Brandt (1977a, 1977b). The distortion of a vortex 
line is described by a two-component displacement field u(X, 2;), where the 
lattice vector X = X„^m = ((2n + m)aA/2, mv^aA/2) denotes the rest posi- 
tion of the vortex line in the plane perpendicular to H = ffz. In many cases 
one can go over to the continuum description: u(X, z) — > u(x, z) = u(r). On 
large scales L ^ A the elastic energy of the vortex-line lattice is then given 

by 

^el = lj d^xd''-h {cn(Vx ■ U)2 + C66(Vx A U)2 + C44(V||U)2} , (2.16) 

where Cn ^ C44 B^/47r. The Abrikosov phase is characterized in particular 
by a non-zero shear modulus cee ~ j^^^ ^1 — "^mf j , which vanishes both 
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at Hc2 and i?ci, reaching a maximum of cqq ~ ifci-f^c2/567r in between. One 
should however take into account that the elastic free energy is in general 
non-local, which is expressed in a strong dispersion of Cn and C44 on scales 
smaller than A (see, e.g., Brandt (1991)). 

In the absence of pinning centres, the system in the Abrikosov phase 
behaves superconducting only for currents parallel to the magnetic field. For 
currents with components perpendicular to the field the Lorentz force drives 
the vortex-line array, which leads to metallic behaviour with resistivity p ^ 
Pi^B / (Bardeen and Stephen 1965). Here pn is the resistivity of the normal 
phase. 

Before we come to the discussion of fluctuation effects, we want to con- 
sider a possible extension of the model (2.1) which describes an isotropic 
superconductor. High-Tc superconductors, however, are characterized by 
a pronounced layer structure, which results in an inhomogeneity and in a 
strong spatial anisotropy of the effective mass of the Cooper pair such that 
m is replaced hy M — m/e^ for electrons moving perpendicular to the lay- 
ers. Typical values for e are eyBCO ~ 0.16 and eBscco ~ 10~^ for the two 
high-Tc materials YBCO and BSCCO. It was shown by Blatter et al. (1992) 
(for a more detailed discussion, see Blatter et al. (1994)), that in the case 
of K 3> 1 the result for a thermodynamic quantity Q{i9, H,T,^, X,e, A) of 
an anisotropic superconductor (where denotes the angle between the mag- 
netic field direction and the xy plane, and A the strength of disorder) can 
be obtained from the corresponding result Q for the isotropic system by the 
relation 

Q{^, H, T, e. A, e. A) = sq Q{e^H, T/e, A, A/e), (2.17) 

where e| = e^cos^-i? + sin^-i?, sy = Se = St = ^ for volume, energy and 
temperature, and sb = sh = l/e?? for magnetic fields. Since T and A are 
increased by a factor 1/e with respect to the isotropic system, it is clear 
that fluctuation effects, which will be considered in the following sections, 
are drastically enlarged (by a factor up to 100) in these materials. 

If the spatial anisotropy is so large that the coherence length = 
in z direction becomes of the order of the layer spacing s, the discreteness 
of the layer structure becomes relevant. In this case, new effects such as 
a decoupling of the layers [for B > B2B ~ (e/s)^ ln(s/^)] may occur. The 
appropriate description is then the Lawrence- Doniach model (Lawrence and 
Doniach 1971). We will not attempt to cover in this review also these par- 
ticular features of strongly layered materials, instead we restrict ourselves in 
the following to the discussion of the isotropic superconductor, knowing that 
the results for the anisotropic case can be found from the relation (2.17). 
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The neglect of layer-effects is also supported by the following argument: 
Our theoretical analysis will be mainly based on the elastic description of 

the vortex lattice and our main interest concerns features on large length 
scales. Sufficiently weak disorder will indeed effectively couple to the vortex 
lattice only on very large length scales, where the elasticity of the lattice may 
be described by local elasticity theory. Since within the London approach 
all information about anisotropy and even about the layered structure is 
encoded in the dispersion of the elastic constants, we expect that the large- 
scale properties are independent of these details (provided that one is in the 
parameter regime where the London approach is valid and that disorder is 
sufficiently weak). 



2.3 Thermal fluctuations 

So far we have ignored the influence of fluctuations, i.e., of conflgurations 
which do not fulfill the GL equations. These can be taken into account if we 
interpret the GL free energy as an effective Hamiltonian from which the true 
free energy J-{T, H) has to be calculated as 

JP(T,H) = -Tin (^J P^DAe~^«L/^^ . (2.18) 

Since the only material independent common feature of type-II superconduc- 
tors is flux quantization it is natural to build from $o ^-nd T a characteristic 
length scale (Fisher et al. 1991a) 

At = ~ 2 • lO^-p^, (2.19) 

^ IGn^T r[K]' ^ ' 

which is the same for all materials. Since the energy per unit length of a 
vortex line (and hence also its stiffness constant e\\, see below) is of the order 
Eq = ($o/47rA)^, At denotes the length scale on which the mean squared 
displacement of a vortex line is of order A. Since At is so large, thermal 
fluctuation effects are expected to be small (however, see our remark about 
strongly anisotropic systems in the previous chapter 2.2). In d — 3 dimen- 
sions the Ginzburg number can be expressed as Gi ~ (kA(0)/At)^. As can 
be seen directly from (2.4) and (2.18) the contribution from fluctuations in 
^(r) and A(r) will indeed be small, if both Gi <C l^col'^"'^ and 1. 



2.3.1 Zero external field 

In zero external field, H = 0, to begin with, it was shown by Halperin et al. 
(1974) that for type-I superconductors, where k, < 1/ \pl and typically Gi ^ 1 
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(e.g., Gi ^ 10"^^ for aluminium), fluctuations in the vector potential A(r) 
render the transition first order. 

For typc-II superconductors, on the other hand, the situation is less clear. 
It has been argued that the transition remains second order (Hclfrich and 
Miiller 1980, Dasgupta and Halperin 1981). In the high-Tc compounds with 
large values for k (kybco ~ 100, kbscco ~ 60) and large Ginzburg numbers 
(GivBCO ~ 10~^, GiBscco ~ 1) fluctuations in the vector potential are weak 
compared to those of the order parameter. Then there exist two critical 
regions. In the outer critical region 

|Te| < Gi^/^^-*^), |Te| > (Gi/«^) ^(^-'i) , (2.20) 

where Tc = (Tc — T)/Tc now denotes the reduced temperature with respect to 
the true transition temperature Tc, fluctuations of the order parameter lead 
to an Xy-hke critical behaviour. Fluctuations in the vector potential can 
be neglected in this regime. Since the coherence length ^{T) ps ^(0)|rc|~^ 
with u = uxY ~ I in = 3 dimensions increases now more strongly than 
the penetration depth A ^ Ao|tc|~^"^''''/^, the effective value of k ~ ^-(4-'i)/2 
decreases until both lengths are of the same order. This signals a cross-over 
to a second critical regime with (probably) inverted XY behaviour (Helfrich 
and Miiller 1980, Dasgupta and Halperin 1981). In this asymptotic regime 
A and ^ scale in the same way with the correlation exponent Uxy (Olsson 
and Teitel 1999). It should be mentioned, however, that other scenarios have 
been proposed (for a recent discussion of earlier results see Kiometzis et al. 
(1995), Radzihovsky (1995a), Herbut and Tesanovic (1996), Herbut (1997), 
Folk and Holovatch (1999), Nguyen and Sudb0 (1999)). 

In d = 2 dimensions fluctuations prevent the formation of a long-range 
ordered phase. As was shown by Pearl (1964, 1965), the effective London 
penetration depth Lg = 2}? /s diverges with decreasing s, where s denotes the 
fllm thickness. Therefore fluctuations in the vector potential can be neglected 
and the system in zero external field shows a Kosterlitz-Thouless transition to 
a quasi-long-range ordered phase (Doniach and Hubermann 1979, Halperin 
and Nelson 1979). 

2.3.2 Finite external field 

Next we consider the case of finite external field. The most obvious effect of 
thermal fluctuations on the vortex-line lattice is melting (Eilenberger 1967, 
Nelson 1988). Melting has been seen experimentally in YBCO (Safar et al. 
1992, Kwok et al. 1992, Charalambous et al. 1993, Kwok et al. 1994, Liang 
et al. 1996, SchiUing et al. 1996, Welp et al. 1996) and BSCCO (Pastoriza 
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et al. 1993, Zeldov et al. 1995, Hanaguri et al. 1996). To estimate the 
melting temperature one may use the phenomenological Lindemann criterion 



(2.21) 



where u denotes the displacement of a vortex line from its rest position, 
(. . .) the thermal average and Cu ~ 0.1... 0.2 is the Lindemann number. 
Since the shear modulus Cqq vanishes at and H^^, melting will occur 
by approaching both critical fields. Close to Hd the melting line H^i{T) is 
roughly given by (Fisher et al. 1991a) 



The region Hd < H < i^mi, where the vortex lines form a liquid, is ex- 
tremely small, except for the vicinity of Tc{H = 0) where A diverges. H^i 
is reduced with respect to due to fluctuations (Nelson 1988, Nelson 

and Seung 1989). We note, however, that for a proper calculation of the 
melting curves the dispersion of the elastic constants has to be taken into 
account. In anisotropic and layered superconductors, Blatter and Gcshken- 
bein (1996), following an earlier suggestion by Brandt et al. (1996), found 
an additional fluctuation-induced attractive van-der-Waals interaction be- 
tween vortex lines, which may lead at very low temperatures to a first order 
transition between the Meissner and the Abrikosov lattice phase. 

At large fields the melting line H-^2{T) is reached if (Brandt 1989, Hough- 
ton et al. 1989) 



For ifmi < H < H^2 the vortex lines form a solid (cf. figure 1). 

Alternatively, one could start directly from the GL-Hamiltonian and 
study the fluctuation corrections in the vicinity of (see, e.g., Lee and 
Shenoy (1972), Bray (1974), Thouless (1975), Ruggeri and Thouless (1976), 
Ruggeri (1979), Brezin et al. (1985), Affleck and Brezin (1985), Brezin et al. 
(1990a), Radzihovsky (1995b)). As first observed by Lee and Shenoy (1972), 
the fluctuations in a d-dimensional superconductor in an external fleld are 
like those of a (d — 2) -dimensional system in zero fleld, suggesting that the 
mean-field phase transition to an Abrikosov phase with ODLRO is destroyed 
by fluctuations in dimensions d < rf^^^^*-* = 4 (the upper critical dimension 
of the system is now 

^ODLRO ^ rpj^.g soj^g^J^at surprising conclusion is 
in agreement with calculations of Moore (1989, 1992) and others (Glazman 




(2.22) 




(2.23) 
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Figure 1: Schematic illustration of the phase diagram of type-II super- 
conductors in the (T, H) plane. In mean-field theory, vortices form the 
Abrikosov vortex lattice between the critical fields and . Due 
to thermal fluctuations the vortex lattice is solid only between the melt- 
ing fields ilmi and ilin2 (dark shaded area) along which it melts into a 
vortex liquid. 

and Koshelev (1991a, 1991b), Ikeda et al. (1992)), who found the destruc- 
tion of ODLRO by strong phase fluctuations below a lower critical dimension, 
d < 0?°°'"^°, starting from the existence of a periodic solution, i.e. of a vortex 
lattice within the GL-theory. 

Detailed considerations show that the lower critical dimension 0?^°^^*^ = 
4 for systems with screening and d^^^^'^ = 3 for systems without screening 
(Moore 1989, 1992). Since the aforementioned calculations on fluctuation ef- 
fects close to use the LLL approximation which usually neglects screen- 
ing, it is clear that a simple dimensionality shift by 2 does not work here. 
Qualitatively, it is plausible that screening as an additional source of fluctu- 
ations increases d^^^^*"*. A similar effect is also observed for gauge glasses 
(see section 2.4). Quantitatively, the shift of d^^^^'-* can be traced back to 
the strong dispersion of the tilt modulus for k > (Moore 1992). (Note 
that the compression modulus does not affect the formation of ODLRO.) 
Breakdown of dimensionality reduction by 2 was also shown by Radzihovsky 
and Balents (1996) in a layered superconductor in a parallel fleld, where 
^d^^^° — 2.5 but the upper critical dimension 0?^°^^° = 5, and the vortex 
lattice is still stable in d > 2 dimensions. 

ODLRO is conventionally deflned by a non-vanishing limit of the gauge 
invariant pair correlation function 

C2(ri,r2;r) = (^(ri)^*(r2)e^(''^/*°)-''r*-^) (2.24) 
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for |ri — — > oo. ((•) denotes the average over thermal fluctuations.) 
Note that C2(ri, F) itself depends on the path F between ri and r2 along 
which the vector potential is integrated. The proposal of Moore (1992) to 
keep only the longitudinal component of A to make C2 path independent 
but preserve its gauge invariance corresponds in the London gauge to the 
complete neglection of the phase factor e*(^'^/*o) ir ^ in equation (2.24). 

In fact, a non- vanishing asymptotic expression for C2(ri,r2;F) if |ri — 
r2| ^ 00 may not be the appropriate definition for the existence of long- 
range order in cases in which topological defects arc forced into the system 
by external boundary conditions or (as in type-II superconductors) external 
fields. The most simple counter example is an Ising magnet below Tc with 
anti-periodic boundary conditions, which force a domain wall into the sys- 
tem. Wall fluctuations will then suppress magnetic correlations. The loss 
of ODLRO due to thermal fluctuations in type-11 superconductors - if con- 
cluded from the asymptotic behaviour of C2 - is related to the fact that phase 
fluctuations 50(r) of the order parameter are related to (shear) distortions 
u(r) of the vortex-line lattice via (Moore 1992, Ikeda et al. 1992, O'Neill and 
Moore 1993, Baym 1995) 

50(k) = g![iMi/(k)- V.(k)_ ^2.25) 



a 



Here we used the Fourier transforms (50(k) and u(k) of the flelds 6(p{r) — 
0(r) — (f)o{r) and u(r) and k = {kx, ky^lcw). (j)o{r) denotes the phase in the 
ground state. Note that equation (2.25) is not in contradiction to equation 
(2.10) since V^0v(r) 7^ only at the vortex sites. Equation (2.25) implies 
that on large length scales phase fluctuations are amplifled displacement 
fluctuations of the vortex-hne lattice. Destruction of ODLRO below d — A 
dimensions does not mean however - this is the point of view we will take here 
- the absence of any ordered phase. This can be concluded from the London 
approximation (2.15): Expressing B(r) with the help of (2.14) in terms of the 
vortex degrees of freedom, the only phase fluctuations left are spin-wave like 
which lead to a destruction of conventional long-range order only in d < 2 
dimensions. The energy of the vortex lattice can be approximated at low T by 
an elastic Hamiltonian (2.16) which may be supplemented by a contribution 
from dislocations. For d > 3 such a system shows a phase with translational 
long-range order (TLRO) for H-^aiiT) < H < H-^taiT). (In d = 2 dimensions, 
elastic fluctuations reduce the order to quasi-TLRO characterized by a power- 
law decay of correlations.) In the elastic description, which we will use in 
most parts of this review, the appropriate order parameter for TLRO is 

^/;Q(r) =e^*^-"W (2.26) 
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with (■0Q(r)) 7^ in the Abrikosov phase. Q denotes a reciprocal lattice 
vector of the Abrikosov lattice. From this point of view the loss of ODLRO 
has no physical significance. C2 is therefore not the appropriate quantity to 
define order. This conclusion is in agreement with a number of numerical 
investigations in which a clear indication for a transition to an ordered phase 
was seen in two and three dimensions. Hu and MacDonald (1993) and Kato 
and Nagaosa (1993a) have considered the pair correlation function of the 
superfluid density 

C4(ri,r2) = (|^(rOn^(r2)P). (2.27) 

Using the LLL approximation both groups found a clear indication for a 
first order melting transition in 0? = 2. Their work was extended by Sasik 
and Stroud (1993, 1994a) to three-dimensional systems. They considered the 
helicity modulus 

= ^'^(r - v')5^^ - ^Oa(r)iMr')>c, (2.28) 

A'=0 

where A' is an additional vector potential and j the current density. In the 
LLL (where A — > 00) they found a rapid increase of 7^^ at the freezing tem- 
perature, indicating the formation of a vortex-line lattice. Sasik et al. (1995) 
have subsequently shown that the vortex liquid-to-solid transition is not ac- 
companied by a divergence of the correlation length for phase coherence. In 
their simulation C2 decays exponentially even in the solid phase in agreement 
with the predictions of Moore (1992), Ikeda et al. (1992), and Baym (1995). 

It has to be mentioned, however, that the point of view we adopt in this 
article - namely that the loss of ODLRO docs not rule out the existence of a 
vortex lattice - is not shared by all authors. In particular Moore (1989, 1992, 
1997) argues that there is no mixed phase in type-II superconductors at finite 
temperatures. The observed effects in the behaviour of resistance and mag- 
netization are explained as cross-over phenomena which would disappear in 
the thermodynamic limit. Moreover, the vortex lattices found in the Monte- 
Carlo simulations mentioned before is considered to be an artifact produced 
by the quasi-periodic boundary conditions used in these simulations. Moore 
and co-workers (O'Neill and Moore 1992, O'Neill and Moore 1993, Lee and 
Moore 1994, Dodgson and Moore 1997, Kienappel and Moore 1999, Moore 
and Perez-Garrido 1999) have tried in their own Monte-Carlo simulations to 
avoid these effects (which they consider to be crucial) by placing the two- 
dimensional superconductor on the surface of a sphere. In their studies no 
freezing transition to a vortex-lattice state was observed. It should however 



7a/3(r,r') 



d{3a{v)) 
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be mentioned that zero-energy modes connected with the rigid rotations of 
the film on the sphere and dischnations arising from topological constraints 
for a triangular lattice on the sphere may obscure the transition. 

To conchidc: Although the absence of ODLRO at finite temperatures 
in the mixed phase of type-11 superconductors (see, e.g., Houghton et al. 
1990) and its consequences for the existence of a vortex lattice are still under 
debate, the most simple scenario is the existence of a vortex lattice in the 
absence of ODLRO below d = A dimensions. The continuous translational 
invariance of the system is reduced to finite translations by lattice vectors X 
of the Abrikosov lattice. In contrast to ODLRO, the lower critical dimension 
for TLRO in pure systems is = 2. 

Since the critical regime is enlarged in a non-zero external field (in d — ?>) 
to (Ikeda et al. 1989) 



with respect to the zero-field critical region, but is still too small to describe 
the difference between H^{T) and the melting line, Feigelman et al. (1993) 
have argued that between the Abrikosov lattice and the normal phase there 
may be an intermediate liquid phase which still shows longitudinal super- 
conductivity. We will not follow this idea here since more recent extensive 
simulations show only a single transition between the Abrikosov and the 
normal phase (Hu et al. 1997, Hu et al. 1998, Nguyen and Sudb0 1998, 
Nordborg and Blatter 1997, Nordborg and Blatter 1998, Olsson and Tei- 
tel 1999), although further proposals for an intermediate phase were made 
recently (Tesanovic 1999, Nguyen and Sudb0 1999). 

In the rest of this article we will ignore the existence of these exotic 
phases in pure superconductors (although their existence cannot be ruled 
out completely), but concentrate instead on the influence of randomly dis- 
tributed frozen impurities on the vortex-line lattice phase. The impurities 
will be assumed to be completely uncorrelated as already mentioned in the 
introduction. Columnar or planar defects lead to very different physics (for 
a recent review on the so-called Rose glass, see Tauber and Nelson (1997) 
and references therein) and will not be discussed in this paper. 

2.4 The influence of disorder 

Disorder can be introduced in the Ginzburg-Landau model (2.1) by assuming 
small random contributions to all parameters which characterize the system, 
i.e. a, P, m and H. We will restrict ourselves here to the case of systems 




(2.29) 



19 



with random mean-field transition temperature Tco, i.e. we substitute 



TeO^TeO + 5Teo(r) 



(2.30) 



with 



<5Tco(r) 5Teo(r')=e'5€(r-r) 6T^, 



(2.31) 



where the overbar denotes the disorder average and S^{r) is a (5-function of 
widths ^. 

In mean-field theory, a randomness in Tco would correspond to a smear- 
ing out of the transition. However, as first shown by Khmelnitskii (1975), 
thermal fluctuations permit the occurrence of a sharp phase transition. Ac- 
cording to the so-called Harris criterion (Harris 1974) weak randomness does 
not change the critical behaviour if the exponent of the specific heat of the 
pure system ctpure is negative, or, equivalently, if the correlation length expo- 
nent i/pure obeys 1/ > 2/d. Since for the XY-m.o<\e\ in three dimensions a is 
negative (Lipa et al. 1996), the zero-field critical behaviour for type-II super- 
conductors should be unchanged in the outer critical region. This would also 
apply in case that behaviour of the pure system in the inner critical region is 
also of (inverted) XY-typc. On the contrary, for type-I superconductors it 
was shown by Boyanovsky and Cardy (1982) that a sufficiently large amount 
of disorder may convert the first order transition into a second order one. 

If an external field is applied, the situation is quite different. Because 
of the problems connected with ODLRO even in pure systems, which were 
discussed in the previous chapter, it seems to be expedient to first consider 
the influence of disorder on the structural properties of the mixed phase. 

Inside the Abrikosov phase, disorder leads to a destruction of transla- 
tional long-range order as first shown by Larkin (1970). This follows from 
the fact that a randomness in the local value of Tc(r) leads to a random 
potential acting on the vortices. The specific form of the resulting coupling 
of the disorder to the vortex displacements will be given in the following 
chapters. The disorder averaged order parameter (e*Q ") for TLRO vanishes 
since disorder leads to diverging displacement u in the limit of large system 
sizes. However, it will be shown below that the correlation function 



may still obey an algebraic decay, which shows up in Bragg peaks in the 
structure factor (Giamarchi and Le Doussal 1994, 1995, Emig et al. 1999). 
Moreover, in analogy with spin glass theory (see Binder and Young (1986)) 
one may consider the positional glass correlation function 



5(Q,r) = (eiQ-[u(r)-u(o)]^ 



(2.32) 



2 



5pg(Q, r) = (eiQ-["W-"{o)]) 



(2.33) 
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as another signature of the existence of some residual, 'glassy' order of the 
Abrikosov lattice. Below we will call a system a positional vortex glass if 
5'pG(Q,r) decays not faster than a power law for |r| — > oo. For T ^ 0, 
S'pG(Q,r) approaches one. At non-zero temperature, (2.33) measures the 
strength of thermal fluctuations of the vortex lines around the ground state. 
Two hmiting cases seem to be conceivable: If the disorder acts effectively as 
a random force on the vortex lines, then the thermal fluctuations around the 
disordered ground state are identical to those of the pure system. In this case 
(2.33) is non-zero above d = 2 dimensions. In the opposite case of very strong 
pinning the vortex lines can be considered to fluctuate thermally in a kind 
of narrow parabolic potential and (2.33) may be finite even in d < 2. The 
relevance of these structural correlation functions for the glassy behaviour of 
the mixed phase will be further discussed in the following chapters. 

A complementary discussion of the glassy behaviour was proposed by 
M. P. A. Fisher (1989) and Fisher et al. (1991a). These authors use the 
correlation function for ODLRO as the starting point. Clearly, because of 
the relation (2.25), (^(r)^'*(0)) will vanish for r — > oo in d < 4 dimensions. 
In analogy to spin glasses they therefore introduce the correlation function 
for phase- coherent vortex-glass order (using the London gauge V • A = 0): 

CvG(r) = |(M/*(r)il/(0))^ (2.34) 

It is instructive to consider (2.34) in the London limit where CvG(r) 

l^oPl (e*^'^"^*-'"''^'^"^'-^-") ^0(r) describes the fluctuations around the ground 
state pattern 4>o{y). If the disorder is of random-force type, according to 
(2.25) thermal fluctuations should destroy C\g{y) below d — A. On the other 
hand for thermal fluctuations in a parabolic potential around the ground 
state, CvG(r) should be flnite. 

Dorsey et al. (1992) and Ikeda (1996a, 1996b) have calculated the vortex 
glass susceptibility xvg = / '^'^^ C'vg(i") a GL-Model with random transi- 
tion temperature in the LLL approximation. Dorsey et al. (1992) found in a 
mean-field calculation a second-order transition with a diverging vortex-glass 
susceptibility by approaching the vortex-glass transition temperature Tg(if), 
which is slightly below Tc2{H). In mean-field theory, the vortex-glass corre- 
lation length (^vG diverges with a mean-field exponent z/yj = 1/2. From this 
one can expect the existence of a non-zero limit of Cvg(i") for |r| — > oo below 
Tg{H). Taking critical fiuctuations into account, which can be done within 
a d = 6 — e expansion, the model is found to be in the some universality 
class as the Ising spin glass. In particular, the dynamical critical exponent 
is found to be ^ = 2(2 — r/), where rj ~ — e/6. It should be observed that a 
similarity between the Ising spin glass and vortex glasses was also mentioned 
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for so-called gauge-glass models (see Huse and Seung (1990)). On the other 
hand, the results of Dorsey et al. (1992) and Ikeda (1996a, 1996b) - if ap- 
plied to d = 3 dimensions - cannot be so easily reconciled with the findings 
from the clastic description of a vortex lattice in a random potential. As we 
will discuss in chapter 6, the vortex lattice exhibits indeed a glassy phase 
(the Bragg glass) for 2 < d < 4, in which the correlation function Cvg(i'), if 
calculated to lowest order in e = 4 — ci, vanishes for large |r| exponentially. 
The reason for this decay consists in the strong thermal fluctuations of the 
phases around the (distorted) ground state. To order e, these phase fluc- 
tuations are only weakly suppressed by disorder and hence Cyq{t) decays 
exponentially. However, higher order terms in e may still change this result. 
In principle, there could be two glassy phases, which show a non-vanishing 
S'pG(Q,r) and Cyc, respectively. For the moment we consider it to be more 
likely that there is only one glassy phase and the discrepancy between the 
results follows from the use of different approximations valid in o? = 6 — e and 
d — 4 — e dimensions, respectively. 

So far we assumed that the disorder is weak. On the other hand gauge- 
glass like models were proposed for the description of granular supercon- 
ductors or systems with strong disorder (Ebner and Stroud 1985, John and 
Lubensky 1986). Each superconducting grain of centre position is de- 
scribed by the phase (pi of the order parameter which is assumed to be con- 
stant within a grain. The Hamiltonian then reads (see e.g. Wengel and 
Young (1996)) 

n = -jY, cos {<j)i - <j)j - Aj - X^'aij) + ^ ^(V A a)=^ , (2.35) 

<ij> □ 

where Uij — J a(r) dr. ^ is the sum of all nearest neighbors of a cubic 

Vj <ij> 

lattice and is the sum over all plaquettes. a denotes the fluctuations of 

□ 

the vector potential which are limited by the bare screening length Aq. The 
influence of the external fleld as well as the contribution from randomness are 
assumed to be included in the Aij which arc taken to be independent random 
variables with a distribution between and 27i. A detailed discussion of the 
relation between the vortex glass (the expression is here understood in the 
sense that it describes the glassy phase of an impure type-II superconductor 
in an external fleld) and gauge glasses is given in Blatter et al. (1994). The 
model (2.35) is in particular isotropic in contrast to the GL-Hamiltonian 
which shows a pronounced anisotropy due to the presence of the external 
fleld H. In addition, the disorder of the gauge glass has a nature which is 
completely different from the one in equation (2.30). The flrst one couples 
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to the vortices via the phase of the order parameter, while the latter one 
couples only via the amplitude. As a consequence, the gauge-glass disorder 
distorts the vortices much more than local variations of T^o- 

For Aq oo, i.e. in the absence of screening, the gauge glass was investi- 
gated numerically by a number of authors (see Fisher et al. (1991b), Gingras 
(1992), Cieplak et al. (1991), Reger et al. (1991), Cieplak et al. (1992), 
Hyman et al. (1995), Maucourt and Grempel (1998), Kosterlitz and Akino 

(1998) , Olson and Young (1999)). While in two dimensions a transition to 
a glass phase is found to be only at T = 0, this transition takes place at 
finite temperatures in three dimensions. Huse and Seung (1990) found at 
the transition a diverging gauge-glass susceptibility 

XGG = J]CGG(ri-r,), (2.36) 

3 

where Cqg, in analogy with (2.34), denotes the gauge glass correlation func- 
tion 

^^GG(r.-r,) = |(e*(^-^.))|2. (2.37) 

The divergence of Xgg may signal the transition to a phase with non-zero 
limit of CGG(r) for |r| oo. If one assumes (2.34) as the definition of vortex 
glass order, then the gauge glass would be a vortex glass. 

More recently, the case of a finite Aq was considered (Bokil and Young 
1995, Wengel and Young 1996, Wengel and Young 1997, Kisker and Rieger 
1998). It turns out that screening seems to destroy the gauge glass transition 
in three dimensions. More recently Kawamura (1999) and Pfeiffer and Rieger 

(1999) have attempted to include the effect of anisotropy into the gauge glass 
model by assuming an extra contribution to A^j arising from the external 
field. However, also in this case no finite temperature glass transition was 
found. 
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3 Directed elastic manifolds in a random po- 
tential 



In the course of this article we will restrict our consideration of fluctuations 
to vortices and we will ignore other independent fluctuations such as those of 
the condensate amplitude or of the magnetic field. Vortices will be treated 
within the London picture, where vortex lines are represented as string-like 
objects. For low temperatures and weak disorder the vortex lines will fluc- 
tuate only weakly around the ground state of the vortex-line lattice (VLL), 
the Abrikosov lattice, where all lines are directed along the magnetic field. 

Before we analyze the VLL as an ensemble of vortex lines in its full com- 
plexity it is instructive to study a single vortex line. In general, a vortex can 
be characterized by two dimensions which depend on the physical realization 
under consideration: the 'internal' dimension D of the vortex considered as a 
'directed manifold', and the number of its displacement components. For 
example, a vortex line in a bulk superconductor has {D, N) = (1, 2), a single 
vortex line in a superconducting film has {D,N) = (1, 1), and a point vortex 
in a film corresponds to {D,N) — (0,2). This concept of elastic manifolds 
also applies to other physical systems such as interfaces between magnetic 
domains, for which {D,N) = (2,1). In all these examples the spatial di- 
mension of the system is d = D + N, since displacements are possible only 
in directions orthogonal to the z direction(s) along which the manifold is 
spanned. 

It is worthwhile to mention at this point that the analysis of random 
manifolds applies not only to single vortex lines, but to a certain extent (i.e., 
within a regime of length scales) also to vortex line lattices. For example, 
the VLL in a weakly disordered bulk superconductor resembles over a large 
range of length scales an elastic manifold with (D, N) — (3, 2). In this case, 
where the elastic manifold is spanned in all spatial dimensions including 
the displacement directions, d = D. Nevertheless, there is a fundamental 
difference between manifolds and VLLs, which is crucial for the physics on 
large scales: VLLs have a periodic structure in contrast to manifolds. 

To parameterize vortex conformations, we use the D-dimensional coordi- 
nate z along the direction of the magnetic field and a A^-dimensional coordi- 
nate X = x(z) in transverse directions, see figure 2. 

In this chapter we are mainly interested in qualitative aspects of vortex 
fiuctuations in the presence of disorder. To this end we describe the elastic 
energy of a vortex line, which arises from the kinetic energy of the super- 
currents and the magnetic field energy, in a harmonic approximation (i.e., 
we keep only terms of second order in the displacement) and we ignore non- 
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Figure 2: Illustration of the manifold coordinate system. The manifold 
conformation (solid line) is described by the A'^-component vector x as 
a function of the £)-dimensional coordinate z. The displacement u(z) = 
x(z) — X is defined with respect to a fiat reference position X (dashed 
line). 



localities or possible anisotropies in the elastic stiffness constant £||. Then 
the elastic energy can be written as 

= y"d^z |(Viix)^ (3.1) 

Therein the gradient Vy = (^, . . . , acts along the longitudinal z direc- 
tions. 

Pinning of vortex lines due to the presence of impurities, grain boundaries 
etc. can be described by a potential ^(x, z) which is 'quenched', i.e., frozen 
on the relevant time scale for vortex fluctuations. It yields a contribution 

Hpin = J d^'z V^(x(z),z) (3.2) 

to the total energy H = Ticx + Tipin of a vortex line. For simplicity, we will 
discuss here only point-like disorder. We will assume that this potential is 
Gaussian distributed with zero average and variance 

y(x,z)l^(x',z') = A(x-x')5(z-z'). (3.3) 

For simplicity we will restrict ourselves to disorder with short-ranged and 
isotropic correlations. Then the correlator A(x) = A(|x|) can be character- 
ized by its integral 

Ao = A(0) = / d^x A(x) (3.4) 
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and a correlation length. Such a pinning potential can arise from point-like 
impurities that locally suppress the density of the superconducting conden- 
sate. The correlation length of the disorder then also coincides with the 
superconducting coherence length ^. In principle, one can retain a finite cor- 
relation length in the x and z directions. However, it turns out a posteriori 
that - as long as the correlation length does not exceed the smallest scale of 
vortex conformations, which also is ^ - correlations in z direction are irrel- 
evant. Subsequently, we will often write A(x) = Ao5g(x), where denotes 
a Dirac delta-function smeared out on a scale ^. For semi-quantitative pur- 
poses we occasionally use ^^(x) = (27r^^)~'^/^ exp(— x^/2{^). If we denote the 
pinning force of an individual impurity by /pin, then 

Ao «^ /5nr^S'^e^^+^ (3.5) 

where n['^^ denotes the concentration of the impurities in the d dimensional 
space. 

Directed manifolds in disorder represent a paradigm for systems which 
are dominated by disorder. Since disorder leads to a frustrating competition 
between elastic and pinning energies, the structural order can be reduced 
substantially. In addition, the dynamics of the system can get extremely 
slow due to the quenched nature of disorder. Therefore such a system can be 
called a 'glass'. Because of the simplicity of the manifold model, it plays a 
paradigmatic role for glassy systems and it allows the identification and the 
understanding of many characteristic features of a more complicated 'vortex 
glass'. 

3.1 Equilibrium properties 

As we have seen, vortices in impure superconductors are a realization of 
'elastic manifolds' in 'random media' which are paradigms for disordered 
systems. We now proceed to summarize briefly some key properties of the 
latter model system. For more detailed presentations about these random 
manifolds the reader is referred to review articles (Kardar 1994, Halpin- 
Healy and Zhang 1995, Lassig 1998a), in particular to Blatter et al. (1994) 
for a discussion in the context of vortex systems, as well as to Nattermann 
and Rujan (1989), Belanger and Young (1991) and Young (1998) in the 
context of disordered magnets. Subsequently we will focus our attention on 
the question, in what physical quantities 'glassiness' appears in these systems. 

From a principal point of view, it is important to note that disorder 
breaks symmetries of the Hamiltonian of the pure system: translation in- 
variance x(z) — > x(z) -|- Xq for a constant shift Xq of vortex lines, as well as 
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an analogous rotation symmetry x(z) — > Tt ■ x(z), for a rotation matrix 7?,. 
Therefore it is interesting to examine to what extent the physical state of the 
manifold in disorder actually reflects these broken symmetries, or whether 
these symmetries are restored due to thermal fluctuations. 



3.1.1 Structure 

One quantity of primary interest is the structure of the manifold in disor- 
der which can be described in terms of displacement correlations. In the 
absence of disorder the vortex line has a flat ground state, x(z) = X for all 
z. Introducing the displacement u(z) = x(z) — X, shape fluctuations can 
be described by the relative displacement at points separated by a distance 
z — z' parallel to the magnetic field: 



W{z-z') = ([u(z) - u(z')]2) ~ |z - zf ^. (3.6) 

As already introduced above, (. . .) and 777 denote the thermal and disorder 
average respectively For large distances this quantity typically follows a 
power law with the roughness exponent ( — C{D, N). The manifold is called 
flat, if W is finite for |z — z'| — > cxd (then the convergence of W to its 
asymptotic value can be described by an exponent ( < 0), whereas it is 
called rough if W diverges (i.e., C ^ 0) where C = typically corresponds to 
l^(z) ~ In" z with some power a). 

In the absence of disorder thermal fluctuations are described by an expo- 
nent 

Cth = (3.7) 

The presence of disorder may increase the roughness and lead to a larger 
exponent C > Cth- In particular, disorder always induces roughness in dimen- 
sions 2 < D < 4, as we will discuss in more detail below. 

For D < 2N/ (2 + N) disorder is irrelevant at sufficiently high tempera- 
tures. There is a phase transition (see e.g. Imbrie and Spencer (1988) for 
the case D — 1) between a low-temperature phase, where the manifold is 
disorder dominated and essentially shows the structure of the ground state, 
and a high- temperature phase, where the manifold is entropically driven out 
of the ground state and shows a structure as in the absence of disorder. 

Actually, the physical situation can be more complicated if an upper criti- 
cal dimension A^cu exists, above which the low-temperature phase is governed 
by the Gaussian exponents. The existence of an upper critical dimension is 
still controversial. Moore et al. (1995) argue for N^u = 4 and Lassig and 
Kinzelbach (1997, 1998) argue for Acu < 4 in D = 1. We will not further 
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discuss this possible complication here, assuming A^cu to be large enough such 
that the vortex systems of physical interest are not concerned. 

For D > 2N / (2 + A^) one can think of the manifold as having a unique 
disorder-dominated ground state. For all temperatures entropic effects are 
too weak to detach the manifold from its ground state. 

3.1.1a Structural order parameter 

In order to have a tool to quantify the structural order of the manifold, 
we introduce 

M^) = e*''-"(^) (3.8) 

as order parameter. If the manifold performs only weak thermal fluctuations 
around its ground state, (■?/'k(z)) 7^ and disorder actually breaks the trans- 
lation symmetry of the manifold. Otherwise, if thermal fluctuations detach 
the manifold from its ground state, (■0k(z)) = 0. The correlation function of 
this order parameter, 

S{k, z - z') = {r^{z)M^')) ^ exp (^-ik • W(z - z') • , (3.9) 

is related to the displacement correlation 

W^0{z-z') = {[u^{z) - u^{z')][u0{z) ~ up{W)]). (3.10) 

The function W, previously introduced in (3.6), is simply the trace of the ma- 
trix Wai3- The approximate relation in (3.9) neglects higher cumulants of the 
displacement distribution and holds in general only for small displacement 
fluctuations. It even holds for large displacement fluctuations and actually 
is an identity, if the fluctuations of u have Gaussian distribution. 

3.1.1b Perturbative analysis 

A first qualitative insight into the relevance of disorder can be obtained 
from an elementary perturbative analysis. Such an analysis was performed 
originally by Larkin (1970) for vortex lattices and by Efetov and Larkin 
(1977) for the closely related charge-density waves. 

In this approach the manifold (at T = 0) is considered in an abso- 
lutely flat reference state x(z) = X, for which the pinning force FP''^(z) = 
— V_Ly(x(z), z) is calculated. We denote by V_l = ■ ■ ■ , g~) gradients 
in X directions in contrast to for z directions. From this force the disor- 
der induced manifold displacement is obtained in linear response theory as 
Ua{ci) = / d^z e~^'^'''Ua{z) — Gap{(l)Fp^^(q) using the response function 

Ga,{z - z') ^ G.,(q) = (3.11) 

z' m 
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of the free manifold. Thus the displacement correlations are given in linear 
response by 



(uMM-^)) = Ga-yicL) Fr(q)Fr(-q) ^^^(-q) = ^S^f,, (3.12) 

where we introduced the variance of the pinning force A^^^ via 

A(2)5„^ = -a«a^A(x)|,^o- (3-13) 

Its value is related to the variance and correlation length of the pinning 
potential approximately through 

A(2) = ^ . (3.14) 

(27r)^/2^2+iV W--^'*^ 

Here we have assumed a Gaussian form for A(x). 

The perturbative correlation function (3.12) is characterized by a rough- 
ness exponent 

Crf = (3.15) 

which we refer to as the 'random force' value. This exponent characterizes the 
actual correlation function only on sufficiently small scales |z| < below the 
Larkin length (Larkin 1970), since the perturbative treatment is justified 
only as long as W{z) < The Larkin length can be estimated by equating 
the elastic energy E^i ~ e^^L^^^^"^ with the pinning energy in the random- 
force approximation Epi^ ~ (L^^^A*^^))^/^ as (Larkin 1970, Bruinsma and 
Aepph 1984) 

L, « [eie/A^'r-, (3.16) 

where € = A — D (to be distinguished from sy). Prom this rough analysis the 
manifold is found to be flat in D > 4 (since Crf < 0), logarithmically rough in 
D = A (with (j-f = 0), and rough with an exponent Crf > in < 4. Although 
the perturbative analysis can provide a good approximation only on small 
length scales |z| < (however, in certain cases intermittency can be relevant 
on small scales (Bouchaud et al. 1995)), the roughness of the manifold in 
D < 4 persists in more sophisticated approaches (such as a self-consistent or 
renormalization-group analysis), which go beyond a perturbative approach. 

Thus, a single vortex line {D = 1) in a bulk superconductor {d = 3), 
which is roughened by pure thermal fluctuations (cf. equation (3.7)), is also 
roughened by disorder at T = 0. In contrast, a VLL in a bulk superconductor. 
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which can be considered as a manifold with D — d, is not roughened by pure 
thermal fluctuations but by pinning. 

3.1.1c Flory analysis 

As already mentioned, the above perturbative analysis breaks down on 
length scales |z| ^ because perturbation theory does not adequately treat 
a system with many minima in the potential energy (for a pedagogical exam- 
ple see Villain and Semeria (1983)). Although perturbation theory could be 
continued to higher orders, it will never be able to describe the displacements 
on largest scales |z| — > oo, where the multi-stability (existence of many local 
minima) of the potential energy landscape is crucial. 

For the further analysis it is convenient to start from the replica Hamil- 
tonian (Edwards and Anderson 1975) 

'Hn — ^el,n ~l~ ^pin,n 

= E/^""^ ||r''(V.u'^)2-^A(u»(z)-u^(z))|, (3.17) 

a, 6=1 ^ 

which is obtained from replicating the original system n times and performing 
a disorder average. Note that Greek lower indices denote transverse spatial 
components, whereas Roman upper indices denote replicas. The analysis 
of this Hamiltonian is highly non-trivial, since the displacement enters the 
argument of the disorder correlator A. A dimensional analysis shows that A 
has to be retained in its full functional form and may not be represented by 
a truncated Taylor expansion in D < 4 (Balents and Fisher 1993). 

Physically, the structure on large scales is determined by a competition 
between elastic and pinning energies. The Flory argument (Imry and Ma 
1975, Kardar 1987, Nattermann 1987) allows one to obtain an improved 
value for the roughness exponent ( by requiring both energy contributions 
to scale with the same exponent on large scales |z|. To be more precise, we 
rescale z = Lz' and u(z) = L''u'(z'), where L denotes the (variable) length 
scale on which we consider the system. Then Ti.ei,n ~ with an energy 
scaling exponent 

e = 2C + D-2. (3.18) 

The fluctuations of u' arc determined by an effective temperature T', which 
has to be rescaled like the Hamiltonian, 

T = T'L^ (3.19) 

in order to keep the Boltzmann factor e~^/^ scale invariant. 
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In pure systems it is possible to achieve scale invariance not only of the 
Boltzmann factor e~^«i/^ but of both cHei and T by the choice C = Cth = 
and ^th = 0. 

At zero temperature, weak disorder is a relevant perturbation in D < 4. 
Since a short-ranged disorder correlator should essentially scale like an N- 
dimensional 5-function, the replicated pinning energy scales as Hpm/T ~ 
j^D-NC-2e ^ rpj_^g exponent of the Boltzmann factor includes now two 
terms, which after rescaling behave as Tid/T ~ l-L'^/T' = 0(1) and 7ipin/T = 
0(L^~^''~^^). In the limit of vanishing temperatures, a finite width of the 
distribution of {u'} is only possible if ij^-^^C-^^ = 0(1), i.e., if 

which imphes 9 = 2( + D - 2 = ^{D - D^), where = Ahhough 
(p is an improvement over this result is (in general) not exact, since the 
scaling behaviour of Hpin was over-simplified. 

On the other hand, at finite temperatures e~"^p'"'"/^ is a relevant pertur- 
bation to the Boltzmann factor e~^«i/^ if D — NQ^ — 2^th > 0, i.e., for 

2N 

which is equivalent to Cf > Cth- Note that < 2 ior N < oo and D^o — 2. 
Prom this observation we conclude that for D < Djv weak disorder is irrel- 
evant. Since on the other hand disorder will certainly become relevant if it 
is sufficiently strong, A > Ac(T), one expects in this transition from 

an unpinned phase for weak disorder to a pinned phase at strong disorder 
(which is equivalent to a thermal depinning transition for increasing temper- 
atures). Contrary to the result for (f, which is an approximate expression 
for the true roughness exponent the result for is exact. 

3.1. Id Renormalization group analysis 

The algebraic roughness (3.6) of the manifold means that it is scale invari- 
ant on large length scales. Hence it does not have a finite correlation length 
and the manifold can be considered as being in a critical state with ( corre- 
sponding to a critical exponent. In analogy to ordinary critical phenomena, 
a renormalization group (RG) analysis is suitable to describe the large-scale 
features of the system going beyond perturbation theory and self-consistency 
arguments. 

In Z) < 4 it is not possible to describe pinning by a finite set of parameters 
since C > and a Taylor expansion of A(u) would yield terms the relevance 
of which increases with increasing order of the expansion. Therefore one has 
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to use a functional renormalization group analysis for the present problem, 
which was established by Fisher (1985a, 1986a) and Balents and Fisher (1993) 
to first order in e = 4 — D. On increasing length scales L oc e' the system can 
be described by a renormalized temperature and disorder correlator, which 
flow according to (Fisher 1986a, Balents and Fisher 1993): 

diT = -er, (3.22a) 
diA{u) = (e-4C)A(u) + Ciia5„A(u) 

+^dadpA{u)dc.dpAiu)-dadpA{u)dad^A{0). (3.22b) 

Here we introduced a rescaled correlator A/ca ^ A with some non- universal 
constant ca that depends on the short-scale cutoff and the stiffness constant 

We will see that 9 — D—2+2( > such that equation (3.22a) imphes that 
the effective temperature vanishes on large scales. The system is therefore 
described asymptotically by a 'zero temperature' fixed point. The fixed-point 
correlator and exponent have to be determined numerically from equation 
(3.22b). 

The actual value of the roughness exponent ('random manifold' value) can 
be represented introducing a correction factor i^(i?, N) in the Flory expression 
(3.20), ^ ^ 

^'"^ ^ A + u{D,N)N- ^^-^^^ 

In all known cases this correction factor lies in the range < v{D,N) < 1. 
The findings of the functional RG analysis to order e — A — D (Fisher 1986a, 
Balents and Fisher 1993) are equivalent to i/(4, 1) 0.800(3). It is interesting 
to note that the roughness can be determined exactly for special dimensions: 
z/(D, oo) = 1, i.e., the Flory exponent becomes exact for an infinite number 
of displacement components (Mezard and Parisi 1990, Mezard and Parisi 
1991), and i^(l, 1) = 0.5 corresponding to C = f (Huse et al. 1985). From 
a more involved field-theoretical self-consistency argument Lassig (1998b) 
obtains z^(l,2) = | corresponding to C = | and i^(l,3) = ^ corresponding 
to C = ^. 

Remarkably, there is no renormalization of the elastic constant e\\ in the 
Hamiltonian. The renormalized elastic constant of the manifold can be deter- 
mined from its response to tilt field /i*^, which is coupled to the displacements 
through 

n^^- j d^'z iJ,ld,u^iz). (3.24) 
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In the absence of a pinning potential a constant tilt field leads to a manifold 
displacement ii*(z) = iJi,*^Zi/e\\ and to a change of the free energy by an 
amount 

^ -m--Jd''z (3.25) 

Even in the presence of the pinning potential with a stochastic translation 
symmetry, this identity holds exactly. This is due to a 'statistical symmetry' 
of the pinning energy under a transformation u(z) u'(z) = u(z) + Su{z) 
for an arbitrary function 5u{z) (Goldschmidt and Schaub 1985, Schulz et al. 
1988, Balents and Fisher 1993), provided disorder has a vanishing correlation 
length in z directions, as assumed in equation (3.3). In particular, the choice 
Sua{z) = /j*^Zi/e\\ transforms the pinning potential V^(X + u(z), z) — >■ V(K + 
u(z)+(5u(z), z) = y (X+u(z), z), which has the same statistical properties as 
the original potential. This stochastic symmetry is reflected most obviously 
by 'Hpin,n in equation (3.17), which is invariant under the transformation 
u"(z) — > u"(z) + 6u{z) that is identical in all replicas. The non-renormali- 
zation of e\\ follows directly from (3.25), since renormalized elastic constants 
are defined by 

the change of the free energy and the dependence of ^ on is independent 
of disorder. Here — J dPzl is the system size. 

If one starts from a disorder with a finite correlation length also in z 
direction, there will be a finite renormalization of the elastic constants due 
to fluctuations on these small length scales. Beyond this scale, the manifold 
behaves as for vanishing correlation length. Therefore the asymptotic be- 
haviour on large scales (and in particular the roughness exponent) will not 
be affected by such correlations. 

The roughness exponent is, however, sensitive to the behaviour of A(x) 
on large scales, if this correlator is long ranged. Although this scenario is not 
of immediate interest for the present consideration of pinning by point like 
defects, it is realized, for example, for pinning of vortex lines by columnar 
defects or for magnetic domain walls with random fields. In this case (where 
N — 1), A(x) — A(0) oc \x\ and the roughness exponent has the 'random-field' 
value (Villain 1984, Grinstein and Fernandez 1984, Bruinsma and Aeppli 
1984) 

Crfi = ■ (3.27) 
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Surprisingly, this exponent characterizes the manifold in the presence of a 
driving force at the depinning threshold, as we will discuss below. 

S.l.le Pinning vs. thermal fluctuations 

In section 3.1. Id, we have described the structure in the disorder-domi- 
nated regime, for the 'zero-temperature' fixed point. Now we reconsider the 
effect of thermal fluctuations in the presence of disorder. For this purpose 
Fisher and Huse (1991) separated the displacement correlation function (3.6) 
into two contributions, the first one due to disorder 

W^pin(z - z') = (u(z) - u(z'))2, (3.28) 

and a second one 

Wth(z-z') = 



describing thermal fluctuations. In the absence of disorder, l^pin(z — z') = 0. 

In particular Wth is found to be independent of disorder due to the sta- 
tistical symmetry (Schulz et al. 1988, Fisher and Huse 1991) already men- 
tioned above. Thus VFth(z) ~ \z\'^'^*^ with the exponent (3.7). Consequently, 
Wpin(z) ~ |zp^"" in the low-tcmperaturc phase, where (^rm > Cth- 

It was argued (Fisher and Huse 1991, Hwa and Fisher 1994a, Kinzel- 
bach and Lassig 1995) that the manifold would have (with probability 1) a 
unique ground state. Then VFpin essentially characterizes the ground state 
of the manifold (Nattermann 1985a, Huse and Henley 1985) and Crm would 
be the roughness exponent thereof. However, a small fraction of samples 
(or, small areas within a sample) will have nearly degenerate excited states 
with large displacements relative to the ground state. Although such exci- 
tations are rare, due to the large displacements they can dominate several 
disorder-averaged thermodynamic quantities (Nattermann 1988, Hwa and 
Fisher 1994a). In particular, these rare fluctuations are responsible for the 
growth of the thermal width Wth in D < 2. Although Wth shows exactly the 
same behaviour as in the absence of disorder, it is important to emphasize 
that the distribution of thermal fluctuations is highly non-Gaussian (Fisher 
and Huse 1991, Hwa and Fisher 1994a). 

For D > 2 thermal fluctuations have only a finite width and the equi- 
librium state of the manifold globally reflects the ground state. The broken 
translation of the Hamiltonian are reflected by the equilibrium state of the 
manifold. 



Wpin{z - Z 



([u(z) - u(z') - (u(z) - n{z'm 

T . 



— z 



J\2-D 



(3.29) 
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In summary, the structure of elastic manifolds in weak disorder strongly 
depends on its dimensionalities and shows the following behaviour (see also 
figure 3): 

• For D > A the manifold is flat (C < 0) at all temperatures provided 
disorder is weak. Sufficiently strong disorder will always induce rough- 
ness. 

• For 2 < D < A disorder roughens the manifold {( — Crm) at all tem- 
peratures. The manifold stays close to its ground state; the second 
moment of thermal displacements is finite. 

• For2>i:)>i:)jv = 2N/(2 + N) disorder still roughens the manifold 
(C = Crm) at all temperatures. Now also the second moment of thermal 
displacements is infinite. 

• For D < Djv and sufficiently high temperatures, the manifold is en- 

tropically driven out of the ground state and shows a structure similar 
to that in the absence of disorder {( = (ti^). However, there is also 
a low-temperature phase where disorder is relevant (Cook and Der- 
rida 1989). In other terms, the manifold shows a temperature-driven 
depinning transition at a finite temperature. 

For more details, the interested reader is referred to Halpin-Healy and 
Zhang (1995) and Lassig (1998a). 

3.1.2 Thermodynamics 

So far we have described the relevance of disorder for the structure of the 
manifold. In particular, the roughness is qualitatively increased in the dis- 
order dominated phases. Could such an increased roughness be taken as an 
unambiguous sign for glassiness of the state? To understand this problem 
better, we study the possibiUty of anomalous scaling properties of thermo- 
dynamic quantities such as energy, entropy, and free energy. 

Since the elastic energy and the pinning energy are local quantities, en- 
ergy, entropy, and free energy scale extensively with the system size L: 

Er^L^, 5~L^, Tr^L^. (3.30) 

Due to the randomness of pinning there are important sample-to-sample 
fluctuations with a scaling (for T > 0) 

(A^~L^, (A5)2~L^, {E^^L'^ (3.31) 
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one flat phase: C=C,r{ 



one rough phase: ^=^, 




two rough phases: 
high T : disorder irrelevant: ^=^th 



N 



Figure 3: Illustration of the relevance of disorder for the roughness of 
the manifold in the {D,N)-plane (as described in the text). 



which has been confirmed numerically (Fisher and Huse 1991). Here 6 is 
again the energy scaling exponent (3.18). The fluctuations of AS = S — S 
and AS = S — S are normal and dominated by fluctuations on small scales: 
each sub-volume of the size contributes independently. However, the 
fluctuations of AT = T — T are anomalous and governed by large-scale 
contributions. Clearly, for T ^ there exists a diverging length scale below 
which the fluctuations of AE scale as those of AT . 

In thermodynamic equilibrium the manifold minimizes its free energy. 
Therefore energy and entropy are not minimized independently, their leading 
order cancels (at T > 0), and the free energy fluctuations can be smaller than 
those of the energy. One expects 6* < y in all dimensions, which corresponds 
to C ^ ^i-^'^) > in equation (3.23). Indeed, the exponent ^im = 

(4 — (i)/4 follows from an Imry-Ma type argument (Nattcrmann 1985a) and 
is considered to be an upper bound for ^rm in random-bond systems (Fisher 
1993). At strictly zero temperature E — T and both quantities must scale in 
the same way, AE = AT ~ . Therefore temperature plays a crucial role, 
although it seems to be irrelevant in the flow equation (3.22a). This is because 
temperature is a dangerously irrelevant variable (Fisher 1986b). In very 
rare samples there are excited states which are nearly degenerate with the 
ground state but which deviate from it by large displacements. Such rare but 
large fluctuations still give dominant contributions to many thermodynamic 
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disorder-averaged thermodynamic properties, for a more detailed discussion 
see Fisher and Huse (1991) and Hwa and Fisher (1994a). 



3.1.3 Susceptibility 

In order to find an unambiguous signature of glassy phases, Hwa and Fisher 
(1994b) proposed to examine the susceptibility of the system with respect to 
a tilt field. Although in this section wc are dealing with directed manifolds, 
we estabhsh the analogous susceptibility for these simpler systems and show 
how it is related to the displacement correlation functions discussed above. 

For simplicity we restrict the discussion of susceptibilities to = 1; the 
general case > 1 follows from a straightforward generalization. Wc couple 
the manifold to an inhomogeneous tilt field /x(z) via an additional energy 
contribution as in (3.24), 

= - / d^'z [/u* + /u(z)] • V\iu(z). (3.32) 

As in section 3.1. Id, n* is a constant field applied to measure the tilt response 
of the system. 
The system 

= |d^^||(vp-l[/x* + /x]y-^[M* + /i]^|, (3.33) 

without random potential V but with a random tilt /x(z), has a ground state 
u*{z) which couples only to the longitudinal part Hi of the tilt field and 
which is determined by fi* + fj,{z) = V||U*(z). For a random tilt with zero 
average and long-ranged correlations 



lii{z)iij{z') ~ |z - z I"", a<D, (3.34) 

the manifold has in its ground state an exponent ( = |(2 — «) and will 
thus be rough if a < 2. Nevertheless, since this toy model is harmonic and 
trivial (it is bilinear in u, and H' is even translation invariant), it cannot be 
considered as glass, which motivated Hwa and Fisher (1994b) to search for 
unambiguous signatures of glassiness. 

They proposed to examine the sample-to-sample fluctuations of the linear 
response susceptibility Xap of the system for spatially constant tilt /Lt*, which 
can be obtained from the free energy of a particular sample by 

1 d^J^lfi* 



X 



(3.35) 
lx*=o 
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Then it is interesting to consider its sample-to-sample fluctuations 

(A^ = P-X'. (3.36) 

The toy model Ti' = Tici + that has only random tilt fi but no random 
potential, has % = l/ey for all realizations of the random tilt because of the 
statistical symmetry discussed after equation (3.25). Hence the vanishing 
of (Ax) 2 = reflects its trivial nature, i.e., that random tilt fields only 
deform the ground state, and thermal fluctuations around the ground state 
are similar to those in the absence of disorder. 

Now we come back to the manifold model in a random potential (we drop 
the random tilt but we still assume N — 1 ior simplicity) to demonstrate 
that (Ax)2 > 0, i.e., that susceptibility fluctuations are a less ambiguous 
signature of glassiness than mere roughness. Due to the statistical symmetry, 
— JF[0] = — L-^27^(/x*)^ and x = l/£\\ exactly. We can calculate its 
susceptibihty fluctuations perturbatively from the free energy fluctuations at 
T = 0, 

A:F[/li*]A:F[/li*'] = J d^'zA (^^(m* - A**') • , (3.37) 

involving the potential correlator A(x) evaluated for the tilted manifold in 
the absence of pinning. This results in 

-dA(^)(0) 



(Ax)2 oc L^-^^^, (3.38) 

where we have dropped numerical factors and A*^^)(0) = d^A{u)\^^Q. The 
dependence of this result on the system size L shows again that (A^)^ ^ 
in Z) < 4. Certainly, we cannot expect this pcrturbativc result to be 
quantitatively correct for large L, but qualitatively it reflects the relevance 
of disorder and the glassiness of the manifold. In D < A one expects (Ax)^ 
to be finite. Its correct value has to be calculated within the RG scheme (see 
Hwa and Fisher (1994b)). 

We will not further pursue an accurate calculation of the susceptibility 
fluctuations here. Instead we establish its connection to the displacement 
fluctuations. To this end we introduce the response of the local tilt dau{z) 
to the held Hpiz!) (which is now considered as an external control parameter 
rather than a form of disorder) 

, _ 5(9„ii(z)) 52 

= ^[{dauiz)dpuiz'))-{do.uiz)){dpuiz'))], (3.39) 
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which is now exphcitly expressed as a displacement correlation function and 
can be related to correlations of the order parameter ip exploiting 

(3.40) 

k=0 

The total susceptibility (3.35) discussed above can be obtained as integral of 
the local susceptibility (3.39), 

X-^j d'^zd^z' Xaa(z,z'). (3.41) 

(Note that here is no factor 1/D in contrast to (3.35) since the displacement 
responds not only to the longitudinal component of /^(z) but to the entire 
/i*.) For a particular sample the local susceptibility explicitly depends on 
two space coordinates z and z'. We introduce a susceptibility Xapi'^ — z') = 
J d^Zo Xapi'^ + Zo,z' + Zo) averaged over the volume Zq G L^, which 
now depends only on the coordinate difference z — z'. Then x is conveniently 
related to the Fourier transform of x(z) through 

X= Xaa(k)|k=o- (3-42) 

(For example, the toy model has Xa/diX) = (V^||)-Pi/3(k) with the longitudinal 
projector P^^(k) = kak^/k"^ and hence x = V^llO 

The fact t hat x has sample-to-sample fluctuations in a glassy phase, 
lim^^^oo (^x)^/x^ 7^ 0, means that the susceptibility is not self-averaging 
(Aharony and Harris 1996). Since x is obtained from a volume average 
sample-to-sample fluctuations also mean that the translational average is 
not equivalent to the average over many samples. 

From the expression for the susceptibility in terms of the displacement 
fleld, equation (3.39), one recognizes that the susceptibility essentially mea- 
sures fluctuations around the ground state of the system. More precisely, one 
can relate the susceptibility to the thermal displacement correlation (3.29) 

by 

d^d'f.Wti.iz - z') = -2T Xap{z-z'). (3.43) 

Thus the disorder-averaged susceptibility is related to Wth- The fact that 
X is independent of disorder is related to the disorder independence of VFth 
discussed in section 3.1.1e. A signature of glassiness can therefore appear 
only in the sample-to-sample fluctuations of x- 

To conclude the discussion of susceptibility fluctuations as characteristic 
for glassiness, we wish to point out that they cannot be taken as sufficient 
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criterion for glassiness. This can be demonstrated by a counterexample with 
a 'pinning' energy 

7ipin = j d^z^[u-Vi\u{z)]\ (3.44) 

where u is a. vector of fixed length \ u\ but with random orientation. This type 
of disorder actually represents randomness in the elastic constants and the 
model is not a glass since it is Gaussian in the displacement. Nevertheless, 
it has sample-to-sample fluctuations 

(A^=[(l + i.Vqi)^/2-l]x' (3.45) 
that are finite for i/^ > 0. 



3.1.4 Barriers 

The glassy nature of a system is generally related to an extremely slow dy- 
namics that is dominated by thermally activated processes in a complex en- 
ergy landscape with many meta-stable states. Glassy systems have not only 
a disorder-dominated ground state, but also a huge number of meta-stable 
states. This section is devoted to the description of the energy landscape. 

The dynamical behaviour is determined by 'neighbouring' meta-stable 
states, which arc related to each other by the slip of a restricted part of the 
manifold over a certain barrier. The part of the manifold is assumed to have 
a certain length in z-directions, the displacement in this region is of a 
magnitude L-^ and the barrier height is denoted by U, see figure 4. 

A rigorous characterization of the statistics of barriers is very intricate. 
Therefore, a scaling picture (Villain 1984, Nattermann 1985b, Huse and Hen- 
ley 1985, loffe and Vinokur 1987) has been put forward, where it is assumed 
that the statistics of barriers is essentially identical to the statistics of the 
free energy fiuctuations. Therefore the barrier height should scale with the 
barrier length like 

U{L,)^U^(^^^ , (3.46) 

where = U{L^) ^ is the typical height of the smallest possible 

barrier of the size of the Larkin length. The scaling exponent for the barrier 
height %Ij is assumed to coincide with the exponent from equation (3.31) 
that describes the free energy fiuctuations: 

ijj^e = 2C + D-2. (3.47) 
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Figure 4: Illustration of two meta-stable configurations of the manifold 
(solid lines) that are separated by a barrier (shaded area). The charac- 
teristic length scales of the barrier are in longitudinal direction and 
L-r in transverse direction. The height of the barrier is U . 

This exponent is related through equation (3.18) to the roughness exponent 
C of the manifold, which describes the scaling relation 



between the lateral and transverse sizes of the barrier region. 

The basic assumption ip = 9 has been confirmed by analytic arguments 
combined with numerical simulations (Mikheev et al. 1995, Drossel and Kar- 
dar 1995, Drossel 1996). However, the statistics of barriers and free energy 
fluctuations have turned out to possibly be not strictly identical: U{Lz) and 
{[AJ^{Lz)]'^Y^'^ can differ by factors that are powers of InL^ which do not 
modify the exponent relation (3.47). 

It is clear that not all barriers of a given size have strictly the same 
height but that they are distributed according to a certain distribution 
Vl^{U). The above scaling relations are valid for the 'optimal' (dynamically 
relevant) barrier height, which is the average barrier obtained from VlX^) 
with an additional weighting factor for the density of barriers of the size con- 
sidered (Vinokur et al. 1996). In particular in low dimensions the dynamics 
of the manifold is not necessarily dominated by the average barrier height, 
but by the largest barriers. Therefore it is also of fundamental interest to 
characterize the distribution Vl^{U) for large U. 




(3.48) 
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For a string, i.e. {D, N) = (1, 1), Vinokur et al. (1996) found 



VlAU) ^ -^^e-^/^« exp[-(L,/L^)e-^/^«] (3.49) 

from a combination of extreme value statistics (Gumbel 1958, Galambos 
1978) and a coarse graining approacli. Large barriers are exponentially rare 
since Vl^U) ~ c^i^e-^/^f for U ^ oo. Gorokhov and Blatter (1999) found 
that the free-energy distribution obeys a similar exponential decay at large 
negative values of the free energies J-'. The free-energy distribution decays 
much faster at large positive JF than at large negative JF, (Gorokhov and 
Blatter 1999), which one can expect because the manifold minimizes the free 
energy in equilibrium. 



3.2 Transport properties 

The glassy nature of systems does not manifest itself only in equilibrium 
properties but, maybe even more unambiguously, in its dynamic behaviour. 
The existence of many meta-stable states and the thermally activated nature 
of transitions between these states makes the dynamics of the glass extremely 
slow. We describe now the transport properties of manifolds in a stationary 
driven state, leaving aside interesting topics such as relaxational dynamics. 

Vortices move dissipatively in the superconducting condensate. For small 
velocities the friction force is proportional to the velocity with a constant 
viscous drag coefficient tiq. For a single vortex line in a bulk superconductor 
(Bardeen and Stephen 1965) 

rjo ^ —, 3.50 

where pn is the normal state resistivity and c the velocity of light. On large 
time scales inertial forces can be neglected in comparison to the friction force. 

The equation of motion is then obtained from balancing the friction force 
(more precisely: the force density in a I?-dimensional space; for brevity forces 
and force densities are not strictly distinguished) with the forces arising from 
the interaction with other vortices, the pinning force, the driving force F and 
a thermal noise 

^x=-^ + F + C (3.51) 
The noise is taken as Gaussian distributed with zero average and variance 

iUz, t)C)3{z', f)) = 2r]oTS^fsS{z - z')6{t - t') (3.52) 
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such that the equation of motion properly describes thermodynamic equilib- 
rium according to the fluctuation-dissipation theorem in the absence of the 
driving force. 

At finite temperatures the manifold will respond to the driving force with 
an average velocity v = (x). For calculational simplicity it is advantageous to 
consider v as prescribed and to calculate the driving force F = F(v) required 
to maintain this velocity. In the driven case we define the displacement 
u(z) = x(z) — X — such that (u) = 0. This displacement follows the 
equation of motion 

r/ou = — i + F-?7ov + C 

ou 

= £|,V^u- Vx^(X + vt + u,z) + F-r;oV + C, (3.53) 



which serves as basis for the following analysis [recall that Vy 
■ ■ ■ ' sfe) ^-L = (afl' ■ ■ ■ ' d^)^- 



3.2.1 Friction 

In the absence of pinning, V = 0, the manifold moves with constant velocity 
V = F/rjo. Pinning tends to slow down the motion of the system, because 
the manifold has to overcome barriers and it loses more time by sliding uphill 
than it wins by sliding downhill. Thus, to establish a certain velocity v in 
the presence of pinning, one has to apply a force F(v) which is larger than 
in the absence of pinning. In general, the velocity-force characteristic will be 
non-linear but one can define an effective zero-velocity friction coefficient by 



dv 



(3.54) 



v=0 



We now show that this 'renormalized' friction coefficient can diverge due to 
the presence of disorder, which is another signature of glassiness. 

The divergence of the effective friction coefficient is examined by treating 
disorder perturbatively. We make use of the linear response function of the 
disorder-free system 

G^^{z,t) - SFp{0,0) - noWe^^t) " II U(t)()„^, (3.55) 

which describes the reaction of the displacement to a locally perturbating 
force. 
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By iterating the equation of motion (3.53) up to 0{V^), the condition u = 
yields the velocity-force characteristic and the effective friction coefficient 

Vap = r]oSai3 + J dtt G^s{'z^O,t)dadpd^dsA{0). (3.56) 

Inspecting the large-time behaviour of the response function (3.55) for z = 0, 
this coefficient diverges in D < 4 as long as dadpd^ds A{0) ^ for certain 
indices. The divergence of Ty^/j resembles that of the susceptibility fluctuations 
(3.38): both couple to the foTirth derivative of A(u). This divergence actually 
does not mean that a stationary state with non-vanishing velocity could be 
established only by an infinite driving force; it means that at small velocities 
v(F) is no longer linear (Ohmic) but sub-linear. 

The quahtative change of the transport characteristic in < 4 is consis- 
tent with roughening as a qualitative structural change. This consistency is 
not just a fortunate coincidence. One can actually show the relation between 
dynamic and static features using a fiuctuation-dissipation relation between 
the correlation and response function of the free system (Scheidl and Vinokur 
1998b). Nevertheless, the toy model studied in Sec 3.1.3 has a strictly hn- 
ear transport characteristic (the same as in the absence of disorder) despite 
its roughness. The transport characteristic is not at all affected by random 
bonds coupling to the displacement as in equation (3.32) because the trans- 
lation symmetry u — u -|- xq is not broken. Thus the qualitative sensitivity 
of the transport characteristic to the nature of disorder suggests its use as a 
supplementary indicator for glassiness. 

If the dynamics of the system is glassy, the transport characteristic shows 
several characteristic features, cf. figure 5. For small driving forces it is 
important to distinguish between the cases of zero and finite temperature. 
At T = the driving force has to exceed a threshold value Fc to set the 
manifold in motion with a finite average velocity, since the driving force has 
to overcome the average pinning force. (Strictly speaking, this is true only 
in the absence of quantum fiuctuations, which can lead to quantum creep 
at T = and F < Fc, see Blatter et al. (1994) and references therein.) 
This threshold phenomenon is called depinning. For T > a finite velocity 
is found even for |F| < Fc due to thermal activation, and the transport 
characteristic will be sub- linear {creep regime). At large driving forces |F| ^ 
Fc the effect of pinning diminishes and the free differential mobility is reached 
asymptotically, dFa(y)/dvij — > SapTjQ. This latter regime is called the flow 
regime. 

The divergence of rjap implies the break-down of perturbation theory for a 
calculation of the transport characteristic. Although it is difficult to calculate 
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Figure 5: Illustration of the dynamic regimes (shaded areas) in the 
transport characteristic v{F) of a manifold driven in disorder. In the 
absence of disorder the characteristic is linear (dashed line). In the 
presence of disorder it is non-lineax (solid lines). At T = the manifold 
depins as a critical force Fc. At T > thermally activated creep occurs 
already for F Fc. For F S> Fc the manifold flows essentially as in the 
absence of disorder. 

the global form of the characteristic, it is possible to achieve an analytic 
description for small and large velocity. 

3.2.2 Depinning 

The value of the critical force F^. where motion sets in at T = can be esti- 
mated as follows: On length scales below the Larkin length the manifold 
performs only displacements smaller than ^, such that the dependence of the 
pinning force on the displacement can be neglected. In such a volume the 
local pinning forces add up to = I\z\<l^^^'^ [— V_Ly(X, z)]. The forces 
from different volumes are essentially uncorrelated. Thus depinning, where 
the driving force is balanced by pinning force, is found at (Feigelman 1983, 
Bruinsma and Aeppli 1984) 

Fc ^ L-^ {fif' (A(^))V2L-^/^ ^ silCLf ~ (A(^))^A, (3.57) 

where again e = 4 — D and the pinning force variance was given m 
equation (3.14). For fixed but weak strength of A^^^ the critical force Fc 
vanishes as D y 4 and it is zero for D > A. 

Thus for weak disorder (de)pinning occurs in D < 4 dimensions only, to 
which we restrict our analysis. The velocity is expected to increase continu- 
ously when the driving force exceeds the threshold value. This observation 
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was made originally by Middleton (1992a) in the context of charge-density 
waves. Thus it should be possible to consider depinning as a continuous 
transition, at which the manifold exhibits critical fluctuations (Nattermann 
et al. 1992). The depinning transition was studied theoretically for directed 
manifolds (Nattermann et al. 1992, Narayan and Fisher 1993, Erta§ and Kar- 
dar 1996, Leschhorn et al. 1997) and for the closely related charge-density 
waves (Fisher 1985b, Narayan and Fisher 1992a) within a renormalization 
group approach. Kardar (1998) and Fisher (1998) review depinning and its 
relevance in a much wider variety of phenomena. 

The phenomenology of the transition is described qualitatively by the 
following scenario: as long as the manifold is pinned, it will be in a rough state 
because of competition between elastic forces, pinning forces, and the driving 
force. In the absence of the driving force, the roughness is described by the 
exponent d-m = Cmi(O) of manifolds in a random potential at equilibrium, 
which was described in equation (3.23). A driving force polarizes the state 
of the manifold and turns the situation into a non-equilibrium one, where 
the exponent is modified to Crm(-^)- Preliminary numerical simulations (Leaf 
et al. 1999) indicate that for forces |F| < Fc below the depinning transition, 
Cmi(-^) continuously increases with the driving force, until it reaches the value 
Crm(-fc) = Crfi of maulfolds iu a random force field at equilibrium, which was 
described in equation (3.27). This means that the effective force acting on 
the manifold has qualitatively changed its character. Above depinning, where 
the manifold moves with a finite velocity, the force acts on large scales like an 
effective thermal noise (although we consider the situation at T = 0, without 
thermal noise C, in the equation of motion). Accordingly, the manifold will 
be flat on large scales, as long os D > 2. 

For small velocities, the saturation of the displacement correlation func- 
tion occurs on a large length scale £z, which acts as a dynamical correlation 
length that diverges at the transition as (Middleton 1992a, Nattermann et al. 
1992, Narayan and Fisher 1993) 

4~(F-Fe)-^ (3.58) 

with a certain exponent u. For this zero-temperature non-equilibrium transi- 
tion, the force acts as a control parameter that plays the role of temperature 
in a usual finite-temperature equilibrium transition. Associated with the 
longitudinal scale 1^ there is a transverse scale 

4 ~ (3.59) 

with being the roughness exponent at depinning. As the transition is 
approached from above, F \ F^., there is also a diverging time scale t^, the 
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correlation time of the pinning force acting on the manifold, 

ty ^ —. (3.60) 

V 

The dynamical exponent z then can be defined from the scaling relation 
between and ty, 

ty ~ il. (3.61) 

A further scaling relation describes the continuous onset of motion, 

vr^{F- F^f. (3.62) 

Prom a balance of the elastic and driving force in the equation of motion 
(making use of the non-renormalization of the elastic constant) one derives 
(Nattermann et al. 1992, Narayan and Fisher 1993) 

ly = (3.63) 

The consistency of equations (3.60), (3.61) and (3.62) requires the scaling 
relation 

P^{z- Cy (3.64) 

In the general case with more than one displacement component (for 
N > 1), the velocity selects one particular direction out of the A^- dimensional 
space and we may choose the coordinates such that v = v6a^i. Due to 
this selection of a direction by velocity one has to expect the displace- 
ment correlations to be uniaxially anisotropic. We subsequently distinguish 
Wa{z,t) — {[ua{z,t) — iiQ,(0, 0)]^), where a = 1 for the parallel component 
and a = 2, . . . , A^ for the perpendicular component. The anisotropy is not 
a minor quantitative effect, it concerns the scaling exponents as shown by 
Ertas and Kardar (1996). The anisotropic scaling laws for the displacement 
components are 

W^{z,t) = |z|2C«^„(V|zr«), (3.65) 

with g'a(O) finite and gaiv) ~ for large arguments y. The above scaling 

relations (3.63) and (3.64) then hold with the identification 

C = Ci, z = zi. (3.66) 
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The values of the exponents which can be calculated analytically within a 
renormalization group approach, using an expansion in small e — 4 — D 
(Nattermann et al. 1992, Erta§ and Kardar 1996) are: 




(3.67b) 



(3.67d) 



(3.67a) 



(3.67c) 



and all directions a >2 scale identically. Equation (3.67a) means that the 
manifold in a random potential at depinning has the roughness of a manifold 
in random force field in equilibrium, cf. equation (3.27). Narayan and Fisher 
(1993) proved that equation (3.67a) is correct to all orders in e. 

For the depinning of periodic media such as charge-density waves the 
same scenario holds as for directed manifolds. However, they belong to a 
different universality class with different values of the exponents (Middleton 
1992a, Narayan and Fisher 1992b, Narayan and Fisher 1992a). 

So far, we have discussed depinning as a phenomenon at T = 0. The 
effect of finite temperature is that potential barriers can be overcome by 
thermal activation. Therefore a finite velocity is expected for every finite 
driving force. Fixing the force at the T = threshold value F = Fc, velocity 
increases with temperature according to a power law 



The exponent r depends on the type of the pinning potential in a non- 
universal way (Fisher 1985b, Middleton 1992b). For the 'ratcheted kick 
model' Middleton (1992b) argued for r = 2 in D = 2, 3. The actual value of 
r is still controversial and subject to investigations (Nowak and Usadel 1998, 
Roters et al. 1999). 

3.2.3 Creep 

Below the depinning threshold, for |F| < F^, the dynamics of the manifold is 
in the creep regime, where motion is possible only through thermal activation. 
At |F| = Fc, barriers of all sizes are equally relevant and it is precisely for this 
reason that depinning appears as a critical phenomenon. Below depinning, 
barriers of a typical size are prevalent and dominate the thermally activated 
dynamics. In the presence of a driving force the manifold experiences an 



rp(3/T 



(3.68) 
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effective potential Ves^x., z) = F(x, z) — x • F. For weak forces the manifold 

still has meta-stable states in this effective potential, which are separated by 
barriers that also depend on the driving force. 

Wc now focus on the creep regime of very small forces |F| <^ F^, where 
the barriers are nearly identical to those in equilibrium and where a scaling 
picture for the typical, dynamically relevant barrier has been developed. The 
starting point of this picture is the assumption that barriers of a length 
and free energy fluctuations of a system of size have identical scaling as 
described in equation (3.46). The typical lengths and ~ are related 
by equation (3.48). When the manifold overcomes such a barrier, it gains an 
energy 

Ef ^ FL.Lf ^ F^Lf i^-^J (3.69) 

in the field of the driving force. Therefore, in the effective potential Ves only 
those barriers are still effective for which U < Bp. From U{Lp) = Ep follows 
(loffe and Vinokur 1987, Nattermann 1987, Feigelman and Vinokur 1988, 
Feigel'man et al. 1989, Nattermann 1990) the size of the largest effective 
barrier 

/FX -i/{D+c-e) 

Lp ^ L^{F)^L^i^yJ ^F-VP+<-^) (3.70a) 

U{F) = U{Lp) ^ ~ F-" (3.70b) 



with the creep exponent 



We have used the abbreviation F^ = U^L^^^~^ and the identity (3.18). The 
effectiveness of barriers of a given size Lz is illustrated in figure 6. If one 
denotes by Fb{Lz) = U{Lz) / {L^L^) ~ -^^^^ the threshold force density 
which is required to overcome barriers of size L^, the size of the largest 
effective barrier is determined alternatively by F = FB(LF)■ 
Since the time t to overcome a barrier by thermal activation increases 
exponentially with the height C/, 

t{U) ~ e^/^, (3.72) 

it is the largest effective barrier that dominates the dynamics and determines 
the creep velocity 

v{F) ~ (F/r/o)e-^(^)/^ ~ (F/r?o)e^(^^/^)(^/^^)-^ (3.73) 
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Lf L, 

Figure 6: The height of barriers scales Uke U ^ with the barrier size 
Lz (in general, > at a zero-temperature fixed point). The effective 
pinning force of these barriers scales like Fb ~ Lz (the manifold model 
is valid only as long as (" < 1). In the presence of a driving force F only 
barriers of a size L^ < Lp are still effective. 

The pre-exponential factor has been inserted by hand for dimensional reasons 
and to match the flow regime. However, in the creep regime this factor may 
be of more complicated form (corresponding to a logarithmic contribution to 
U{F)) which is not captured by this scaling argument. 

Prom equation (3.73) we see that the dynamic response will be expo- 
nentially small below depinning, since a disorder-dominated roughness with 
exponent ( > Cth = (2 — -D)/2 implies /i > 0. The characteristic dynamic 
exponent /i depends only on the dimension and the static roughness exponent 
of the system. 

In order to provide a deeper understanding of these scaling relations sev- 
eral authors (Radzihovsky 1998, Balents et al. 1998b, Chauve et al. 1998) 
have recently developed a renormalization group analysis for the whole range 
of velocities. For T = their approach correctly reproduces the properties of 
the depinning transition and for T > they obtain a more precise character- 
istic v{F), which confirms the exponential dependence of the characteristic 
(3.73) to leading order for small forces but contains also subleading orders. 

Experimental observations on magnetic domain walls, charge- density 
waves, and VLL are widely consistent with this scaling picture of barri- 
ers. The most explicit measurement of the roughness and creep exponent 
examining magnetic domain wall creep (Lemerle et al. 1998) was made only 
recently and is in good agreement with the theoretical results. 
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3.2.4 Flow 



The flow regime at large driving force |F| 3> Fc does not usually receive 
much attention since it is considered trivial and is expected to be captured 
within perturbation theory. However, recently the flow regime has attracted 
some attention for driven VLL, where it turned out to be quite non-trivial 
and inaccessible to a perturbative analysis (see section 7) due to an intricate 
interplay of pinning forces and driving force. Therefore we will briefly sketch 
this regime also for the manifolds. 

Let us consider again zero temperature. There the displacements are 
expected to be very small for large driving forces: the manifold can respond 
to the pinning force only with certain (scale-dependent) relaxation times and 
the pinning force changes rapidly in time for large F. Therefore the effect of 
pinning effectively gets averaged out for F — > oo. 

Since the disorder induced displacements are small at large F, one may 
again perform a perturbative analysis as a first step. According to the equa- 
tion of motion (3.53) the effect of disorder is given to lowest order by a 
pinning force 

FP^'^(z, t) = -V±V{X + vt, z), (3.74) 

which shakes the manifold as it moves over the pinning potential. Its average 
is zero and its variance is 



Fr(z,t)F7(0,0) = -9„9^A(vt)<5(z). (3.75) 

For qualitative purposes this shaking force might be compared to an effec- 
tive thermal noise. Because the velocity selects a particular space direction 
this noise in no longer isotropic. Effective shaking temperatures (which were 
introduced by Koshelev and Vinokur (1994) for vortex lattices) for the di- 
rections parallel and perpendicular to v = {v,0) can be identified from 
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(3.76a) 
(3.76b) 
(3.76c) 

These shaking temperatures are to be added to the physical temperature 
to give the total effective temperature, which is finite. Thus, for D > 2 the 
manifold should be flat with exponent ^ = (^th = as in the absence of 
disorder. This flatness is expected to persist at any finite velocity, i.e., for all 
forces above Fc. However, it is interesting to observe that in the flow regime 
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the manifold has a larger width in the direction perpendicular to v since 
> Tf^ (with a > 2). On the contrary, at depinning the manifold has a 
larger width in the direction parallel to v since < Ci, see equation (3.67). 
This means that at a certain velocity the manifold has to reverse the aspect 
ratio of its widths. Apart from this interesting observation, perturbation 
theory gives no hints of qualitative changes of the state with changing velocity 
such as dynamic phase transitions. However, such transitions exist in periodic 
structures such as VLL (section 7). 

3.3 Summary 

In this chapter we have reviewed various properties of manifolds in order as to 
identify a glassy state of the system. Of course, disorder always quantitatively 
modifies the state on small length scales. Some large-scale properties of the 
system can be modified qualitatively by disorder and are hence suitable for 
the characterization of glassiness: 

(i) Disorder reduces the positional order and it increases the roughness. 
Typically, ( > Cth but such an increase is also possible for ( = (^th if W 
is increased only by logarithmic factors. 

(ii) Disorder leads to large sample-to- sample fluctuations, which can be 
described by susceptibilities. 

(iii) In a renormalization-group description of the large-scale physics, 
there is a disorder-dominated fixed point. Typically, this is a zero- 
temperature fixed point. In low dimensions (D < 2), this may also be 
a finite-temperature fixed point. 

(iv) Disorder induces thermally activated sub-linear response to a driving 
force. 

We consider all these features as generic to identify the glassy nature of a 
system. The insufficiency of individual features can be shown by examples: 
An increased roughness alone is not significant as emphasized by Hwa and 
Fisher (1994b) using the bond-disorder model, which has no broken trans- 
lation symmetry. Sub-linear response is a signature of broken translation 
symmetry but can also be achieved by 'ordered' potentials. Similarly, a zero- 
temperature fixed point can, in general, represent a low-temperature ordered 
phase which is not necessarily a glass. The sine-Gordon model provides an 
example for the last two statements. 
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4 Superconducting film in a parallel field 



In treating a many-vortex system we begin with the technically most simple 
problem of a superconducting film in a magnetic field parallel to the film 
plane. The mean-field phase diagram of this problem was considered by 
Abrikosov (1964), who found for Hd < H < Hc2 a solution with an equidis- 
tant vortex lattice of spacing a. Here Hd — In | and Hc2 = where 
s is the thickness of the film and we assume A > s/vr > ,^ (see Sonin (1992) 
for films described by the Lawrence-Doniach model). Such a geometry can 
be realized experimentally and the vortex physics was observed recently on 
mesoscopic scales (BoUe et al. 1999). 

In such a planar geometry vortices are not strictly confined within the 
film. In principle they may leave and enter the film on a side. However, such 
events, as well as a crossing of vortex lines within the film, always cost energy. 
In our theoretical analysis we will neglect such events. Then, by construction, 
vortex lines leave the superconducting plane only at its boundary, the vortex 
lattice does not exhibit topological defects like dislocations and, hence, it 
cannot melt unless we reach Tc. Thermal fiuctuations lead therefore to a 
quasi-long-range ordered lattice as in other 2D crystals (Mermin and Wagner 
1966) which persists roughly up to the mean-field transition temperature T^q. 
In the following we consider the influence of frozen-in disorder on this system. 
It turns out that the problem can be mapped onto a magnetic model, the XY 
model in a random field, for which many results are known. In addition to the 
phase with quasi-long-range order there is a low-temperature glassy phase, 
which is the result of the balance between thermal and disorder fiuctuations 
(in this sense, the situation is different from that for bulk materials, where 
thermal fluctuations have a much weaker effect on the formation of the glassy 
state). Drift of the vortices due to an external current (perpendicular to the 
film plane) leads in the glassy phase to a vanishing linear resistivity. 

4.1 Mapping on the XY model in a random field 

The mapping of this model onto a two-dimensional XY model in a random 
field was first considered by Fisher (1989) in his seminal paper on vortex 
glasses. Here we will follow a slightly different derivation. 

To formulate the model, we assume that the superconducting film is in 
the {x, z) plane and that the vortex lines are directed along the z axis with 
an average line spacing a. The position of the nth vortex line is denoted 
by Xn{z) = Xn + Un{z), whcrc Xn = na and Un{z) = u{Xn,z) are the rest 
position and the displacement, respectively. The Hamiltonian can then be 
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written in the form H. — Hei + "Wpin, where 

= ld'r{'f{d^ur+^f{d,ur}, (4.1a) 
n^in = / d^r Pu{r)V{r). (4.1b) 



Here r = {x,z), cn and C44 are the compression and tilt elastic constant, 
p„(r) the vortex line density and V^(r) the random pinning potential. For the 
latter we assume for simplicity 



V{r) = 0, (4.2a) 



V{r)V{r') = A{x - x')S{z - z') (4.2b) 

and, as before, A(a;) ki Ao/("\/27r^)e~^^/^^^, but Aq ~ npVpin^^*) where np^ 
and /pin denote the impurity density per unit area and the individual force 
of a single impurity on the vortex line, respectively (Blatter et al. 1994). 

In general, the elastic Hamiltonian will be non-local, reflecting the disper- 
sion of the elastic constants for wave vectors |k| > A~^. Since we, however, 
are interested in the disordered case, which is dominated by large scale fluctu- 
ations, wc ignore this fact here. This approximation is justified in particular 
for weak disorder, a case in which all relevant length scales are much larger 
than A. Thermal and disorder fluctuations will considerably renormalize cn 
and C44 with respect to their mean-field values. 

If one starts with isolated vortex lines of displacements w„ which have 
a stiffness constant and a mutual short-ranged interaction energy U{a + 
Un+i - Un), then 

Cii = all" (a), C44 = (4-3) 

On scales a ^ A mean-field theory yields U{a) = UMpio) ~ e""/'*'. 

Thermal fluctuations on scales |Aa;| < a lead to a steric repulsion (Po- 
krovsky and Talapov 1979) between the vortex lines, which reads for A — > 

rp2 2 

UM - (4.4) 

such that {cnCuy/'^/T — n/a^ becomes independent of temperature. As we 
will see below [cf. (4.25)], this relation is exactly the condition for the glass 
transition temperature Tg = (cuCuY^'^a^ /n. In other words, the relation 
(4.4) maps the system at a// temperatures to the glass transition temperature 
Tg. However, in deriving (4.4) only a hard-core repulsive interaction between 
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vortex lines was assumed (indeed, Pokrovsky and Talapov (1979) mapped 
their problem onto non-interacting fcrmions). If we include an additional 
interaction energy of the order of Uc per contact point of two fluctuating 
vortex lines, we get 



This leads to Tg = {cuC^^.f/'^a^/Ti = T(l + Ujn'^Tf/'^, which corresponds to 
T ^ Tg for C/e ^ 0. 

Disorder leads also to a stcric repulsion (Kardar and Nelson 1985, Nat- 
termann and Lipowsky 1988) of the form (we drop here coefficients of order 
unity because they cannot be determined very accurately) 

t/pin(a) ^ Tf;^-. (4.6) 
where Ta ~ (£||^Ao)^'^^ (Nattermann and Renz 1988). This gives 

The vortex-hne distance 

(4.8) 

denotes a cross-over length scale from a thermal steric repulsion (for a < Oc) 
to a disorder dominated steric repulsion (for a > Cc) . The formulas for [/pin 
and Qc are not exact l)ut represcut crude interpolation formulas between the 
limiting cases T <^ Ta and T ^ Ta. 

Finally, a contact interaction can also be included in [/pin(a) such that 



Here we have written c^i{o) in a form similar to c^i{[a). Comparing these 
two quantities, it is easy to see that c^i{a) > cfl{a) in the region a > 
and hence we are below Tg (provided [4 > 0) . 

One should, however, take into account that all these expressions were 
derived under the assumption that the bare interaction between vortex lines 
is short ranged, i.e., for the case A <^ a. In the opposite case the calculation 
is more involved and will not be considered here further. 
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Next we consider the pinning energy (4.1b). Using the Poisson- 
summation formula, wc may rewrite the density Pu(i') ^ (Nattermann 1990, 
Nattermann et al. 1991; see also appendix A) 

oo 
n=—oo 

= ^ / dX S{X-X^)S{x-X-Ur,{z)) 

n=-oo '^-'^ 
1 /"fo 

^ - J2 dX e'^'-''S{x-X -u{X,z)) 

° m=-oo'^-°° 

= - / dX {1 + 2Y^ cos{QmX))5{x -X- u{X, z)) 

" -^-^ m>l 

f« ^|l-9^^x(r) + 2 J]cos(g^[x-ii(r)])}. (4.10) 



m>l 



= 27rm/a is a reciprocal lattice vectors of the 2D line array. The pinning 
energy (4.1b) can hence be written as 

+ J^[^im(r) cos(mv?(r)) + /i2m(r) sin(mv9(r))] j, (4.11a) 

m>l 

2 2 

^im(r) = -l^(r) cos(Q^a;), /i2m(r) = -V^(r) sin(Q„a;), (4.11b) 
d a 

where we introduced the phase field ip = 2nu/a [which should not be con- 
fused with the phase (j) of the superconducting order parameter; see equation 
(4.38) below]. The disorder fields ham {ot^P = 1,2) are Gaussian distributed 
with 



hnJr) = 0, (4.12a) 

ham{r)hpn{r') ~ 2-^fam,l3niQmX,QnX')S^{r -r'). (4.12b) 

a 

On scales Ax ^ a, fam,i3n{QmX, Qn^') ~ ^apSran siucc the rapidly fluctuating 
contributions with m ^ n or a ^ (3 average to zero. On these scales, Ttpin 
can also be written in the form 

n^in ~ fd^r[- ^a,(^(r)y(r) + ^ -A^' cos(m(^(r) - a„(r))}.(4.13) 

^ m>l 
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The random phases ^^(r) obey 



0, 



(4.14a) 
(4.14b) 



gi[am(r)-a„(r')] _ 



The effect of fluctuations of individual vortex hues on scales Ax < a was 
discussed briefly in the first part of this section. Since we study the film in a 
parallel field mainly because of the possibility of obtaining explicit results for 
the vortex-glass state in a case which we understand to a large extent (and 
not so much because of its experimental relevance), we will ignore collective 
fluctuation effects on small and intermediate scales but concentrate directly 
on the large-scale behaviour where typical phase fluctuations are of the order 
27r. On these scales terms with m > 1 in (4.13) can be ignored since they 
are less relevant than that with m = 1. 
It is convenient to rewrite (4.1a) as 



by rescaling the z coordinate and introducing the stiffness constant J accord- 
ing to 



This rescaling leaves Tipin form- invariant apart from a multiplication with a 
factor (c44/cii)^/''. Thus Kei + ^pin now equals the Hamiltonian of a random- 
field (RF) XY model without vortices, which was considered in the context 
of magnetic models by a number of authors (Cardy and Ostlund 1982, Gold- 
schmidt and Houghton 1982, Goldschmidt and Schaub 1984, Villain and 
Fernandez 1984). 

After rescaling (4.16a), Tipm can be written as 




(4.15) 




(4.16b) 



(4.16a) 




(4.17) 



with 



V{ip,r)V{^',r') 



y = o, 

cos{ip — ip')6^{T — r'). 



(4.18a) 
(4.18b) 

(4.18c) 



57 



where we introduced the variances a and /i^. Their bare (unrenormahzed) 
values are given by 




1/2 



(4.19a) 



(7o = 



87rJ2' 



(4.19b) 



This relation between a and Ir? in (4.19) holds only for the unrenormahzed 
quantities. However, under renormalization these quantities flow indepen- 
dently. Although (4.13) originally contained only a linear gradient term in 
the X direction, we rewrote it here in an isotropic form, anticipating the 
generation of dz'^ terms under the renormalization-group transformation. 

Before we come to the renormalization-group calculation for the model 
(4.17), we will apply a Flory-type analysis to it. To be more general, we 
will present this analysis in d dimensions. One complication with respect 
to the problem considered in section 3.1.1c is the fact that the correlator of 
the pinning Hamiltonian now contains the oscillatory term (4.18c). For weak 
disorder /la -C J, a given configuration of the random phase ai(r) = a{r) 
implies a certain ground-state configuration (po{r) = v5o(i", {«}) which de- 
pends on the value of the random phase everywhere in the system. Go- 
ing through the different configurations of {a} will generate a distribution 
of ground state configurations {v^o} which we assume to be Gaussian dis- 
tributed. Since ha -C J, the correlation of the value of a(r) and (po{r) at the 
same position r will be weak. In averaging over the distribution of a{r) we 
can therefore neglect the local correlations between a{r) and </^o(r) in a first 
approximation. The typical energy gain from the pinning term then becomes 



Depending on the dimension d, the leading contribution to E'pin can be the 
random-tilt or the random-field contribution. The dependence of ip'^ on the 
system size L follows from equating the absolute value of Epi^ to the averaged 
elastic energy 



For 2 < d < 4 one finds that the random-field contribution dominates in 
-Spin, which results in 




f« — max 



{ JV^L('^-2)/2(^)V2^ ^^d/2g-,,V2| 



(4.20) 



E^i JL'^-'^cp^. 



(4.21) 



(^2 ~ ^4 _ In 



L 



21n[(^2]^ 



(4.22) 
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where we introduced the Larkin length ^ ( J/Zi)^/^^"''). Thus, the Flory 
argument gives at T = a logarithmic increase of (/^^ in all 2 < (i < 4. 
We will see in section 6 that a renormalization-group calculation confirms 
essentially this result, apart from a (small) modification of the coefficient 
of the logarithm. For d < 2 on the other hand, Epi^ is dominated by the 
random-tilt contribution and the balance of Epi^ and Eei ai T — results in 
aLP'~'^ in agreement with renormalization-group calculations (Villain 
and Fernandez f984). We will show below that in ci = 2 dimensions a is 
renormalized to large values such that Ep^^ and ip^ are dominated by the 
random-tilt contribution. 



4.2 Renormalization 

After replication, we get from (4.17) and (4.18) 



= E / - ^)^^" • -^^^^^^ - /)} • (4-23) 

ah ^ ^ 

The RG flow equations of this model, first found by Cardy and Ostlund 
;i982), are then 



dJ 


= 0, 




H 


da 




H 




dh 








H 





with the glass transition temperature 



Tg = 471 J. 



(4.24a) 
(4.24b) 
(4.24c) 

(4.25) 



The integration of (4.24c) yields 







1 













(4.26) 



with / = ln(L/a), the fixed-point value h*"^ = T^rg/c2 and the reduced tem- 
perature 



-1 ^ 

-L CF 



(4.27) 
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h 




g 



Figure 7: RG flow in the {h, r)-plane. The random-field amplitude h 

flows to h* = above the glass transition temperature (T > Tg) whereas 
it flows to h* oc (Tg - T) below (T < Tg). 



It can be shown quite generally that the flow of a does not feed back into that 
of h (Hwa and Fisher 1994b). The constants Ci and C2 are cut-off dependent, 
but their ratio 

| = ^ + OW (4.28) 

is universal. The whole calculation is valid only for Tg <^ 1. Considering 
equivalent models, equivalent RG equations were found by Goldschmidt and 
Houghton (1982), Rubinstein et al. (1983), Goldschmidt and Schaub (1984), 
and Paczuski and Kardar (1991). The parameter flow is shown schematically 
in flgure 7. 

The quantity h"^ approaches a fixed point h*^ T^Tg/c2. Integrating 
(4.24b) for /i /i*, we find on the length scale L — ae'- 

a{L) ^ 47rT^\n{L/a). (4.29) 

If we go back to the renormalized but unrescaled quantities - these are the 
physical quantities we would measure in an experiment on this scale - we 
obtain for the effective parameters on scale L 

hes ~ h{L)^, (4.30a) 
(Teff = a(L). (4.30b) 

Thus, while the effective random-field strength vanishes as L~^, the variance 
of the coefficient of the linear gradient term grows logarithmically with length 
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scale L. As a side remark we mention here that the effective value of the 
coefficient of the non-linear term in 7Y„ scales as L^^ in agreement with the 
general result for A^s ~ h'^g ~ L'^^^, cf. equation (6.28) below. 

It is interesting that Villain and Fernandez (1984) considered this problem 
at T = in a complementary study and found the following set of RG flow 
equations (we rewrite here their discrete recursion relations in a differential 
form) : 

with some numerical constants Ci, C2. h flows now to inflnity, h{L) ^ 



2/c2JlnL, such that da/dl ^ Ci/c2. Thus both calculations, for T < Tg 
and for T = 0, give a logarithmic increase of a^siL) and a vanishing heg{L) 
for L — * oo. This gives further credibility to the existence of a unique low- 
temperature phase. 

A central role for the characterization of this phase plays the displacement 
correlation function (Goldschmidt and Houghton 1982) 

T , r „ o , 9 r 



W{r) = {[ip{r) - <^(0)]2) = ^ In - + 2t^ In^ -. (4.32) 

Quahtatively, the same result was found by Villain and Fernandez (1984) 
at T = 0. From (4.32) it is possible to calculate the correlation function 
of the order parameter for translational order ■0fc(i") — e**^"'-'"-' [cf. equation 
(3.8)]. For k = Qm — 27rm/a one has ku — rrup and hence (Goldschmidt and 
Houghton 1982) 



S{Qm,r) = ^e*<3mWr)-«(0)]^ ^ ^gim[¥'(r)-<^(0)]\ 



where 



Thus, in the glassy phase, at T < Tg, the correlation function S{Q, r) decays 
slightly faster than with a power law. However, the In^ r-behaviour dominates 
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W{r) only on scales |r| ^ Lj-g ~ ae^/'^s which is large within the range of 
vahdity Tg <^ 1 of the RG equations (4.24). The structure factor 




behaves for small q = {qx, \/ c/^j c\\qz) near the first reciprocal lattice vector 
Q — ±2n/a, therefore, as 



S{Q + qx,qz 



) ~ q 



2Tg[l-rgln(ag)] 



(4.36) 



Contrary to the pure system, where quasi-long-range order is accompanied by 
algebraic Bragg peaks, these arc smeared out here for q < L~\ For T > Tg, 
the term oc Inr in (4.34) is absent and the power law decay of correlations is 
regained. 

For completeness, and since the pair correlation function S{Q, r) vanishes 
faster than a power law, we also consider the positional glass correlation 
function S'pG(Q,r)- cf. equation (2.33). This was also calculated already by 
Goldschmidt and Houghton (1982), who found 



for |rg| ^ 1 on both sides of the transition. Note that to the given accuracy of 
order Tg this glass correlation function decays as a power law in both phases. 
In the high-temperature phase, where disorder is essentially irrelevant, the 
lack of a mechanism for melting in the model - in particular the absence of 
dislocations - prevents S'pg from vanishing. In the low-temperature glassy 
phase, on the other hand, thermal fluctuations are still too strong to allow for 
true long-range order of the glass correlation function in this two-dimensional 
system. Note, however, that the relative magnitude of the correlation func- 
tions is qualitatively different in both phases, since Spq/S'^ — > const, for 
T > Tg but Spg/S^ ^ oo for T < Tg and |r| oo. 

We flnally remark here that for this simple model, which has no shear 
modes, also the phase-coherent vortex glass correlation CvG(r) (2.34) shows 
an algebraic decay. Indeed, since in this layered geometry the phase (p of the 
superconducting order parameter changes by tt at a vortex line, the latter is 
related to the displacement by 



^PG(Q,r) 



~ r 



-4T/Tg 



(4.37) 




(4.38) 



With (2.34) and (4.37) we find therefore 




T/T, 



(4.39) 
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i.e., the system shows quasi-long-range vortex-glass order. 

Our result (4.32) for the glassy phase is in contradiction with a number of 
results by other authors using Bethe Ansatz (Tsvelik 1992, Balents and Kar- 
dar 1993), a variational treatment with one-step replica-symmetry breaking 
(Korshunov 1993, Giamarchi and Le Doussal 1994), and variational methods 
without rephcas (Orland and Shapir 1995). In these studies W{r) oc Inr 
was also found for the low-temperature phase, but with a temperature in- 
dependent coefficient for low T. We believe that these results are incorrect 
and in particular that they demonstrate the flaws of the Gaussian variational 
method, which gives only Flory-like results for the correlation function. In- 
deed, Bauer and Bernard (1996) showed for an A'"-component version of our 
model (4.23), that the coefficient of the In^r-term in W{r) vanishes as 
for large N. This explains the absence of the In^ r-term in the variational 
calculations which give exact results only in the limit N ^ oo. 

Numerical studies give a controversial picture: while in the investigations 
of Batrouni and Hwa (1994) and Cule and Shapir (1995) only a In r-behaviour 
of l^(r) was found, later studies (Lancaster and Ruiz-Lorenzo 1995, Marinari 
et al. 1995, Zeng et al. 1996, Rieger and Blasum 1997) were able to detect a 
In^ r-behaviour. However, the coefficient of the In^ r term found in the finite- 
temperature simulations is much smaller than the RG prediction (4.32), even 
if one takes into account that the true coefficient 2t^ is smaller by a factor | 
than assumed originally. The difficulty in observing the latter behaviour may 
partially be explained by a large crossover length. More recently, Zeng et al. 
(1999b) were able to con&m in their simulation both the In^r dependence 
and the Tg prefactor of the correlation W{r). 

The faster decay of correlations of the order parameter for translational 
order in the low-temperature phase clearly indicates the influence of disorder. 
However, the resulting phase is not necessarily a glassy phase. Indeed, as was 
remarked by Hwa and Fisher (1994b) and mentioned in connection with our 
discussion of manifolds, a harmonic toy model (4.15), (4.17) with l^((^, r) = 
but [compare (3.34)] 

//,(r)//,(rO = 5,,JV/(r-r'), /(r) ~ |r|-", a<d (4.40) 

already results in 

W{v) = ([(^(r) - (^(0)]2) ~ (7|r|2-". (4.41) 

Although the roughness exponent C, — \{2 — a) can be larger than zero in 
d = 2 dimensions, and hence S{Q, r) would decay like a stretched exponential 
function, the system is certainly not glassy since it has indeed a continuum 
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of ground states which differ by the value of a constant </7i but have the same 
energy. Thus, there is no pinning of the vortex-hne lattice in this toy model. 

We conclude from this exercise that a stronger decay of the correlation 
function S{Q, r) than that resulting from thermal fluctuations is a necessary 
but not sufficient requirement for the existence of a glassy phase. It is clear 
that the existence of many meta-stable states and, hence, the anharmonicity 
of the model is decisive. 

4.3 Susceptibility 

As a better signature of glassiness Hwa and Fisher (1994b) proposed to look 

at the response to a change SH = (SHx.SHz) in the applied field. This 
change leads to a new term — ($o/87r^) J (Pr SH ■ in the Hamiltonian. 
One can now consider the susceptibility 

X^J ^ (4-42) 
In isotropic systems Xij — X^ij the magnetic permeability is found as 

/^mag-X$o/(167r')- (4-43) 

It was then shown that the sample-to-sample fluctuations of x or A*mag in the 
low-temperature phase fulfill the relation 

{A^^{x^xr = CT,x\ {T<T,), (4.44) 

where x = 1/J is independent of the disorder as a result of an underlying 
statistical tilt symmetry (Schulz et al. 1988). C is a universal boundary- and 
geometry-dependent coefficient and Tg ~ C2(/i*/Tg)^ [cf. equation (4.27)] is 
a measure of the effective non-hnearity of the model on large length scales. 
On the contrary, (Ax)^ ^ for T > Tg and for large systems. The non- 
vanishing sample-to-sample fluctuations of the susceptibility for T < Tg, 
which may be tested already at a simple sample by rotating the direction of 
the magnetic field, are therefore a better way to define glassiness. Although 
early numerical simulations by Batrouni and Hwa (1994) were not able to 
confirm prediction (4.44), more recent simulations by Zeng et al. (1999b) are 
in agreement with (4.44). 

For completeness we mention in addition a study by Sudb0 (1993), in 
which he considers the response D ~ Xxx of the vortex-line array to a tilt 
of the magnetic field. By mapping the system of hard-core repelling vortex 
lines onto non-interacting ID fermions in a time dependent random potential. 
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he identifies a possible glassy phase of vortex lines with the localization of 
fermions. For electrons in a time independent random potential one expects 
a vanishing average charge stiffness D in the localized phase, in contrast to 
his Bethe-Ansatz calculation, which gives D ^ 0. However, as mentioned 
before, D 7^ is the consequence of an underlying statistical tilt symmetry 
(which exists only for a fermion potential which is random both in space and 
in time direction) and not a signature of the absence of a glassy phase (Hwa 
and Nattermann 1993). 



4.4 Dynamics 

The physically most convincing way to demonstrate glassy properties of a 
system is to consider the dynamics. Before we begin this topic, we will make 
a small digression to discuss the appearance of electrical resistivity from the 
motion of vortex lines. Under the influence of an external current density j 
a Lorentz force with density 

F = -j A B (4.45) 

c 

acts on the vortex-line array. B denotes its magnetic induction and c the 
velocity of hght. The force density leads under dissipative conditions to a 
steady-state motion of vortex hues of velocity v = v(F), which generates an 
electric field E = BAv/c. Ifv and F arc parallel to each other, E and j are 
as well. The resistivity of the vortex-line array follows, therefore, from 

A simple situation exists if the relation v(F) —F/rj is hnear with a friction 
coefficient rj, for which p = B'^/{c'^r)) holds. An example for a linear relation 
is the Bardeen- Stephen flux flow where rj — BHc2/{c^Pn) with the normal 
resistivity p„ (Bardeen and Stephen 1965). 

We consider now the influence of disorder in the 2D case with the magnetic 
field parallel to the film plane. The equation of motion of the vortex- line array 
under the influence of a driving force density F in x direction (the external 
current is assumed to be perpendicular to the film) reads 

W = + / + (4-47a) 

(C(r,i)C(r',0> = 2rj^T5{v-v')5{t-t'), (4.47b) 
with f ^ F and 770 ~ BHc2/pn- 
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Prom a dynamical RG one finds (Goldschmidt and Schaub 1985, Tsai 
and Shapir 1992, Tsai and Shapir 1994) a renormalization of the friction 
coefficient according to 



dr] 



"g 



(4.48) 



wliere Cg = e^'''^"/4 and 7eu is Euler's constant, i.e., 2Jc^ ~ 1.78. C2 was in- 



troduced in equation (4.24c). The integration of (4.48) gives for the effective 
friction constant on the scale L — ae' 



r]{L) f» 77o 



1 + 



C2 / ho 



V^g 



{{L/ar^ - 1) 



-z-2 



(4.49) 



for T < Tg, where we used (4.24c). Because of the absence of a trivial 
rescaling in (4.48), r]{L) coincides with the effective friction coefficient r]es{L). 
The relation 77 (L) ~ L^~^ defines the dynamical critical exponent z with 



2 + 4^rg for T < Tg 
2 for T > Tg. 



from (4.49). 

For T > Tg, Tg < and the asymptotic (L — > 00) value ?7oo of V is 



(4.50) 



(4.51) 



Prom this one finds for the linear resistivity close to T„ 



B 



Pn 



H. 



c2 



1 + C2;^|rg 
g 



|2v^ 



(4.52) 



In other words, we obtain for the depinned solid phase a linear (Ohmic) 
resistivity which vanishes at Tg as a power law in |Tg| with the exponent 
2^. 

For T < Tg, Tg > and 1]{L) increases with increasing length scale until 
L ~ Tmax ~ a{J/Fy^'^ is reached where the RG flow is stopped (Tsai and 
Shapir 1992). The response is therefore non-hnear with (F <S J) 



77eff(Tg > 0,F) ^ r]Q 



1 + 



C2 / k 



g \-^g 



{{J/FP - 1) 



2,/c 



(4.53) 
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which gives a non-linear resistivity 



PeflfU) ~ Pn 



c2 



1+^ ^1 -1 



(4.54) 



where has the meaning of a zero-temperature critical current density, C3 = 
0(1) is a constant and j jc- Apparently, the linear resistivity Pes{j — >■ 0) 
vanishes, i.e., the system in the glassy phase is a true superconductor. 
For Tg \ we obtain in particular at the vortex-glass transition 



-He 



c2 



In^ 

. j 



(4.55) 



i.e., the system is a true superconductor also at Tg. 

The transition to the high-temperature phase is therefore continuous at 
Tg where p(0) = 0. The results of the RG calculation were derived under the 
assumption |Tg| <^ 1 and no definite conclusions can be drawn for T <^ Tg. 
However, it is very likely that the dynamical critical exponent z diverges as T 
goes to zero. The increase of z with decreasing T was confirmed numerically 
by Lancaster and Ruiz-Lorenzo (1995). 

It is interesting to remark that one can obtain the current-voltage re- 
lation ~ £■ ~ j^/^ also from an argument which is a modification of a 
scaling argument first used by Fisher et al. (1991a). Fisher et al. (1991a) 
assume that the vortex-glass transition is continuous and characterized by 
the divergence of the characteristic length and time scales, and t^^ ~ 

respectively. Their original argument is as follows: We express E as — A/c. 
Since in the GL-Hamiltonian A appears in the combination (V — (e/c) A), A 
should then scale at the glass transition as an inverse length and hence 
E as A ~ ig^^^g ~ ^g^^~^- critical properties are described by a finite- 

temperature fixed point, the free energy density j ■ A scales as ^g"'' and hence 
j ~ ^g'^^^- The current- voltage relation has, therefore, in the critical region 
the form 

Eeg^'^^dijii-')- (4.56) 
Since ^g has to drop out of this relation for T — > Tg, we obtain $d(ii) ~ 

E ~ (4.57) 

at Tg, which is the result of Fisher et al. (1991a). 

In the present case of a two-dimensional system with B parallel to the 
plane this argument has to be modified. Since only the z component — 
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dxAy — dyAx of the magnetic field is non-zero and because of the absence 
of any y dependence in this geometry, the relevant component of A is Ay. 
In particular, there is no {dy — ^Ay) term in Hgl and hence we conclude 
Ay E J ~ ^ which results in 

EC^^Mje,)- (4.58) 

At the transition, this gives E ~ J^^^, in agreement with the RG result 
E i (apart from a logarithmic correction) at the vortex glass transition 
temperature, where z — 2. 

In the high-temperature phase (rg < 0) where $2 ~ J^g and 2; = 2, we 
get E j, such that the correlation length drops out of the result. The 
same is true for T < Tg, where the RG result gives E ~ J^^^. From this one 
has to conclude that both phases are critical such that the transition is not 
accompanied by a diverging correlation length. 

For completeness, we remark here that Toner (1991a) and Nattermann 
et al. (1991) tried to calculate the current-voltage relation from an esti- 
mate of the energy barrier EbU) between different mcta-stablc states. Both 
authors use the result from the statics for ct(L), equation (4.29). While 
Toner (1991a) finds EbU) ~ rg[ln(JVi)]^^^ Nattermann et al. (1991) get 
EbU) ~ ''■g lii(ic/i)- Together with the general creep formula (3.73), the 
latter result gives a power law for the current-voltage relation, but with a 
wrong result for the exponent z — 2, which is proportional to Tg and not to 
Tg. This discrepancy is related to the fact that for T < Tg the physics is 
described by a finite-temperature fixed point at which energy barriers are 
not well defined. 

We conclude from this section that below the temperature Tg the vortex- 
line array in an impure superconducting film in a parallel field exhibits a 
glassy phase with vanishing linear resistivity. This phase is the most ac- 
curately studied example of a vortex glass (although of mainly academic 
interest). The spatial correlations of the vortex-line array exhibit a decay 
slightly faster than algebraic. The vortex-glass transition is not accompa- 
nied by a diverging correlation length (as assumed in the scaling argument 
of Fisher et al. (1991a) and as seen in many experiments for bulk materials) 
since both phases are critical. Because of the absence of shear modes, below 
Tg the system is both a positional and a phase-coherent vortex glass (in the 
sense discussed in section 2.4). 
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5 Superconducting film in a perpendicular 



field 



So far we examined in chapter 3 a single vortex as an example of a directed 
manifold and in chapter 4 a periodic array of vortex lines in a film. Both cases 
are simple in the sense that the vortex conformations can be described by 
a single- valued displacement field, or in other terms, that these systems are 
free of topological defects. The possible presence of such defects implies that 
the elastic approximation, which was the basis of our analysis in the previous 
chapters, may break down. We now turn to the simplest system that allows 
for topological defects: an array of point-vortices in a superconducting film 
induced by a perpendicular magnetic field. 

We consider here only the case of a thin film, where the film thickness s is 
much smaller than the bulk correlation length ^. Then the condensate wave 
function can be considered as constant in the direction perpendicular to the 
film, i.e., it is a two-dimensional (2D) degree of freedom. Nevertheless, the 
screening mechanism remains three-dimensional, since the electromagnetic 
field still propagates into the third dimension. 

5.1 Pure system 

To set the stage for an analysis of pinning effects in this two-dimensional 
system, we briefiy summarize important properties of the pure system (see in 
particular Doniach and Hubermann (1979), Hubermann and Doniach (1979), 
Fisher (1980)). 

The structure and interaction of vortices in thin films were first studied 
by Pearl (1964, 1965). Several distinct features are to be mentioned, (i) 
Because of screening, the creation of a vortex costs a finite amount of energy. 
This is in contrast to vortices in 2D superfluids, the energy of which diverges 
logarithmically with the system size due to the absence of screening, (ii) 
Every vortex has a magnetic moment M that diverges proportional to the 
system size L. Thus the lower critical field vanishes in infinite films, Hd oc 
1/L. (iii) The vortex pair interaction potential is 




(5.1a) 




x < 1 
X > 1 



(5.1b) 



where 




(5.2) 
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is the energy scale and A is the bulk penetration length. A further important 
length scale is the screening length 



L, = 2AVs, (5.3) 

which may be macroscopically large (in the range of millimeters). Screening 
sets in only on scales beyond Lg, i.e., this scale is intimately related to the 
inhomogeneity of B. 

The importance of the screening length becomes evident already at zero 
magnetic field {H — 0). Only in the hmiting case Lg ^ oo the vortex self 
energy diverges logarithmically with the system size and individual vortices 
cannot be created by thermal activation. Then vortex- ant ivort ex pairs inter- 
act logarithmically on large scales and dissociate according to the entropy- 
driven Kosterhtz-Thouless (KT) mechanism (Kosterhtz and Thouless 1973) 
at a temperature 

Tkt - ^ = (5.4) 

In this expression we neglect renormalization effects which lead only to a 
quantitative reduction of the transition temperature. However, in the limit 
Lg — > oo the transition temperature vanishes like Tkt oc l/Lg and the system 
becomes critical only at zero temperature. At large but finite Lg a crossover 
in the vicinity of the small but finite Tkt survives. 

The electric resistivity of the film can be understood in terms of vortex- 
antivortex pairs under the influence of a transport current which drives vor- 
tices and antivortices in opposite directions (Halperin and Nelson 1979). For 
finite Lg and T > the film has Ohmic resistivity at sufficiently small cur- 
rent densities, since vortex-antivortex pairs have to overcome only a finite 
energy barrier to dissociate. Thermal activation therefore leads to a finite 
dissociation rate and a resistivity proportional to the driving current. At 
larger current densities, where the vortex interaction is probed on scales 
r <S Lg (i.e., on this scale the attractive vortex interaction is balanced by 
the driving force), the logarithmic vortex interaction implies a power law for 
the current- voltage relation, which was predicted (Halperin and Nelson 1979) 
and observed experimentally (Epstein et al. 1981, Kadin et al. 1983) already 
long ago. More recent experimental (Repaci et al. 1996), numerical (Simkin 
and Kosterhtz 1997), and analytical (Pierson and Vails 1999) work was de- 
voted to examining finite-size effects on the KT transition, which includes 
the case of a superconducting film regarding Lg as 'intrinsic' finite-size scale. 

Although the case of zero field would be of interest on its own, we sub- 
sequently consider only finite fields, H > 0, which induce a finite average 
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vortex density. Then one can think of the vortex system as a superposition 
of a neutral subsystem composed of vortex dipoles and a subsystem of excess 
vortices oriented in the direction parallel to H. 

We now focus on the subsystem of excess vortices and assume that the 
effect of the dipole subsystem amounts to a screening of the interaction be- 
tween the excess vortices. This interaction is expected to be purely repulsive 
even in the presence of screening by dipoles such that the ground state of 
the excess vortices should be a triangular lattice. Thermal fluctuations lead 
to the creation and annihilation of extra excess vortices and thereby create 
vacancies and interstitials in the lattice. Such fluctuations cost finite mini- 
mum energy in contrast to vortex displacements. Therefore we restrict the 
subsequent analysis to displacement fluctuations and examine their influence 
on the positional and orientational order of the lattice. 

In the literature the vortex lattice is, in general, not examined on the 
basis of the physically given vortex interaction (5.1). Instead, simplified 
model interactions are used. Since certain physical features depend on the 
choice of the model interaction, we attempt to classify them by distinguishing 
the following cases according to the behaviour of the potential on large length 
scales: 

'Short-range' case: U is short-ranged, such that J (Px f/(x) < oo. Only 
in this case the interaction energy of the vortex lattice is extensive and 
all elastic constants are finite. 

'Long-range' case: The interaction U decays on large length scales but it 
is long-ranged, such that J (fx C/(x) = oo and the interaction energy 
is not extensive. This case includes the Pearl interaction. The system 
is incompressible, i.e., it has a sm^'uZar compression modulus. 

'Logarithmic' case: The interaction is logarithmic on large length scales, 
C/(x) oc Inx. This case, which also represents an incompressible system, 
includes the so-called two-dimensional Coulomb gas and the supercon- 
ducting film in the limiting case when is larger than the system 
size. 

We now briefly summarize the elastic properties of the system and the pres- 
ence of phase transitions in these cases. 

Short-range case. Among these three cases, this one has been examined 
most since it can be represented by a well-defined elasticity theory with 
finite elastic constants. In general, the elastic Hamiltonian can be written in 
Fourier space as 

^ei[u] = 11^^ ^[cn(q)g>L(q)r + C66(q)g>T(q)r]. (5.5) 
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The displacement field is defined in position space only as a function of the 
discrete positions X of the undistorted lattice. In Fourier space it reads 

u(q) = a'^e-^'i-^ u(X) (5.6) 

X 

with the area = ^o/B per vortex. The displacement can be decomposed 
into a longitudinal and a transverse component 

Ui(q) = P^(q)-u(q), UT(q) = P^(q) • u(q) (5.7) 

by means of the projectors P^/^iq) = ga^/j/g^ and Pj^(q) = 6^,^ - P^f^iq). 
The energy of longitudinal and transverse displacements is given by the com- 
pression modulus Cii and the shear modulus Cqq, respectively. In general both 
moduh depend on the wave vector q because of the non-locahty of the inter- 
action. 

The clastic moduli are known explicitly for the special interaction 
U^'^\x) = Kq{x) with the Bessel function Kq, which decays exponentially be- 
yond the screening length. This is the interaction of straight vortex lines in a 
three-dimensional superconductor apart from a replacement of the screening 
length of the bulk material, A, by the screening length of the film, Lg. Thus 
one finds from the bulk moduli (Blatter et al. 1994) within the continuum 
isotropic approximation 

C66(q) = 2s£o^, (5.8a) 

5<i>0 

c„(q) = (5.8b) 

In this chapter we use B to parameterize the vortex density S/0o per unit 
area rather than as average magnetic induction. 

In the case with finite elastic constants, the pure system can have two 
distinct phase transitions (for a review, see e.g. Strandburg (1987)): a melt- 
ing transition between a solid and a hexatic liquid phase, and a transition 
between the hexatic and an isotropic liquid (Kosterlitz and Thoulcss 1973, 
Halperin and Nelson 1978, Nelson 1978, Young 1979, Nelson and Halperin 
1979). In this scenario, called the KTHNY scenario after the aforementioned 
authors, the solid has quasi-long-range positional order and long-range ori- 
entational order; the hexatic liquid has short-range positional order but still 
quasi-long-range orientational order, and finally in the isotropic liquid orien- 
tational order is also short-ranged. These transitions have been discussed in 
the particular context of superconducting films by Doniach and Hubermann 
(1979), Hubermann and Doniach (1979) and Fisher (1980). 
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The melting transition is driven by the entropic unbinding of dislocation 
pairs. In a triangular lattice, which is the type of lattice formed by vortices, 
dislocations have a pair interaction 

t/disioc(x) c ln(x/a) (5.9) 

7rV3 cii 

(Young 1979, Fisher 1980), where cu = Cii(q = 0). The given expression 
represents only the leading term for large distances; at smaller distances 
additional terms are present. Thus, according to the Kosterlitz-Thouless 
scenario, dislocation pairs unbind at a temperature 

T ^ 1 ^0 C66(cii - cee) , . 

where we have ignored the renormalization of the elastic constants. The 
quantitatively correct value of Tm is hard to determine analytically, since it 
depends on the dislocation core energy, i.e. the discrete vortex structure of 
the core, which is difficult to include in analytic treatments. 

It is interesting to note that the melting transition can be captured 
by a modified Lindemann criterion. The traditional Lindemann criterion 
(u^(X)) < CLja^ for the stability of the lattice is not applicable, since in 
two dimensions (u^(X)) = oo at any finite temperature. However, consider- 
ing the relative displacement of neighbouring vortices (separated by a basis 
vector b with |b| = a a) instead of the absolute displacement, one obtains 

cl,al < ([u(X)-u(X + b)]^) 
q 

(5.11) 



2 



alT 



where we have assumed cn ^ cqq (Scheidl and Vinokur 1998a). A compari- 
son of (5.10) and (5.11) allows the determination of the Lindemann parameter 



1/2 



This value may be a good reference value to estimate the stability of higher- 
dimensional systems with quasi-long-range positional order (such as the bulk 
vortex lattice with pinning) to the proliferation of topological defects. 
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The second transition from a hexatic liquid to an isotropic liquid is driven 
by the unbinding of disclination pairs. Such pairs also have a logarithmic 
interaction, 

^disci(a;) = ^^^ln(a;/a) (5.13) 

with the Frank constant Ka (Nelson and Halperin 1979). Thus disclination 
pairs unbind according to the KT mechanism at a temperature 

T. - (5.14) 

Again, the quantitative value of Th is very difficult to determine analyti- 
cally. Partially, the difficulty is due to the strong temperature dependence 
of Ka = Ka{T), which diverges above the melting transition. However, this 
divergence ensures that Th > Tm, i.e., that the lattice is more stable against 
disclinations than against dislocations. Thus, in order to study the stability 
of the elastic (solid) phase, it is sufficient to focus on the proliferation of 
dislocations. 

Long-range case. The actual interaction (5.1) between vortices in a 
film is long-ranged. If we ignore the regime on length scales below Lg and 
take U{x) — 2seoLs/x, the interaction is the electrostatic one for a Wigner 
crystal. In this case the elastic moduli are (Bonsall and Maradudin 1977) 

/ ^ \ 3/2 

C66(q) ~ 0.25-2s£o^J — j , (5.15a) 

27r / 5 ^ ^ 



cii(q) = 2seoL,— — ■ (5.15b) 

As in the short-ranged case, the shear modulus is essentially local and, in 
particular, C66(q = 0) is finite. However, cii(q) oc q~^, which means that 
the vortex lattice is incompressible. Thus conventional local elasticity theory 
does not apply. Nevertheless, by formally taking Cn oo, most of the 
results obtained for the short-ranged case can be transferred to this long- 
ranged case. Since the dislocation interaction remains logarithmic on large 
scales one still expects a melting transition at a finite temperature according 
to equation (5.10). 

Logarithmic case. Finally, the logarithmic case is of interest, since it 
describes the superconducting film in the limit Lg oo. In this limit, the 
elastic moduli can be obtained directly from equations (5.8), 

C66(q) = 2seo^, (5.16a) 
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cii(q) = 2seo——^. (5.16b) 

Again, cqq is finite on large scales, but the longer range of the interaction 
has even further increased the incompressibility. Nevertheless, a melting 
transition should be present according to the formal argument used already 
for the long-ranged case. 

The analytical description of the two distinct transitions (solid-to-liquid 
and hexatic-to-isotropic) assumes that the topological defects are very dilute 
at the transition, i.e., that they have a sufficiently large core energy. How- 
ever, for small core energies the KT transition - which is thermodynamically 
of infinite order - may become a first-order transition (FOT) (Minnhagen 
1987, Thijssen and Knops 1988). Thus it is not clear whether the 2D melt- 
ing must have a unique nature or whether there can to be two separate 
transitions for the loss of positional and orientational order. There has been 
controversial evidence for both the KTHNY scenario and a FOT scenario. 
Whereas earlier simulations always seemed to support the FOT scenario, 
it is now believed that this might be an artifact for finite-size effects and 
more recent simulations (Bagchi et al. 1966, Perez-Garrido and Moore 1998, 
Alonso and Fernandez 1999, Jaster 1999) provide evidence for the KTHNY 
scenario, which seems to be particularly well established for particles with 
short-ranged interactions (see, e.g., Chaikin and Lubensky (1995)). 

For the long-ranged case and even more for the logarithmic case, the 
thermodynamic nature of the transition is not unambiguously understood. 
Conventional thermodynamics already runs into trouble since the energy of 
the system is no longer extensive. To close this section on the pure sys- 
tem, we quote some references to provide the reader with some references 
into the literature concerning the long-ranged systems rather than to give a 
comprehensive overview. 

CaUiol et al. (1982) examine the long-ranged case with U{x) <x 1/x and 
find crystallization at T ~ -^2seQ (without being able to resolve the question 
of whether it is a single transition or a sequence of two transitions). This 
value of Tm is smaller than the one given in equation (5.10) by a factor of 
approximately 5, which in principle can be ascribed to a renormalization of 
elastic constants that is not taken into account in equation (5.10). 

Recent simulations (Moore and Perez-Garrido 1999) of the logarithmic 
case suggest that the solid might be destroyed at any finite T. Many simu- 
lations have been performed on the model in the lowest Landau-level (LLL) 
approximation, which also results in a logarithmic vortex interaction due to 
the absence of screening. On the one hand there is evidence suggestive of 
a first-order melting (Hu and MacDonald 1993, Kato and Nagaosa 1993b, 
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Sasik and Stroud 1994b) using the LLL approximation, but on the other 
hand there are also claims of the absence of the solid phase at finite temper- 
atures (O'Neill and Moore 1992, 1993, Yeo and Moore 1996a, 1996b, 1997). 
Moore and Pcrez-Garrido (1999) argue that in the logarithmic case dislo- 
cations screen the interaction of disclinations such that disclination pairs 
unbind and destroy crystalline order at any finite temperature. However, in 
the long-range case numerical simulations are likely to be affected by finite- 
size effects and the choice of boundary conditions even for very large system 
sizes. 

On the experimental side there is evidence for the KTHNY scenario in 
superconducting films (Wordenweber et al. 1986, Berghuis et al. 1990, 

Theunissen et al. 1996). In other systems experimental evidence supports 
the FOT scenario (see Chaikin and Lubensky (1995) and Pcrez-Garrido and 
Moore (1998) for a summary), which partially may be ascribed to the pres- 
ence of additional symmetry-breaking fields. 

5.2 Disordered vortex lattice without dislocations 

We now address the question again of how quenched impurities affect the solid 
vortex lattice. Hereby we focus on the case with short-ranged interactions. 
Then the vortex interactions can be represented by the elastic Hamiltonian 
(5.5), provided displacements are small (i.e., the phase is solid). Since we 
aim mainly at the properties of the system on length scales larger than Ls, 
we may ignore the dispersion of the elastic constants, the value of which 
will be renormalized by fluctuations on small scales, as argued in section 
4. Although the incompressible vortex lattice in a film does not strictly 
belong to the short-ranged case, it should be well captured by the formal 
limit cii — s> oo. 

The effect of disorder on the solid vortex lattice of the film in a perpendic- 
ular field turns out to be very similar to that of the film in a parallel field, as 
treated in the previous chapter. The main complication consists in the fact 
that the displacement field has two (instead of one) components and, conse- 
quently, the different elastic response of shear and compression modes. Our 
presentation of this generalization closely follows Carpentier and Le Doussal 
(1997) and Carpentier (1999). 

According to the lines sketched in Appendix A the effective contribution 
of pinning to the Hamiltonian can be written in the form 

^pin = J]nX + u(X)) 
X 
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= J £x |v(u(x),x) - ^/i«/j(x)[aaM/3(x) + a/3Ma(x)]| (5.17) 

with an effective potential y(u(x),x) = pu(x)V"(x) of zero average and 

y(u,x)y(u',x') A(u-u')5(x-x'), (5.18a) 
A(u) = po5]^(X + u), (5.18b) 

X 

where the effective correlator A is periodic due to a sum over the undistorted 
lattice positions X and A(x) = V{^)V{0) is the correlator of the original 
pinning potential. 

From the correlator of V we have separated a random field /x, the corre- 
lation of which can be written as 

+Crmcl6{SajSf3S + Sas5f3j)]Q^S{x - x'). (5.19) 

For the unrenormalized model a^e — and an — Q\p1/S.q/T'^ arises from the 
coupling of the random potential V to the divergence of the displacement field 
(cf. equation (A.4)). The length Qa of the smallest non- vanishing reciprocal 
lattice vector (RLV) is related to the vortex spacing via Q\ = IGn^/Sa^. 
Equation (5.17) is the generalized analog of equation (4.17) for the film in a 
parallel field. 

Carpentier and Le Doussal (1997) and Carpentier (1999) studied this 
model retaining only the contribution from the smallest non-vanishing RLV 
to the correlator, i.e., approximating A(u) ^ Pq X]q-|q|=Qa ^(Q) cos[Q ■ u]. 

Then the Fourier coefficients A(Q) act like a random-field amplitude, cf. 
equation (4.18c), for which we introduce an effective variance with bare value 
g = poA(Qa)/T^ 

Under renormalization, not only the field amplitude g evolves, but the 
random potential induces also bond disorder, i.e., a flow of an and ctqq. Car- 
pentier and Le Doussal (1997) and Carpentier (1999) found the flow equa- 
tions 



dicu = 0, (5.20a) 

dicee = 0, (5.20b) 

dian = bn{a)g^ (5.20c) 

diaea = beeia)g^, (5.20d) 

dig = 2(^l-^^g-b,{a)g', (5.20e) 
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with the coefficients 
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bn{a) = p^[2Io{a) + h{a)], (5.21a) 

4 

666(«) = ^^^[2/o(a)-/i(a)], (5.21b) 

^ ^66 

b^{a) = 27T[2Io{a/2) - Io{a)], (5.21c) 
that depend through the Bessel functions /o,i on the parameter 

„ . ^ (± _ ±) . (5.22) 

47r \ce6 cuj 

which measures the difference between the lattice response to shear and com- 
pression. Prom the linear term in equation (5.20e) one recognizes that disor- 
der is relevant only at temperatures below the glass transition temperature 
Tg. Its value is given by 

Stt CiiCge V3$o CiiCee 
ig = —5 — — [0.26) 

Qa '^'n +^66 TT S Cll + C66 

(Giamarchi and Le Doussal 1995). 

In principle, there is a finite renormalization of the elastic constants on 
small length scales because of the finiteness of the disorder correlation length 
^. On large scales, where the pinning energy correlator effectively factorizes 
as in equation (5.2), the renormalization of the elastic constants is absent 
due to a statistical symmetry analogous to the tile symmetry discussed in 
section 3.1. Id. The short-scale renormalization is neglected here, since it 
does not influence the large-scale properties of the vortex lattice. However, 
to be accurate the transition temperature is given by equation (5.23) using 
the renorm,aMzed value of the elastic constants. 

The fluctuations of the displacement fleld can be represented by the cor- 
relation function in the form 



= WL{^)P^f,{^) + WT{^)P^f,{^) (5.24) 

with projectors P^p{x) = XaX^s/x^ and Pj^ = 6af3 - Pap- 

Above the glass transition, T > Tg, the renormalization of the elastic 
constants and of the strengths Un and crge of the random 'tilt' flelds are 
flnite such that W{x) ~ ln(a;/a) and the correlation function 5'(Q,x — x') 
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decays algebraically, 



5(Q,x-x') = (e-»Q-[u(x)-u(xO]) - ( ' ) , (5.25a) 



g2 /TQl TQl \ 
Vq = 77-712- ^ ^ ^11 + ^66 ■ (5.25b) 

Thus there is quasi-long-range order as in the absence of disorder, and 
the effect of disorder (as seen in the static structure factor) amounts to an 
increased effective temperature. 

Below the glass transition, Tg = 1 — T/Tg > 0. g flows to a certain finite 
fixed-point value and generates an unlimited increase an oc Uee oc ln(L/a) 
with increasing length scale L oc eK Therefore the roughness is larger than 
logarithmic and, to leading order in Tg, is given by (Carpentier and Le Doussal 
1997): 



Wt(x)-Wl(x) - ^T|ln2(x/a), (5.26a) 

w{a) 2 



Wt{^)-Wl{^) ~ ^r'Hx/a), (5.26b) 



Q 



( X _ ^ 2/o(a)(l + aV4)-a/i(a) 
«^ « = 6 — , 5.26c 

[2/o(a/2) - Io[a)Y 

~( N _ ^ /i(a)(l + «V4)-2^/o(a) 

w{a) = 6 ; z — , . 5.26d 

^ ' [2/o(a/2) - /o(a)]2 ^ ^ 

Precisely at the transition the parameter 

a{T,) ^ 2^ii^ (5.27) 

Cll + Cqq 

depends on the ratio of the elastic constants only. The main structural 
feature, the increase H^(x) oc ln^(a;/a), is thus similar for magnetic fields 
parallel and perpendicular to the film, cf. equation (4.32). 

A further similarity is found in the dynamical properties, where a dynam- 
ical exponent (Carpentier and Le Doussal 1997) 

^ (2 + a)[(2-«)/(2 + a)](2-")/4 

2/o(q;/2) - /o(q;) ^ ^ ^ 

(with the Euler constant 7eu) is found on the basis of the overdamped equa- 
tion of motion 

■r]QdtUa = -|^-FCa, (5.29a) 
(Ca(x,t)C/3(x',t')) = 2r/oT5,^<5(x-x')<5(t-t')- (5.29b) 
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Consequently, the vortex lattice responds with a velocity 

V oc (5_3Q-) 

to a driving force F. Since F is proportional to the current density and v is 
proportional to the electric field, the linear electric resistivity vanishes in the 
glassy phase because of 2; > 2 in agreement with the film in a parallel field, 
cf. equation (4.54). 

5.3 Disordered vortex lattice with dislocations 

The analysis of the previous subsection was based on the elastic approxi- 
mation, excluding topological defects such as dislocations. In order to check 
whether the glassy features found above actually describe the vortex system, 
it is crucial to examine the relevance of dislocations. 

By comparing the melting temperature and the glass-transition temper- 
atures, equations (5.10) and (5.23), we immediately run into a dilemma be- 
cause Tin/Tg pa [1 — (c66/cii)^]/6 < 1. At T > Tm the pure system is liquid 
and the effects of disorder cannot be studied starting from the elastic approxi- 
mation (disorder will certainly not stabilize the ordered phase) . On the other 
hand, if T < Tm, then also T <Tg and disorder generates W{x.) oc ln^(a;/a). 
Unfortunately, this increased roughness due to the divergent renormalization 
of the bond disorder makes dislocation pairs unbind, because on large scales 
the gain of pinning energy is larger than the attraction energy of dislocation 
pairs. (The situation is analogous to the XY model, where bond disorder 
beyond a finite critical strength induces vortices, see e.g. Nattermann et al. 
(1995).) 

Consequently, on length scales beyond a scale Ldisioo where disorder in- 
duces the dissociation of dislocation pairs, dislocations are relevant at all 
temperatures and the elastic description breaks down and there can be nei- 
ther a true melting transition nor a true glass transition. Only on length 
scales below Ldisioc the system can be described within the elastic approxi- 
mation. This scale is of the order of the vortex spacing slightly above 
but it can become very large below if disorder is very weak. In this case 
melting survives as a crossover phenomenon, i.e., as a rapid decrease of Ldisioc 
if the temperature is increased above a renormalized T^* < Tm. 

The length scale Ldisioc can be estimated on the basis of the above RG 
flow equations (Carpentier and Le Doussal 1998). For qualitative purposes, 
one can examine Ldisioc from the simpler random-field XY model. Using 
further approximations Giamarchi and Le Doussal (1995) and Le Doussal and 
Giamarchi (1998a) found that for temperatures shghtly below topological 
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defects appear beyond 



W~i?ae-[(^-/^-^)^'^(^"/")]^'^ (5.31) 

Here Ra is the length scale on which the displacement becomes of the order 
of the vortex spacing, W{Ra) = a^, and Ci is a numerical constant. The last 
equation overestimates -Ldisioc at low temperatures, where Ldisioc saturates 
at a finite value for sufficiently weak disorder. This value is obtained from 
equation (5.31) after replacing T by an effective value T* < that depends 
only on the disorder strength (Carpentier and Le Doussal 1998, Le Doussal 
and Giamarchi 1998a). 

Given the finiteness of Ldisioc even at very low temperatures, one has to 
draw the conclusion that in a superconducting film with perpendicular field 
an elastic vortex glass, where dislocations would be irrelevant on asymp- 
totically large scales, does not exist. Evidence for this fact stems also from 
numerical calculations (Shi and Berlinski 1991, Gingras and Huse 1996, Mid- 
dleton 1998, Zeng et al. 1999a). The resulting physical situation can be 
thought of as a vortex lattice broken into crystallites of size -/^disioc an 
immediate consequence of this fragmentation, the system has only short- 
ranged positional order (but the correlation length diverges with vanishing 
disorder strength for T < Tm). A finite correlation length also implies finite 
pinning energies for these crystallites. Therefore thermal activation leads to 
a linear current- voltage relation for sufficiently weak currents. 
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6 Bulk Superconductor 



In this chapter we consider the weakly disordered bulk superconductor in an 
external field. We will first treat the system in the elastic approximation 
by neglecting dislocations. Later we show that this is indeed justified in a 
certain region of the phase diagram. 

6.1 Elastic approximation 

For not too large fields, H <^ Hc2, we may use the London approximation to 
describe the vortex- line lattice (VLL). Then the elastic Hamiltonian can be 
written in Fourier space as 

^ei ^Ij^^ (^^t"^^^^' + C44gJ]|ui(q)|2 + [ceegi + C44gJ]|uT(q)P}, (6.1) 

involving only wave vectors q = (qy, q^) within the first Brillouin zone (BZ). 
The displacement field u = + ut has been decomposed into its longitu- 
dinal and transverse component with respect to q_L as in equation (5.7). We 
extended the model from three to general dimension d since we will perform 
below an expansion around d = 4 dimensions below. In this generalization 
q± stands for the two-dimensional component of the wave vector perpendicu- 
lar to the generalized magnetic field and qy stands for the (d — 2)-dimensional 
parallel component. 

In principle, the elastic constants are scale (momentum) dependent. How- 
ever, since we are interested in the case of weak disorder, which becomes 
relevant only on large scales (L ^ A), we may neglect the dispersion of the 
elastic constants. In the case of strong disorder, in which relevant length 
scales (like the Larkin length) are smaller than A, the consideration of the 
dispersion of the elastic constants can be achieved in principle (see Kierfeld 
(1998)). 

Next we include the interaction between the randomly distributed impu- 
rities and the vortex lattice. As in the previous chapters for simplicity we 
assume that the impurities result in a random potential V{r), which has zero 
average and is Gaussian correlated with two-point correlations 

V(r)V(0) = A(x)<5(^-'Hz). (6.2) 

The weight of this correlator is Aq = / d^xA(x) = fpin'n'f^^^, where /pin de- 
notes the pinning force of an individual impurity, n'f''' is the impurity density, 
and ^ the maximum of the coherence and disorder correlation length (Blatter 



82 



et al. 1994). The pinning energy is then 



= J2 [ d'^-^zV{X + u{X,z),z)= f d'^rV{u{r),r) (6.3) 
^ J J 

where Pu(r) = X^x ^^^H^^ ^ ^(-^5 ^)) denotes the vortex-hne density (per 
unit area) and X is a vector of the Abrikosov triangular lattice. V{u{r), r) = 
Pu(r)V^(r) is the pinning energy density (per unit volume). 

It is interesting to note that the energy^ of the distorted VLL has certain 
symmetries. The pinning energy density V{u,r) = p-a{r)V{r) is apparently 
invariant with respect to the transformation 

u(X, z) ^ u(X', z) + X' - X, (6.4) 

i.e., by relabeling the vortex lines. Such relabeling also leaves the elastic 
Hamiltonian invariant and this symmetry is a symmetry of the total Hamil- 
tonian 7i = Tiei + Ti-pin- 

Hei possesses a second symmetry related to a rotation of the displacement 
field and a simultaneous exchange of shear and compression modulus. To 
demonstrate this symmetry, it is useful to introduce 

Cs = V^nCm, 7 = — (6-5) 

Cll 



and to rescale z — > z' z^^/c^Jcu, which leads to 

+7^/2(VxAu)2 + (Vj|U)2}. (6.6) 

As far as the thermodynamic properties of the pure system are concerned, 
these are invariant under the transformation 7 — > 7~^, since the Hamiltonian 
is invariant under the simultaneous transformation (u = (1^1,^2)) 

7 7"\ Ui u[ = -M2, M2 — > ^2 = Ui (6.7) 

(which also amounts to transforming the longitudinal displacement compo- 
nent into the transverse one and vice versa) and the partition function is 
invariant under the rotation of coordinates. The correlation functions will 
change sign if they include an odd power of U2- We note, however, that the 
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symmetry (6.7) breaks down as soon as dislocations are present, since the 
interaction energy of a dislocation pair is not symmetric with respect to the 
transformation 7 7"^ and a full treatment including both Cqq and Cn is 
required for studying the general case. In the presence of pinning this sym- 
metry is absent even as a statistical symmetry in a strict sense. This can be 
seen from the disorder- averaged replicated pinning energy. However, we will 
find below that this symmetry is restored on large scales a. 



6.2 The Bragg glass 

Historically, three different approximations have been discussed in treating 
the pinning Hamiltonian, which turn out to be valid on (i) small, (ii) inter- 
mediate and (iii) large length scales. 

(i) In his pioneering work on pinning in type-II superconductors, Larkin 
(1970) expanded (6.3) for small displacements u 

Hpin ~ f d'^'^z J]{l^(X, z) + u(X, z) • V±V(X, z) + 0{u'')}. (6.8) 

Vj.^ (X, z) is a random pinning force per length which acts on the vortex 
line of rest position X. Although it was proved by Efetov and Larkin (1977) 
that, using perturbation theory, the lowest order (linear) term in u in the 
expansion of Ttpm (under certain conditions) already gives the exact result 
for the thermodynamic quantities, it became clear later that perturbation 
theory breaks down in random systems with many energy minima (Villain 
and Semeria 1983, Fisher 1985a). Thus we can use the expansion (6.8) only 
as long as |u| is smaller than the typical distance between energy minima of 
V{r) which is of the order of C,- This restricts the validity of the expansion 
(6.8) to length scales small compared to the Larkin length. In d = 3 one 
finds from first-order perturbation theory (Larkin 1970) 



W{r) = [u(x,^) -u(0,0)]^ 



^ + ^fy^' + + > (6-9) 

7rC44 \^Cii Cgg J 



where we introduced 



z^^'-^z^ zf = '^z\ (6.10) 
C44 C44 

The coefficient poA^'^^ = — | V^A(x)|^^jj pa /pin^f^C^/^^ has the meaning 
of the density of the fluctuations of the pinning forces in the volume fraction 
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occupied by the vortex lines. As already mentioned, this result is valid only 
as long as Vr(x, z) < which defines the Larkin lengths and parallel 
and perpendicular to the magnetic field respectively: 

r. ~ ttC^ C11C44C66 

In general dimensions perturbation theory gives a roughness exponent [cf. 
equation (3.15)] 

C = Cf ^ (6.12) 

From equation (6.9) one finds for the correlation function of translational 
order 



,S(Q,r) = e^Q-W"-)-"{o)], (6.13) 

where Q is a reciprocal lattice vector, an exponential decay with the corre- 
lation length (a|^/47r^^^)L^ and (aA/47r^^)i?^, respectively. Note, however, 
that the regime of the true exponential decay is not reached because of the 
restrictions x < R^ and z < L^. 

(ii) On length scales larger than and R^, respectively, the increase of 
W{r) will continue, but with a smaller roughness exponent since perturbation 
theory is known to overestimate the influence of disorder. 

As long as the displacement of the vortex lines is much smaller than the 
lattice spacing a of the VLL (but larger than ^), each vortex line sees its own 
random potential which cannot be reached by other vortex lines. Hence, 
since |u — u'| <^ a, 

y(X + u,z)y(X' + u',z') = A(X - X' + u - u')5(z - z') 

^ Ao5x,x'4'^(u-u')5(z-zO. (6.14) 

The pinning energy density V{u, r) fa PoV{k + u, z) then obeys the approx- 
imate relation (po = -S/^o = fl~^) 

V(u, r)V(u', r') ~ poA(u - u')5^^\yi - x')5(z - z'), (6.15) 

which agrees with that of the random-manifold model of chapter 3 
(Feigel'man et al. 1989, Bouchaud et al. 1991, Bouchaud et al. 1992). In the 
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random-manifold regime the roughness exponent Crm cannot be calculated ex- 
actly. A crude estimate is given by the Flory result Cf = (4 — d)/6 ~ 0.167 
according to equation (3.20). [Note that explicit values for exponents, which 
no longer depend on d, are given for d = 3.] Emig et al. (1999) cal- 
culated Crm from a functional renormalization-group treatment and found 
Crm ~ 0.175 which varies only weakly with 7 = cee/cn (see below). The 
decay of S'(Q,r) ^ exp[— Q^l^(r)/2] in this regime is therefore of stretched 
exponential form. 

(iii) Finally, for L :^ La ^ L^^a/ifl^'^^ (or R-> i?^(a/0^/'^™) the 

vortex line displacement becomes of order a. Clearly, in the approximation 
for the perturbative regime and the manifold regime used so far, l/(u, r) 
does not fulfill the invariance property (6.4). However, it is precisely this 
property which determines the physics in the regime of very large length 
scales L ^ La {R ^ Ra)- In this regime the effect of disorder is very weak 
since displacements of vortex lines larger than a are not very favourable 
because there is already a vortex line within a distance a of any impurity. 
To keep this feature in our pinning Hamiltonian, we rewrite the vortex-line 
density as 

p„(r) = J]5(x-X-u(X,z)) = /"rfV^(5(X-x')5(x-x'-u(x',z)) 

X X 

= f d^x' 5(x - x' - u(x', z))po J] e'"^ "' 
^ Q 

^ Pop— ^ ^^g.Q.[x-u(r)]_ (g_lg) 

I ^ ~r 0(xUa I Q 

The pair correlator of the pinning energy density V{u, r) = pu(r)V(r), which 
is the only quantity which enters the following calculation, is then 



y(u,r)V^(u',rO pgA(x - xOc^(z - zO 1 h , 1 

X ^ giQ-(x-x')+i(Q+Q')-x'-iQ-u-jQ'-u'_ (6.17) 

Q,Q' 

To exploit (6.17) further, a few remarks are in order: 

(i) The terms in the sum over Q' with Q -|- Q' on the right hand side 
of (6.17) are rapidly oscillating on scales |Ax'| a and therefore average to 
zero. (For weak disorder, where L^ ^ a, this is even true in the perturbative 
regime.) 

(ii) The denominators lead to terms of the form (T{daU'^){dpn!'^) in the 
replica Hamiltonian (density). In d > 2 dimensions ^ = {2 — d)a + . . . such 
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that these terms renormahze to zero for weak disorder. We therefore omit 
them in the following, 
(iii) We approximate 



V{u, r)V{u', r') ~ A(u - u')5(r - r'). (6.18) 

Hence we get 

A(u) =pIJ2 A(Q)e^^- = po E + X)' (6-19) 

Q X 

where A(Q) is the Fourier transform of A(x), 

A(Q) = J d^x A(x) e-*Q ". (6.20) 

Apparently, the correlator A(u) is invariant under the transformation u — > 

u + X, where X is an arbitrary vector of the Abrikosov lattice. If A(x) ~ 
(Ao/27r^2) gxp(_x2/2^2) ^Yien A(Q) = Aq exp(-Q2^V2) and 



A(^) = -i ViA(u) 



u=o 2 

Q 



J2 A(Q)Q' ~ PoAo/27re^ (6.21) 



To obtain the disorder averaged configuration of the VLL on a particu- 
lar length scale, one has to take into account the rcnormalization of A(u) 
by fluctuations on shorter length scales. This can be done systematically 
by a functional renormalization group (FRG) (Fisher 1986a) for the rephca 
Hamiltonian Hn resulting from equations (2.16), (6.3), and (6.18). After the 
disorder averaging we obtain 

Hn i$^y"(iW-2;2{cn(Vx-u'^)' + C66(VxAu")2 + C44(V||U")2} 

d'^r A(u" -u''). (6.22) 



--E/ 

ab 



Because of the statistical invariance of the replica Hamiltonian under a shift 
of u by an arbitrary vector field inducing a compression, shear and/or tilt 
of the VLL, there is no renormahzation of the elastic moduh (Hwa and 
Fisher 1994b). Therefore the temperature obeys the exact fiow equation 
dT/dl = —{d — 2)T leading to a T = fixed point for d > 2. Notice, how- 
ever, that in the original model, (6.1) with (6.3), the statistical invariance 
is not fulfilled exactly on length scales smaller than a, leading to a small 
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renormalization cu — > ca of the elastic constants which will be considered 
from now on as effective parameters (similar effects have been considered in 
some detail in section 4.1). This renormalization is negligible for the calcu- 
lation of the leading asymptotic behaviour on largest scales. In the above 
approximation, Hn including the pinning energy is invariant under the trans- 
formation (6.7). Therefore this transformation will be present as symmetry 
in all thermodynamic quantities (such as the the displacement correlation) 
calculated from (6.22). 

The FRG was applied to (6.22) in the case of a sca/ar displacement field 
u first by Giamarchi and Le Doussal (1994, 1995). Here we follow closely the 
derivation of Emig et al. (1999) who treat the general case of a vector field 
u. In the present FRG only coordinates are rescaled, r exp(dl)r, to keep 
the cutoff A fixed with Adl the infinitesimal width of the momentum shell. 
Because of the dispersion of the elastic constants on scales smaller than the 
penetration depth A, we have to choose here A ?a 2n/X. Fluctuations on 
smaller scales can be ignored if the Larkin length is much larger than A, 
i.e., for weak disorder. For larger disorder one has to take into account the 
dispersion of Cn and £44, which will result in a more complicated cross-over 
but which will not affect the asymptotic behaviour (see also section IV. D in 
Giamarchi and Le Doussal (1995)). The flow equation for A(u) can then be 
derived along the lines discussed in detail in Refs. (Fisher 1986b, Balents 
and Fisher 1993). In contrast to previous cases, Emig et al. (1999) took into 
account the existence of a longitudinal and a transverse part in the elastic 
propagator. With the replacement 

^ A(u) ^ A(u) (6.23a) 



2a2 



C 



I { 

J|q|=A [ 



1 



{cuqI + C44gJ)^ {cmql + C44gJ)^ 
1 1 + 7 



Stt^ C44C66 



A-^ (6.23b) 



one obtains to lowest order in e = 4 — d a flow equation for the renormalized 
and rescaled correlator on scale L = Ae' (as in most of the previous chapters 
we suppress here the RG-flow variable /, i.e., we express A(u, L) = Ai{u) as 
A(u)) 



dl 



^^^^^ = eA(n) + f[{d!A)\{d!A)\2{d,d.Ay 

(dlA-d^Ay + 4:(did2A 



+2A(2) [dl + di) A - ^ 



(6.24) 
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with the dimensionless parameter 
A(2) = -ai^iAfu) 



u=0 



aa^sAlu) (6.25) 



u=0 



= 

u=0 



and di = d/dui etc.. Both the last equahty as weU as (5i(52A(u 

foUow from the requirement of hexagonal symmetry for A(u). The anisotropy 
parameter is given by 

5{^) = 1 - ^MlL = 1 - I, /" , ^ 2 ~ (6-26) 

7-7" Jq (ciigl + C44gj) (ceegi + C44gj) 

i.e., < 5 < 1 for any ratio 7 = cee/cn. 

In the special case Cn = Cgg (i-e., 5 = 0) the flow equation (6.24) reduces 
to that of Balents and Fisher (1993). If Cn 00, as often assumed for VLLs, 
5 = 1. 

To obtain the renormalized function A(u) on large length scales L includ- 
ing the fixed point Aoo(u) for L 00, one has to integrate equation (6.24). 
With the bare correlator of equation (6.19) showing the full symmetry of the 
triangular lattice - translation, sixfold rotational axis, six mirror lines - and 
the flow of equation (6.24) preserving these symmetries as it ought to, the set 
of possible solutions is restricted to functions with the full lattice symmetry 
on every length scale. One can solve (6.24) in a straightforward manner by 
rewriting the functional flow equation as a set of non-linear ordinary flow 
equations for the Fourier coefficients A(Q), cf. equation (6.20). Rather than 
solving for the fixed point Aoo(Q) directly, Emig et al. (1999) numerically 
integrated the flow equations exploiting the remaining point group symme- 
tries. It turns out that there is convergence to a flxed point (the pinning force 
correlator of which is illustrated in flgure 8) from a large basin of attraction. 
Of special interest is the flow of AJ , since the renormalized propagator is 
proportional to 



where 



^4 (6-27) 



AS(g-^) = A[l(^/^)(A/g)^-^ (6.28) 



which determines the roughness exponent C, of the VLL. As shown in flgure 



9, AL, , exhibits three scaling regions. 



(i) On scales L = q ^ < L^^ Ain(AL) ^ k^^Y which reproduces the 
result of Larkin (1970) for the perturbative regime. 
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Figure 8: Illustration of the pinning-force correlator — V5^Aoo(u) at the 

RG fixed point after Emig et al. (1999). The cusp-like non-analyticity 
at u = is a characteristic of the Bragg-glass phase. The hexagon 
represents the Wigner-Seitz cell of the vortex lattice. 



(ii) In the region < L < La, Af^AL) ~ {ALf'^^'^^''\ where Crm(7) is 
the roughness exponent of the random- manifold regime. Numerically, Crm(7) 
ranges from 0.1737 for 7 = to 0.1763 for 7 = 1 and is continuously increas- 
ing in this interval (see figure 10). 

(iii) Finally, for La < L, A^'^^ approaches a fixed point which de- 

~ (2) 

termines the asymptotic Bragg glass regime. The numerical value of Aoo , 
which is of order e — A — d ^ 1, still depends on 7 = cee/cn. 

With the numerical value for A^^ at hands, the exphcit form of the 
displacement correlations 



in the Bragg-glass phase is given by (Emig et al. 1999) 



(6.29) 



W^ii(r) = 



AS)(7)a 



A 



1 + 7 



In 



+ 7ln 

2 



l-7-4lnCl + 5V7ilnCl + 5 



2/\^^{-f)a\xrX2 f 



X2 \ 



7 



l-74lnfl + ^ 



x^ 



x^ V 



■(6.30a) 



(6.30b) 
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Bragg glass 
regime 



I = ln{L/X) 



Figure 9: RG flow of Ap-* through three regimes (Emig et al. 1999): 
The random- force regime with a roughness exponent d-f = (4 — d)/2 for 
length scales L < L^, the random- manifold regime with a non-universal 
exponent Crm ~ (4 — (i)/6 for length scales < L << La, and the 
asymptotic Bragg-glass regime with logarithmic roughness is reached 
for L > La- For weak disorder and ^ <C a, the width of the crossover 
regions is small compared to the width of the regimes. 




0.174 



1.156 



- 1.152 



1.148 



1.144 



Y 

Figure 10: Non-universal variation of the roughness exponents i^i-m of 
the random-manifold regime and the correlation function exponent fjQ^ 
with anisotropy parameter 7 = cqq/ch (Emig et al. 1999). 
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and ^22(1") follows from VFii(r) by permuting cn and cee- To lowest order 
in e = 4 — these correlations lead to the translational order correlation 
function 



S{Q,r) 



(exp(-^Q-[u(r)-u(0)])) 



gciLy{x'^ + Zt) 2(1+7) + 2;,^) 2(1+1/^) 



(6.31) 



with the non-universal 7-dependent exponent 



(6.32) 



and the geometrical prefactor 



qq = exp 



A£^(7)(aAg)' r 




1 




(6.33) 



which completely describes the angular dependencies of the translational or- 
der. Note that the factor qq goes to 1 in the limit z ^ 00. Therefore, in this 
limit the dependence of S'(Q, r) on the reciprocal lattice vector Q remains 
only in the exponent ffq. Moreover, it is interesting to note that the expo- 
nents of the algebraic decay in equation (6.31) depend on the elastic moduli 
as soon as z is finite even without taking into account the non- universality 
of the exponent f/q itself. If one ignores the non-trivial 7 dependence of f/q, 
in the case z — the above formulas reduce to those found in Giamarchi 
and Le Doussal (1994, 1995). A logarithmic roughness of l^(r) was already 
predicted earlier from scaling arguments (Nattermann 1990) and found also 
from a variational treatment with replica-symmetry breaking by Korshunov 
(1993) and Giamarchi and Le Doussal (1994, 1995). However, this method 
is not able to capture the non- universality of fjQ. 

The 7-dependence of the exponent fjq^ for a smallest reciprocal lattice 
vector Qa is depicted in figure 10 and varies numerically from fjQ^ — 1.143 
for 7 = to fjQ^ = 1.159 for 7 = 1. In isotropic superconductors at low 
temperatures, where vortex lines interact via central forces, one has < 
7 < 1/3. 7 ~ 1/3 for A < a, i.e., for fields close to ifci, and 7 — > 
for H Hc2- For most of the field region 7 ^o/IQttX'^B. Thus, an 
increase of the external field from Hd to Hc2 should result in an increase 
of fjQ and a decrease of Crm- At higher temperatures, where the vortex line 
interaction is rcnormalized considerably by thermal fiuctuations, as well as 
in anisotropic superconductors, the above inequality for 7 may no longer be 
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fulfilled. Clearly, in the latter case our starting Hamiltonian (6.22) would 
also have to be modified. 

The non-universality of fjQ could, in principle, be tested by neutron scat- 
tering experiments at changing external fields. On the experimental side, 
Bragg peaks have indeed been observed in BSCCO for H < 500G by Cubitt 
et al. (1993). More recently, Kim et al. (1999) have used the decoration tech- 
nique to determine the structural properties of the vortex lattice in BSCCO 
and confirmed the existence of the perturbative and the random-manifold 
regime with (^^f ^ 0.22. On larger length scales their data exhibited non- 
equilibrium features such that the true asymptotic regime was not reached. 
So far, the resolution is, however, too weak to determine the 7 dependence 
of 77. Contact to the neutron scattering experiment can be made by the 
structure factor 

dz e^O^-^^^'^''^^^e'~^'^"''^^"|^^'^'^\ (6.34) 



where we have used the Gaussian approximation (which is correct to order e) 
for the distribution of u and the definition (6.29) for Wap- With = Q+q± 
and \q±\ -C |Q| we get 

S{Q + qx, h) ^poj £r e'('i^-"+'=^")5(Q, r). (6.35) 

5'(k) describes the divergence of the scattered intensity with vanishing q in 
the vicinity of a reciprocal lattice vector Q. The above integral is dominated 
for small q by the large scale S'(Q,r), provided fjQ{'~f) < 3. It is thus the 
sub-algebraic growth of the displacements (6.30) that gives rise to the Bragg 
peaks, hence the name 'Bragg glass'. In the special cases 7 = and 7 = 1 
one obtains 

/ r \ [-3+*?q(t)]/2 
^(Q + q±,M~ + — . (6.36) 

To summarize the situation in impure bulk superconductors: it turns out 
that these still show a quasi-long-range ordered 'Bragg-glass' phase which is 
described by a non-universal power-like decay of the order parameter correla- 
tions. In particular, the decay-exponent fjQ depends on the ratio 7 = cee/cn 
of the elastic constants, similar to 2D pure crystals at their melting tempera- 
ture. For weak disorder we find a crossover of the structural correlation func- 
tions S'(Q, r) from a Larkin regime, where perturbation theory applies and 
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S{Q, r) decays exponentially, to the random-manifold regime with a stretched 
exponential decay of 5'(Q,r) and eventually to the asymptotic Bragg-glass 
regime. 

In addition to the disorder-averaged positional correlation function 
5'(Q, r) it is interesting to also consider the glass correlation functions 
'S'pg(Q, r) and Cvg(Q, r). In the framework of the functional RG ind — A — e 
dimensions, discussed above, it is easy to show that to first order e one finds 
(Bogner et al. 2000) 



'PG 



m>l 



2m 



(6.37) 



where 4"(5(Q,r) = [^^(Q.r)]^. 5W(Q,r) = exp[-iQ«Q^l^th,a/3(r)] de- 
notes the pair-correlation function for TLRO of the pure system, which is 
finite for r — > cxD if d > 2 and is reduced with respect to unity merely by 
the standard Debye- Waller factor. For large distances |r| the leading correc- 
tion in (6.37) comes from the term with m — 1. Since Ci > 0, the disorder 
increases slightly the glass order with respect to the pure system as has to 
be expected. With S'pg(Q, r)/[5'(Q, r)]^ ^ oo for r ^ oo, we conclude that 
the Bragg-glass regime is indeed a positional glass. Using the relation (2.25) 
it is also easy to show to order e that CvG(r) decays exponentially in d < 4 
dimensions (and algebraically for d — A). Expanding all corrections from the 
disorder with respect to e = 4 — d we get 



<:^vG(r: 



~ ^-\m(r)-sm?)t^ 



-) +0 



, (6.38) 



where 



((50(r) - bm)\^ 



,'27r 



— cos(q ■ x)] 
q ^''(ceegi + C44g|) 



1, , 

~ - x 

e 



(6.39) 



denotes the phase fluctuations of the order parameter in the pure system (see 
Moore (1992)). Although the leading correction (m = 0) from the disorder 
increases CvG(r), it does not compensate its exponential decay originating 
from strong thermal fluctuations in the pure system. Whether this result 
remains qualitatively correct if higher order terms in e are taken into account 
remains an open question. For the moment we conclude that to order e there 
is no phase-coherent vortex-glass order in the Bragg-glass phase. 

Our findings are not necessarily in contradiction to the result of Dorsey 
et al. (1992) who, starting from a random-Tc GL model and using mean- 
field theory, found a transition to a phase with phase-coherent glass order. 
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However, this calculation - once fluctuations are taken into account - proves 
the existence of a vortex glass transition with a diverging susceptibility Xvg = 
J df^r CvG(r) only in c? = 6 — e dimensions (e ^ 1). In this region {d > 4) 
our calculation also gives both 5'pG(Q,r) and Cvg(i") non-zero for r — > oo. 

6.3 Stability of Bragg glass 

In this chapter we have been treating so far the vortex-line lattice in the 
elastic approximation, i.e., we have disregarded topological defects such as 
dislocations or disclinations. 

In the following we will consider the stabihty of the Bragg-glass phase 
with respect to dislocations. This problem was considered by several authors 
(Giamarchi and Le Doussal 1995, Kierfeld et al. 1997, Gingras and Huse 
1996, Carpentier et al. 1996, Erta§ and Nelson 1996, Giamarchi and Le 
Doussal 1997, Fisher 1997, Kierfeld 1998, Koshelev and Vinokur 1998). Since 
the disorder seen by the dislocation is already partially screened due to the 
elastic deformations (a fact which was overlooked in the early discussion of 
the elastic approximation by Fisher et al. (1991a)), we start this section with 
a brief discussion of the effective disorder strength Aefr(-R) acting on scale R 
in the various regimes. We will thereby ignore all numerical coefficients of 
order unity and restrict ourselves to the case d — 3. 

For simplicity, we restrict ourselves here to the case cn ^ cee such that 
we can assume V_l ■ u = 0. It is then convenient to rewrite the Hamiltonian 
using the coordinate Zt = [c^^j C/^^^l"^ z introduced in equation (6.10). The 
elastic part of the Hamiltonian is now isotropic with C44 and Cee replaced by 
their geometric mean (c44C66)^''^. The same transformation in the pinning 
part changes A(u) into (c44/c66)^/^A(u). It is also convenient to introduce 
the (anisotropic) distance 

R\ = ^^ + zl (6.40) 

^csiR) introduced already in equation (6.28). We use here (6.27) to 
find explicit expressions for the mean square displacement u^lRt). To obtain 
an estimate for u^{Rt) ^ W{Rt)/2 on the scale Rt one equates the elastic 
energy E^i ~ (£44065)^^^ u^-Rt, with the fiuctuation of the pinning energy 

'^.-'■'-'(^f'^'^i^'SW)^)"'- (6-41) 

^ (2) 

Agg (i?t) denotes the curvature of the effective potential correlator at u = 
on the length scale Rt- The estimate for the pinning energy corresponds to 
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the use of perturbation theory apphed on the renormalized random potential. 
In this way we obtain 

^ (2) 

^2 ^ ^-1/2 -3/2 T (2). p V p _ ^2 ^eff W ^ (p,Ar,\ 
U ~ C44 Cgg ^eff - ^ iD-4^j 

where A^^^ = A^'^\R^) denotes the bare value of A^'^\Rt). 

Now we consider the different regimes already discussed in section 6.2. 

(i) In the perturbative regime i?t < i?j pa ^"^01^0%^ / A^'^\ A^g pa A^^^ in 
equation (6.42). 

(ii) In the manifold region, < Rt < Ra, where ^'^(Rt/R^y^^'^, we 
find from equation (6.42) 

(iii) Finally, in the asymptotic Bragg-glass region, < Rt, we get from 
^^aHn{Rt/Ra) 

p 2 2 1/2 3/2 

A2(i?.)-A(^^ff,-^^^- (6.44) 

Thus, Ay(i?t) is independent of the bare disorder strength. Here we have 
to ignore the logarithm since it originates from the summation over many 
different length scales. 

We use now these expressions to obtain an estimate for the energy of a 
dislocation loop. For simplicity, we assume an isotropic (i.e., almost circular) 
loop of linear size Rt- In the pure system, its elastic energy is 

E'^"" « 6'(c44C66)'/'i?t HRt/ai), (6.45) 

where |b| = and b denotes the Burgers vector of the dislocation. The 
logarithm results from the fact that in a distance Rt from the dislocation 
centre the displacement Udisioc produced by the dislocation line obeys |V A 
UdisiocI ~ l!^' Integration over the plane perpendicular to the dislocation 
line then yields ~ a^(c44C66)^^^ ln(i?t/ai) for the energy per unit length, 
where oi acts as a small distance cut-off. 

In a random system, the situation is more complicated. Since the dis- 
placement Udisioc created by a dislocation line obeys the saddle-point equation 
6T-C/6u = as well as f du = b, the elastic energy of the dislocation now de- 
pends also on the disorder. Since the situation apparently is rather involved, 
we will adopt here a simplified picture using the following observation: The 



96 



displacement in the neighbourhood of a dislocation has two sources: Its very 
existence leads in the distance Rt from the dislocation to a (essentially tan- 
gential) displacement of order a. The disorder, on the other hand, creates 
displacements (both in the absence or in the presence of a dislocation) of 
order a or larger only on scales Rf > Ra- In calculating the energy of a loop 
of linear size Rt > Ra in the disordered case we estimate the elastic energy 
of the dislocation loop therefore by 

^ a^cuCeey/^Rt HRJai), (6.46) 

since distortions originating from the random potential are small compared 
to a on scales i? < i?a but are dominating on larger scales. So far we have 
ignored the dispersion of the elastic constant C44, which becomes relevant on 
scales small compared to A. On these scales the lattice is softer by a factor 
{Rt/\y. In order to take this effect into account, we replace the small scale 
cut-off ai in (6.46) by {aj + X^Y^^ 

The distortions on scales larger than Ra are dominated by the disorder 
and lead to an energy gain of the order (again omitting the factor \n{Rt) in 
u^, since it comes from the summation over displacements on different length 
scales) 

4r - - (^) '^'r^ (a2 - -a^(c44Cea)^/^i?* (6.47) 

where we used equations (6.41) and (6.44). The total loop energy of a dislo- 
cation loop can therefore be written as 

« a\c^ce,y/'Rt (in (^2^) " ^i) , (6-48) 

where Ci is a constant of order unity. (It is easy to see by analogous ar- 
guments, but using equations (6.42) or (6.43) for Aef[{Rt), that dislocation 
loops of size Rt < Ra always have a positive energy.) 

From (6.48) we conclude that the system is stable against the formation 
of a dislocation loop as long as 

Ra>C2{al + \'Y/' (6.49) 

with C2 ~ e'^^. As shown by Kierfeld (1998) for strong disorder (i.e., R^ < a), 
this relation can be rewritten in the form of a Lindemann criterion. Indeed, 
in the manifold regime (note that the Larkin regime vanishes for strong 
disorder), where 

^(«»)~^ (6.50) 



u2(a) V«' + A' 
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and u^{Ra) ~ one obtains from (6.49) 

([u(a)-u(0)]2)'/' < cu a, (6.51) 

which is the Lindemann criterion appropriate for disordered systems in which 
the mean-square displacement of a single vortex line diverges, cu = c^^'^™ 
denotes the Lindemann number. The main conclusion from this consideration 
is that the clastic, dislocation-free Bragg glass is stable as long as the criterion 
(6.49) or, for strong disorder, equation (6.51)] is fulfilled. 
A few concluding remarks are in order: 

(i) The criterion (6.51) [or, equivalently, equation (6.49)] was also derived 
by Kierfeld et al. (1997) and Carpentier et al. (1996) via a variational 
treatment for a layered system with the magnetic field parallel to the 
layers. In this model only dislocations with Burgers vectors parallel to 
the layers are allowed. The most important result is the determination 
of the Lindemann number cli ~ 0.14. We refer the reader to these 
papers for further details. 

(ii) Erta§ and Nelson (1996) took equation (6.51) as a starting point for 
a discussion of the onset of irreversibility and entanglement of vortex 
lines in high-Tc materials. 

(iii) Kierfeld (1998), on the basis of equation (6.51), and including the dis- 
persion of the elastic constants in detail, calculated the stability bound- 
aries of the Bragg-glass phase for YBCO and BSCCO. The resulting 
phase diagram is depicted schematically in figure 11. Since the shear 
modulus decreases for large and small fields, respectively, there are two 
corresponding stability boundaries. 

(iv) Fisher (1997) undertook a much more detailed study of the stability of 
a defect loop in the random-field XY model which corresponds to the 
layered model mentioned in (i). Although the details of his energy esti- 
mate for the dislocation loop are slightly different from those presented 
here (using a statistical tilt symmetry, he estimates ^''""^ by equation 
(6.45) but includes also rare fiuctuations in the energy gain from the 
disorder, which also include logarithmic corrections), his final conclu- 
sions concerning the stability of the defect-free phase are essentially 
identical to those presented here. 

(v) Ryu et al. (1996b) numerically investigated the field-driven transition 
from a dislocation-free to a dislocation-dominated phase in an impure 
layered superconductor and found good quantitative agreement with 
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Figure 11: Schematic phase diagram of the vortex system according 
to Kierfeld (1998) and Kierfeld and Vinokur (1999). In the absence of 
disorder the vortex lattice melts into the vortex liquid at low and high 
magnetic fields (dashed line). Due to the presence of disorder the vortex 
lattice becomes the Bragg glass in a reduced stability region (shaded 
area), where the system is stable to the proliferation of dislocations. At 
low temperatures and low or high fields, disorder induces the prolifer- 
ation of dislocations for energetic reasons and the phase becomes the 
vortex glass. The phase transition from the vortex lattice to the vortex 
liquid has a critical end point. 



the estimates made above as well as with the neutron-scattering data 
of Cubitt et al. (1993), who found the disappearance of Bragg peaks 
in BSCCO for H >500G. 

(vi) Most recently, Kierfeld and Vinokur (1999) used a random-stress model 
of the form (3.33) as an effective model to describe dislocation lines in 
an impure superconductor. The correlations of the random-stress field 
/x(r) on different length scales are chosen in such a way that the correct 
roughness exponents in the manifold and the Bragg-glass regime, re- 
spectively, are reproduced. Considering in particular dislocation lines 
with main orientation parallel to the magnetic field, which enter and 
leave the sample at its surface, Kierfeld and Vinokur (1999) found 
a global phase diagram, which includes - besides the elastic Bragg- 
glass phase with vanishing dislocation density - an amorphous vortex- 
glass and a vortex-liquid phase separated by a first-order transition line 
which terminates in a critical point (see figure 11). 
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7 Vortices driven far from equilibrium 



As the main result of the previous chapter we found that weak pinning re- 
duces the positional order of the vortex lattice in a bulk superconductor from 
long-range order to quasi-long-range order. The discussion was restricted to 
thermodynamic equilibrium. However, since the electric resistivity is one of 
the most relevant physical properties of superconductors, it is important to 
study the driven n on- equilibrium situation. Since the efficiency of pinning is 
related to the structure of the vortex system, it is of particular interest to 
characterize this structure. In addition, since the melting transition of the 
vortex lattice shows up as a pronounced shoulder in the transport charac- 
teristic (see, e.g., Safar et al. (1992), Kwok et al. (1992)), it is desirable to 
locate this transition in the non-equilibrium situation. 

The VLL shows the three regimes of creep, depinning, and flow in its 
transport characteristic, which were already described for manifolds in chap- 
ter 3.2 (cf. figure 5). We recall that in a superconductor the vortex velocity 
is proportional to the electric field and the driving force is proportional to the 
electric current density. So far, we have mainly focused on the creep regime 
when we addressed transport properties in the previous chapters. In this 
regime, which is close to equilibrium, the dynamical response is determined 
within the scaling approach by the structure of the VLL in equilibrium. In 
particular, the logarithmic roughness of the elastic vortex glass (described 
by a roughness exponent C — 0) results in a drift velocity 



where the effective barrier height U{F) scales with the driving force F ac- 
cording to 



with the creep exponent (Nattermann 1990) 



[cf. equations (3.71) to (3.73) with the appropriate substitution for the di- 
mension of the manifold]. In d < 2 one thus expects /x < 0, which means that 
the system has a linear transport characteristic at small driving forces. Only 
in (i > 2 we can expect to find true superconductor with a vanishing linear 
resistivity. The {d = 2)-dimensional case is marginal since it has the creep 
exponent n = 0. There the effective barrier height depends logarithmically 
on the driving force and results in a power-law transport characteristic as 



v(F) ~ (F/rjo) e-^(^)/^. 



(7.1) 



U{F) ~ F-" 



(7.2) 
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found in equations (4.54) and (5.30) for the superconducting film in parallel 
and perpendicular field. 

In this chapter we will focus on the flow regime, which is far from equilib- 
rium but accessible to a theoretical analysis because it is, roughly speaking, 
close to the pure case since disorder is dynamically 'averaged out'. We will 
elucidate this point later on. Although the flow regime might appear to be 
of no special interest at first sight, only recently it turned out to be quite 
non-trivial, including phenomena such as non-equilibrium phase transitions. 

To perform the theoretical analysis of a stationary state at average veloc- 
ity V it is convenient to define the displacement u(X, z, t) = x(X, z, t)— X— vt 
by subtracting the average temporal displacement, such that u still can be 
considered as a small quantity. X is then the ideal vortex position measured 
in a frame moving with velocity v. In analogy to manifolds, the equation of 
motion for the vortex lattice, 

rjou = F''^* + FP^'^ + F - r^ov + C, (7.4) 

is over-damped. In comparison to equation (3.53) for a single vortex line, 
equation (7.4) includes the vortex density po as additional factor, such that 
the latter equation is an equation for the force density in the (i-dimensional 
space. Accordingly, the left-hand side is the Bardeen-Stephen friction force 
?7ou (Bardeen and Stephen 1965), with the friction coefficient r)o ~ BHc2l PnC^ 
referring to the force per unit volume, whereas the previous expression (3.50) 
refers to the force per unit length of a vortex line. F"^' is the force acting on 
a vortex due to its interaction with other vortices. In a harmonic approxi- 
mation for the elastic phase one has F'"^* ^ F^'. The instantaneous elastic 
force F^' depends linearly on u and reads 

F^'(q,i) = r,,;3(q)^^(q,^), (7.5a) 
p 

where F represents the elastic dispersion, p — L,T stands for the longitudi- 
nal/transverse polarization, Cp = Cii,Cqq respectively, and = q'^ + q^- As 
far as we consider a general dimension d, x and q_L are two-dimensional and 
z and qy are (d — 2)-dimensional. Fp™ is the pinning force density 

FP"XX, z, t) = -poV^ViX + vt + u(X, z, t),z), (7.6) 

which implicitly depends on the displacement u. The pinning potential V and 
the thermal noise C are assumed to be Gaussian distributed with correlators 
(6.2) and 

(Ca(X, z, t)0(0, 0, 0)) = 2rjoTpo 5x,o 5apS{z)5{t) (7.7) 
in analogy to equation (3.52). 
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7.1 Qualitative aspects 

On the basis of this equation of motion we can become more specific in what 
sense disorder is 'dynamically averaged out' in the limit of large drift veloc- 
ities. In this limit the pinning force FP'"^(X, z, t) acting on a fixed vortex 
element at position (X, z) in the comoving frame changes rapidly as a func- 
tion of time. Since the VLL has finite response times on finite length scales, 
the displacements induced by the pinning force decrease with increasing v 
and the pinning potential is effectively averaged ('washed') out. Therefore 
the effect of disorder vanishes in the limit of very large driving force. 

To substantiate this relation it is instructive to consider a single vortex 
line X moving with strictly constant velocity (i.e., x = X + vt), neglecting 
its response to pinning. In this approximation the vortex line experiences a 
time dependent force Fp^'^(X, z,t) = — poVj_y(X + vt, z). To some extent 
the effect of this force can be compared to an additional thermal noise with 
a certain 'shaking temperature' T^^. We choose the x = {x, y) coordinates 
such that V = {v,0) points into the direction of the first basis vector e^;. 
Then 

Tt = [ dtd'-^z Fr(X,z,i)Fr(X,O,0) 

^° ^ dx, dlA{xie,), (7.8a) 



2r)ov 

Tt = 0, (7.8b) 

where a — x,y is not to be summed over implicitly. These equations are 
strictly analogous to the equations (3.76) for manifolds. 

The concept of the 'shaking temperature' was introduced by Koshelev 
and Vinokur (1994), who considered in particular the case d = 2 and intro- 
duced the 'incoherent' shaking temperature T^^ = ^ "^aTt- qualitative 
purposes they considered the system driven through disorder as being sub- 
ject to a total effective temperature Tes{v) — T + T^^{v). Then the vortex 
lattice can be expected to melt at a velocity-dependent temperature 

T^{v)=T^-T'\v), (7.9) 

where Tm = T„^{v = oo) is the melting temperature of the pure system. 
Equivalently, the inverted function v^j^iT) defines a melting velocity above 
which the vortices freeze into a solid ('dynamic freezing'). T^^ vanishes 
for increasing v and hence both v„i and the corresponding driving force 
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Fm = F{vm} increase with increasing T for fixed pinning strengtfi. Prom 
the velocity dependence of T^^ one expects 



at large driving forces. This relation is consistent with experimental observa- 
tions (Bhattacharya and Higgins 1993, Hellerqvist et al. 1996) and numerical 
simulations (Koshelev and Vinokur 1994). 

A quantitative analysis of the pinning effects requires taking into account 
the interaction between the vortex lines. For this purpose Koshelev and Vi- 
nokur (1994) introduced a 'coherent' shaking temperature, which they found 
to scale like T^^hi^) ~ {vrei/v)T^^{v) with a characteristic velocity scale Vrei, 
taking into account only the transverse response of the VLL in d — 2. How- 
ever, the question of what physical properties of the system can be described 
by such an effective temperature is quite subtle. This is evident from equa- 
tion (7.8), which shows that the shaking effect is anisotropic and that one 
should distinguish the direction parallel to the driving force from perpendic- 
ular directions. 

A more careful treatment of disorder will be presented below to evaluate 
the effects of disorder on a solid VLL in an elastic approximation (section 
7.2) and to characterize the structure of the driven lattice, before one can an- 
alyze the stabihty of the VLL with respect to the prohferation of topological 
defects, from what one can locate the melting transition (section 7.3). 

7.2 Moving lattice 

To characterize the structure of the vortex lattice driven in disorder, we start 
from the assumption of an elastic and topologically ordered phase, as exists 
in the absence of disorder and for low temperatures in d > 2. The first step is 
to treat pinning perturbatively in the high-velocity regime. This analysis is 
most convenient in Fourier representation, which reads for the displacement 



where — J duj /27r and J^ — J d'^q/{2T^Y is restricted to the first Brillouin 
zone. The dynamic response function G of the pure system is (p = L, T) 



- T^{v) oc - oc - 
V r 



(7.10) 




(7.11) 




(7.12a) 



p 



G^(q,a;) 



i7]0UJ + Cpq\ + C44gy 



,21- 



1 



(7.12b) 
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To zeroth order in u the pinning force acting on the VLL in the comoving 
frame has the correlation 

E^p{X,z,t) = Fr{X,z,t)Fp{0, 0, 0) 
Q 

= p^J]AS(k)<5(u; + vk), 
Q 

(2) 

where Q is a reciprocal lattice vector (RLV) and A^^(k) = kakp A(k). The 
wave vector k = Q + q covers the whole Fourier space without restriction to 
the first Brillouin zone (in contrast to q). 

In the following we assume that the vortex lattice is oriented with one 
principal axis parallel to the average velocity. Schmid and Hauger (1973) 
argued that the vortex lattice orients itself in this direction because this is 
the direction of minimum entropy production and of least power dissipation. 
One can arrive at the same conclusion from a stability analysis (Miillers 
and Schmid 1995). In this case there are RLVs Q which are perpendicular 
to the velocity v and which play an important role for the dynamics. An 
inspection of the correlator (7.13) for u = and q = shows that the disorder 
correlator evaluated at these RLVs will determine the lattice distortions on 
large length and time scales. On these scales the system can be described 
well by approximating 

Sa/3(q,^^) ~Sc«^^(^ + v-q) (7.14) 

with 

E^f, ^pIYI ^S(Q) ^ ^'•'^S^ySfSy (7.15) 

Q(±v) 

and A^^) as given in equation (6.25). 

Prom the pinning force correlator (7.13) the zero-temperature displace- 
ment correlations follow via 

W^p{r,t) = (K(r, t) - M„(0, 0)] [Ufs{r, t) - ugiO, 0)]) 

= 2/" /"(l-e^('i-^-'^*))G,^(q,a;)S^5(q,u;)G5/3(-q,-c^)-(7.16) 

Jul J q 

Within this lowest-order perturbative approach one finds explicitly on large 
scales |r| — > 00 that the fluctuations of are finite (at T = 0), whereas Uy 



(7.13a) 
(7.13b) 
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is rough in dimensions d < 3 (Balents and Fisher 1995, Giamarchi and Le 
Doussal 1996, Balents et al. 1997): 

2 poA(0) /27r^ 



W^..(r,i=0) — ] , (7.17a) 



PF,,(r,t = 0) ^ W (po max <|^,y^^z2 1>), (7.17b) 



-.3— 

770f Ceff V L ^0^" " C44 

with an effective elastic constant Ces = cu m d — 2 and c^s = (011044)^/^ in 
d = 3 (in general dimension, Ces — (hi '^^'^c^Ja ^) ^-nd a scaling function with 
W(0) = and 

{const., d > 3, 
ln(s), rf = 3, (7.18) 
^(3-d)/2^ (i<3, 

for s — > 00. 

In comparison to the equilibrium case, the driven vortex lattice is rough- 
ened by disorder only in dimensions < 3. Thus, in (the somewhat academic) 
dimensions 3 < d < 4, the roughness disappears due to the drive, and disor- 
der is washed out substantially by driving the VLL. It is further remarkable 
that the roughness of the driven lattice (as found in lowest order perturbation 
theory) arises from compression modes. In contrast to this, in equilibrium 
even an incompressible lattice is roughened. Although the shaking temper- 
ature (7.8) suggests larger fluctuations of Uy than of in agreement with 
equation (7.17), the roughness of the lattice cannot be described by adding 
the shaking temperature (7.8) to the physical temperature of the system. As 
long as this temperature is flnite, one would expect the lattice to be flat in 
d > 2, which is inconsistent with the roughness found above in 2 < c? < 3. 

In analogy to the equilibrium case we may deflne a dynamic Larkin length, 
beyond which the perturbative approach breaks down, from W{r, t — 0) = 
Wxx{r, t — 0) + Wyy{r, t — 0) — to which the y component of the displace- 
ments gives the dominant contribution according to equation (7.17). Because 
of the anisotropy of Wyy{r, t = 0), the Larkin length strongly depends on the 
orientation of r: 



4'^ - / CmW (7.19a) 



aexp(^ """'^^f^ ), d = 3, 

^-^(Lf)\ (7.19b) 
cii \ W 

^ (7.19c) 

Cii 
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(Giamarchi and Le Doussal 1996, Balents et al. 1998a, Le Doussal and 
Giamarchi 1998b). 

The roughness of the driven vortex lattice (more precisely, of Uy) in c? < 3 
implies that the perturbation theory breaks down on large length scales and 
progress can be made only by means of a renormalization-group analysis. 
This analysis is very involved (Balents et al. 1998a, Le Doussal and Gia- 
marchi 1998b, Scheidl and Vinokur 1998b) and only its main results will be 
presented. 

One main feature is that the disorder correlator gets renormalized in 
a qualitative way. While the original pinning-force correlator A^^(k) = 
kakp A(k) entering equation (7.13) is derived from a random potential such 

that A^^(k) = 0, under renormalization it develops a random-/orce char- 

acter such that A^^^(k) ^ 0. Thereby the large-scale value of the 

pinning-force correlator will be increased to values '^oo,yy > and, most 
important, 'E.^^^x > 0. Balents et al. (1997) pointed out that a roughness 
of the displacement component will be generated as a consequence of 
the finiteness of 'E.^^xx- Thus, on large scales the displacement correlation 
function is not described by the perturbative result (7.17a) but by equation 
(7.17b) after the substitutions cn cqq and Eyy 'Eoo,xx- Since cee <C cn 
one expects the fluctuations of Uy to be larger than the fluctuations of on 
very large scales (Balents et al. 1997, Balents et al. 1998a). 

Giamarchi and Le Doussal (1996) pointed out that the displacement com- 
ponent Uy actually shows a frozen pattern in the laboratory frame. This 
means that the vortices move along 'static channels', see figure 12. These 
channels are oriented parallel to the average velocity without crossing each 
other, although they are rough in o? < 3. 
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Figure 12: Structure of the moving vortex glass (intersection at con- 
stant z). As long as the elastic glass is stable (in d > 2 for sufficiently 
high drive and low temperature), the topological order of the lattice is 
preserved (dashed lines link neighbouring vortices). Although the mov- 
ing vortex glass is rough (in d < 3 both displacement components have 
unbounded disorder-induced fluctuations) the vortices flow along chan- 
nels (shaded lines, corresponding to the time-averaged vortex density). 

Thus, in elastic approximation, the driven vortex lattice looks like a glass 
with respect to the structure, and hence one may call it 'moving glass' (Gi- 
amarchi and Le Doussal 1996). Nevertheless, unlike the Bragg glass the 
'moving glass' does not match all criteria of glassiness as specified in section 
3.3. Since Uy is quenched, one might expect to find a sub-linear transverse 
response, rjyyiy) = dFy{w) / dvy\^ = oo at finite temperature. This is, 
however, not the case: a finite transverse force Fy induces a linear transverse 
velocity Vy (Le Doussal and Giamarchi 1998b, Scheidl and Vinokur 1998b). 
Only at zero temperature there is a finite critical transverse force. But this is 
true for any potential and not indicative for the actual relevance of disorder 
for the dynamics. 

The renormalization-group analysis refiects characteristic features of non- 
equilibrium, including the generation of Kardar-Parisi- Zhang (KPZ) nonlin- 
earities, and an anisotropic renormalization of the elastic dispersion of the 
lattice, of the linear mobility and of the temperature (Balents et al. 1998a, 
Le Doussal and Giamarchi 1998b, Scheidl and Vinokur 1998b). The detailed 
discussion of these aspects is beyond the scope of this article. 
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7.3 Moving smectic and dynamic melting 

In the previous section we discussed the moving vortex system in the elastic 
approximation and found that is has long-range order in d > 3, quasi-long- 
range order in d = 3 and only short-range order in d < 3. In order to 
examine whether this elastic system is actually stable with respect to the 
generation of free dislocations, it is in principle necessary to examine their 
dynamic generation process, since in the non-equilibrium situation the lattice 
stability can no longer be examined using energetic criteria. 

Nevertheless, if we naively carry over our findings from equilibrium to 
non-equilibrium, wc might expect the elastic system (the moving lattice, 
which is topologically ordered and where vortices move coherently) in > 3 
to be stable since it has (quasi-) long-range order, whereas the short-range 
order in d < 3 should imply the instability. One can arrive at the same 
conclusion on the basis of a dynamical scaling analysis which was performed 
by Balents and Fisher (1995) for CDW systems. 

Balents et al. (1997) argued that because of the dominance of fluctua- 
tions of Ux over fluctuations of Uy on largest scales, the instability should be 
generated by dislocations with Burgers vectors parallel to v. Due to the pres- 
ence of these dislocations the channels would be decoupled dynamically, i.e. 
vortices in different channels could move with locally different velocities. An 
analogous decoupling of charge-density waves in a layered model was found 
by Vinokur and Nattermann (1997) from a variational calculation. 

Despite the dynamic decoupling it is possible that a static channel struc- 
ture still persists. This phase is called 'moving smectic' (Balents et al. 1997) 
or 'moving transverse glass' (Le Doussal and Giamarchi 1998b). The vor- 
tex density can still be modulated with quasi-long-range order in = 3. In 
(i < 3 it is possible that even this channel structure is destroyed, since the 
power-law roughness of Uy may also induce dislocations with Burgers vec- 
tors having a component perpendicular to v. Then the vortex system would 
be essentially a vortex liqTiid, which still has a certain anisotropy since the 
driving force breaks the rotation symmetry in the (x, y) plane. Although 
the moving lattice can exist at sufficiently large drift velocities in d = 3, the 
effective strength of disorder becomes larger as the velocity is reduced and 
therefore the moving lattice can decay first into a moving smectic and then 
into a moving liquid. The schematic phase diagrams (Balents et al. 1998a, 
Le Doussal and Giamarchi 1998b, Scheidl and Vinokur 1998b) for d — 2 and 
d = 3 are illustrated in figure 13. 

A theoretical analysis of the large-scale properties of the driven vortex 
system is hampered by the anisotropy of the system, the relevance of disorder 
and the dynamic non-linearities such as KPZ terms that govern the vortex 
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moving smectic ^ ' 




Tra T T 



Figure 13: Schematic illustration of the dynamic phase diagram for 
vortex systems driven in disorder, (a) In d = 2 dislocations are present 
in the vortex system for zero drive as well as for large drive. Thus at 
small drive the vortices form a liquid, for which is pinned only at T = 
and small drive F < F^. At large drive the moving vortices can start 
to follow static channels and to become a moving smectic. (b) In d = 3 
vortices form a Bragg glass below melting (T < Tm) and below dcpinning 
(F < Fc). For F > and T > this lattice will creep until the stresses in 
the lattice become very large such that the lattice breaks into a moving 
liquid. At T = this destruction of the solid presumably occurs very 
close to depinning. At larger drive, the moving liquid develops transverse 
order as in d = 2. Unlike in d = 2, the moving smectic can dynamically 
freeze into a moving lattice at even larger drive. In both dimensions 
the physics at intermediate drive (F ~ Fc) lacks a precise theoretical 
description because in this non-equilibrium regime disorder is effectively 
strong. 
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dynamics on large scales. For these reasons a rigorous analysis has not been 
achieved so far. The conspiracy of these influences may actually lead to 
further phase transitions, such as first-order roughening transitions, which 
were found by Chen et al. (1996) for charge-density-wave systems in = 1, 2. 

There is experimental evidence not only for dynamic melting, but also the 
structure of the moving lattice and moving smectic have been characterized. 
The transition was detected by resistive measurements (Bhattacharya and 
Higgins 1993, Hellerqvist et al. 1996) and the structure of the dynamic 
phases was observed by dynamic decoration techniques (Marchevsky et al. 
1997, Pardo et al. 1998, Troyanovskii et al. 1999). Numerical evidence for 
dynamic melting was found va.d — 2 (Koshelev and Vinokur 1994, Ryu et al. 
1996a, Faleski et al. 1996) and also d = 3 (Dommguez et al. 1997). 

\n d = 2 melting can show up only as a crossover which, however, may be 
very sharp for weak disorder and restricted system sizes. Spencer and Jensen 
(1997) presented numerical evidence for the absence of true topological order 
even at large velocity in. d — 2. The instability of the driven two-dimensional 
system was demonstrated analytically by Aranson et al. (1998) who explicitly 
examined the dynamics of single dislocations and their interaction within a 
driven vortex lattice. They were able to show that already the presence 
of the KPZ terms leads to a screening of the dislocation interaction on a 
finite length scale (which, however, increases exponentially with the drift 
velocity) . Consequently, dislocation-anti-dislocation pairs will unbind under 
the additional influence of thermal fluctuations and even more due to the 
shaking effect of disorder. 

In the remainder of this section we briefly readdress the melting transition 
and sketch how it can be captured by the conceptually simplest approach, 
a dynamical Lindemann criterion. Such a phenomenological approach can 
be useful in order to locate a transition semiquantitatively, but it certainly 
cannot give insight into large-scale properties. The Lindemann criterion will 
not be used for the total fluctuations of the vortex displacement, but for the 
relative displacement of neighbouring vortices. In this modified form it has 
proved to be successful in the static case even for systems which have no 
true long-range order, see equation (6.51). In order to extract information 
about the anisotropic nature of the driven state it is instructive to look at 
the relative displacement fluctuations 

^^(b) = (Kb) - u(0)]2) (7.20) 
= ^im + ^lW (7-21) 

as a function of the orientation of the vector b between the ideal position 
of the neighbours. Scheidl and Vinokur (1998a) evaluated the Lindemann 
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criterion ^^(b) < cl^a^ for the stability of the lattice and determined the 
velocity dependence of the melting transition. In this context we will not 
reproduce the detailed results. But it is interesting to point out that pinning 
results in contributions 

wl oc 



for the bond fluctuations, which imply that for large drift velocities near- 
est neighbours separated in y direction have larger bond fluctuations than 
nearest neighbours in x direction. Hence this phenomenological criterion 
indicates that the moving lattice (provided it exists at large drive) will de- 
cay into decoupled channels. This conclusion from the Lindcmann criterion, 
which examines displacement fluctuations on short scales, is in agreement 
with the analysis of fluctuations on large scales described above. Concern- 
ing the resulting shift of the dynamic melting temperature it agrees with 
the finding (7.10) from the incoherent shaking temperature (7.8c). Thus the 
different approaches, the shaking temperature as measure of the strength of 
the pinning force, the Lindemann criterion as a measure of the displacement 
fluctuations on small scales and the renormalization-group results for the 
displacement fluctuations on large scales give in combination a consistent 
picture of the physics of driven vortices. 



2 ' 



(7.22a) 
(7.22b) 
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8 Summary 



In this article we have reviewed the influence of weak pinning by point-hke 
impurities on the vortex-hne lattice in type-II superconductors. In particular, 
we have addressed the question to what extent these superconductors display 
glass-like properties and by which order parameter or correlation functions 
these properties can be identified. 

Hereby it is important to distinguish two concepts of order: (i) order in 
the position of the vortex lines, which allows for a breaking of the continuous 
translational symmetry that is reflected by a spatial modulation of the mag- 
netic induction; and (ii) the order of the superconducting condensate wave 
function, which is related to the breaking of the U{1) symmetry and which 
manifests itself in phase coherence (ODLRO). 

Pinning destroys the positional long-range order of the VLL in d < 4. 
However, for weak pinning (and for sufficiently low temperature, which is 
always assumed) quasi- long-range positional order persists in 2 < d < 4. 
This is true also in superconducting films {d = 2) in a parallel field (where 
the exponent of the algebraic decay depends only logarithmically on the 
scale), but not in a perpendicular field, where dislocations induce short- 
ranged positional order. 

In 0? > 2 the vortex lattice behaves as an elastic medium for weak disor- 
der, i.e., it is topologically ordered since there are no free dislocations. Then 
the vortex system can be called an elastic vortex glass. (A superconducting 
film in a parallel field can also be considered as an elastic glass since disloca- 
tions are excluded for geometrical reasons) . In elastic vortex glasses vortices 
are collectively pinned such that vortex motion can take place only due to 
thermal activation over arbitrarily high barriers. Such barriers are necessary 
for a vanishing linear resistivity; otherwise the superconductor in the mixed 
phase would actually be an Ohmic conductor. 

In the mixed phase vortex fluctuations are decisive not only for the degree 
of positional order in the vortex system but also for the phase coherence. In 
pure systems thermal fluctuation destroy ODLRO in d < 4 dimensions (if 
screening is taken into account), including bulk superconductors and super- 
conducting films in a perpendicular field. Preliminary calculations to lowest 
order in e for d = 4 — e dimensions suggest that this conclusion remains true 
also for systems with weak disorder. A phase-coherent vortex glass, where 
the condensate wave function is quenched, can exist only in unphysical di- 
mensions d > A. However, higher order terms in e may still change this result. 
This applies in particular to systems with strong disorder (e.g. gauge-glass 
models without screening) that substantially suppresses thermal fluctuations 
and may permit the persistence of phase coherence in d < 4. 
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When the elastic vortex glass breaks down due to the prohferation of 
dislocations, it may still possess a weaker degree of order that might disappear 
only at higher temperatures and/or stronger disorder. Since the positional 
order of the elastic vortex glass in 2 < d < 4 resembles that of the pure crystal 
in d = 2, it is possible that the elastic vortex glass breaks into a hexatic vortex 
liquid with bond-orientational order (Chudnovsky 1989, Toner 1991b) before 
it decays into an isotropic vortex liquid. A more detailed investigation of 
this possibility as well as of further more exotic disordered phases is left for 
future studies. 
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A Pinning of periodic media 



In this appendix we present the transformation of the pinning energy for vor- 
tex hne lattices, which leads to an effective periodic pinning potential. In our 
notation r = (x, z) denotes a point of space, with the D-dimensional compo- 
nent z component along the vortex lines (in general dimension, a manifold) 
and the A^-dimensional orthogonal component x of the displacement. The 
dimension of the embedding space of the VLL then is d = N + D. In the 
undistorted lattice, the vortex lines are located at positions X. The distorted 
lattice can be described by the density (lines per A^-dimensional volume) 

p„(r) = ^5(x-X'-u(X',z)). (A.l) 

X' 

The average density is denoted by po and it is related to the average vortex 
spacing a through po = a'^ ■ In particular, po — B /^q for dimension N — 2. 
The pinning energy reads 

^pin = y"d^^^V(X + u(X,z),z) =y"ci'^r pu(r)V^(r) (A.2) 
We always assume disorder to be Gaussian distributed with correlations 



V{^, z)V{x', tI) = A(x - x')5(z - z'). (A.3) 

Using the Poisson-summation formula (Nattermann el al. 1991, Gia- 
marchi and Le Doussal 1995), we may rewrite the density 

p„(r) ^ J2 f - ^' - ^)) - 

X' 

= / d^'x' (5(x - x' - u(x', z)) poJ2 e'^"' 

Q 

Q 

^ [l-a„«„(x,z)]po$^e^^-['^-"Wl 

Q 

-po<9aMa(x,z) +po^cos{Q ■ [x- u(r)]}, (A.4) 
Q 

where Q are reciprocal lattice vectors (RLV) and da = d/dxa- The displace- 
ment tells us how vortices are shifted form the "ideal" position x' to the 
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"actual" position x. In the beginning the displacement is considered as a 
function of the "ideal" position, x — x' = u(x', z). During these manipula- 
tion we have expressed the displacement as a function of the actual position, 
X — x' = u(x, z). 

The random potential V couples to the divergence of the displacement in 
the density, Eq. (A. 4), as an effective random-compression iorce. This term 
is of particular importance in dimensions d < 2. 

The remaining contributions Q 7^ to Eq. (A. 4) represent an effective 
periodic pinning potential since it is invariant under shifts u(x, z) u(x, z) + 
X with an arbitrary lattice vector X. However, since under this shift x is 
to be held fixed, the shift means x' ^ x' — X. Thus it does not mean 
a translation of the vortex lattice in the laboratory frame but it means as 
relabeling of the vortices, or a shift of the ideal reference positions x'. 

In practice it is more convenient to start with the exact pinning energy 
for a VLL in the replicated system: 

^Pin,n = /d^^^A(X + U»(X,z)-X'-U^(X',z)) 

a,b X,X' 

A(k) e^'*-[^+""(^'^)-^'-"'(^''^)l, (A.5) 

where Jy^ — J d'^k/{2T^Y Fourier transform of the correlator is 

A(k) = j d^x e-*^-"A(x). (A.6) 

The energy (A.5) can then be transformed in the following way: the wave 
vector k = Q + q is split into a reciprocal lattice vector Q and a vector q 
within the first Brillouin zone and = YIq /q' accordingly. The contribu- 
tion Q = represents the pinning energy for a continuous elastic medium. 
This contribution contains a random compression force (the terms of order q'^ 
in this contribution) which we will not discuss further here. We focus on the 
contributions Q 7^ which encode the periodicity of the VLL. For a short 
disorder correlation length ^ ^ a one may approximate A(Q + q) ~ A(Q) 
and then perform the the integration over q, after which only the contribu- 
tions X = X' survive since |u''(X, z) — u''(X', z)| <^ |X — X'| for a roughness 
exponent C < 1- After these approximations one ends up with 

a,b X Q 

a,b X 
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(A.7) 



with the effective periodic correlator 

A(u) = pgJ]A(Q)e^Q-" = poE^(^ + ")' 

Q X 

which is related to the effective potential via 

V{u, r)V{u', r') « A(u - u')5{r - r'). (A.9) 

This means that the periodic lattice effectively behaves like an elastic man- 
ifold of internal dimension Des — D + N in a. effective embedding space 
deff — D + 2N subject to a periodic pinning potential. 

The advantage of the manipulations (A.7) as compared to (A. 4) is that 
we end up with a periodic correlator for the displacement u as a function 
of the ideal reference position (which is the actual degree of freedom of the 
vortex lattice) rather than for u as a function of the actual vortex position. 
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B List of recurrent symbols 



symbol 


meaning 


definition 


a 


effective vortex spacing 




OA 


vortex spacing in triangular lattice 


= (2/v3)a^ 


A 

A 


magnetic vector potential 




B 


magnetic induction 




c 


velocity of light 




Cg, Ci, C2, ... 


numerical constants 




Cll 


compression modulus 




C44 


tilt modulus 




C66 


shear modulus 




fLi 


Lindemann number 


(2.21) 


CvG 


phase-coherent vortex-glass correlation 


(2.34) 


(i 


(total) spatial dimension 


page 24 


i:) 


(internal) spatial dimension 


page 24 


£^ 


energy 




F 


driving force 


(3.51) 




free energy 




G 


displacement response function 


(3.11) 


Gi 


Ginzburg number 


(2.5) 


H 


magnetic field 




He 


thermodynamic critical field 


page 9 


Hci 


lower critical field 


page 9 


Hc2 


upper critical field 


page 9 


Hml 


lower melting field 


page 15 


Hm2 


upper melting field 


(2.23) 


nj 

n 


Hamiltonian 




3 


current density 




J 


effective stiffness 


(4.16b) 


k 


wave vector, without restriction 


k= Q + q 


L 


system size 




La 


crossover length 


(6.2) 




screening length 


(5.3) 




Larkin length 


(3.16), (6.11) 


n 


number of replicas 




N 


number of displacement components 


page 24 


pL/T 


longitudinal / transverse pro j ector 


(5.1) 
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q 


wave vector, restricted to first Brillouin zone 




Q 


reciprocal lattice vector (KLVj 




Q 


thermodynamic quantity 


(2.17) 

Q2, = 16% /3a2, 


Qa 


length ol farst KLV m triangular lattice 




crossover length 


(6.2) 


Rt 


111 J 1 

rescaled length 


(6.40) 




T 1*1 J 1 

Larkm length 


//-> -1 -1 \ 

(6.11) 


s 


nlm thickness 




S 


translational order parameter correlation 


(2.33), (3.9) 




11 1 j_* 

positional glass correlation 


(2.33) 


s 


structure factor 


(4.35) 


s 


entropy 




T 

c 


superconducting transition temperature 






mean-field transition temperature 


page 8 


T 


glass transition temperature 


(4.25), (5.23) 




hexatic/isotropic liquid transition temperature 


(5.14) 




melting temperature 


(5.10) 


Tkt 


Kosterlitz-Thouless transition temperature 


(5.4) 


u 


vortex displacement 




U 


interaction potential 




V 


velocity 




V 


pinning potential 




w 


variance of displacement difference 


(3.6), (3.10) 


X 


component of x parallel to v 




X 


vector component perpendicular to B 




y 


component of x perpendicular to v 




z 


vector component parallel to B 




z 


dynamical exponent 


(4.50) 


Zl, Zt 


rescaled lengths 


(6.10) 


/ 


1 aiilyJ \J1. dfljOtH^ (^^JliOtCXlitO 




7eu 


Euler's constant 


7eu = 0.577- •• 


r 


elastic dispersion relation 


(7.5) 


S 


anisotropy parameter 


(6.26) 


A 


disorder correlator 


(3.3) 


Ao 


integral of disorder correlator 


(3.4) 




moment of disorder correlator 


(3.13) 


A 


periodic disorder correlator 


(A.8) 


A(2) 


moment of periodic disorder correlator 


page 87, (6.25) 
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e 


anisotropy 


page 12 


e 


dimension 


e = 4-D 


£o 


energy scale 


(5.2) 




stiffness of manifold 


(3.1) 


c 


roughness exponent 


(3.6) 


c 


thermal noise 


(3.51) 


Cf 


Flory roughness exponent 


(3.20) 


Crf 


random force roughness exponent 


(3.15) 


Crfi 


random field roughness exponent 


(3.27) 


Crm 


random manifold roughness exponent 


(3.23) 


Cth 


thermal roughness exponent 


(3.7) 


V 


correlation exponent 


(4.33) 


Vo 


unrenormalized friction coefficient 


(3.51) 


9 


scaling exponent of energy fluctuations 


(3.18) 


K 


Ginzburg-Landau parameter 




A 


magnetic penetration length 


(2.2b) 


At 


thermal length 


(2.19) 


/i 


random tilt 


(3.24) 


/i 


creep exponent 


(3.70b) 


A'rnag 


magnetic permeability 


(4.43) 




coherence length 


(2.2a) 




disorder correlation length 






pinning force correlator 


(7.13) 


P 


vortex density 


(4.10) 


Po 


average vortex density 




Pn 


normal state resistivity 


(3.50) 


a 


random tilt strength 


(4.18b) 


T 


relative temperature distance 


r = 1 - T/Tj 


$0 


flux quantum 


$0 = hcjle 




barrier scaling exponent 


(3.46) 




translational order parameter 


(3.8) 



Note that functions and their Fourier transforms have the same name 
although they are different functions. This ambiguity is removed by the 
argument of the function, which is either a length (such as x, z, r etc.) or a 
wave vector (such as k, q, Q etc.). The only exceptions are A and A. Their 
common Fourier transform is denoted by A. 
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